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SUMMARY. 


Previous authors have discussed the propagation of radio-atmospherics 
by assuming that the space between the earth’s surface and the 
ionosphere behaves like a wave-guide. If the surfaces of such a guide 
are treated as perfect conductors, then there is one mode of propagation, 
the “zero-order” mode, which is unattenuated at all frequencies. In 
such a guide the wave-form of the received e.m.f. due to one mode at a 
great distance from a lightning flash can be calculated, and for modes 
of non-zero order it would be of an oscillatory type similar to that 
observed in many atmospherics from distant sources. But it would also 
include a large undistorted impulse due to the energy propagated in the 
zero-order mode, and this is not observed in practice. Moreover, recent 
e amplitudes of the component frequencies of 
that frequencies below about 7:5 ke./s. are 
propagation from the source to the receiver, 
h the presence of a zero-order mode which is 


measurements of th 
atmospherics have shown 
heavily attenuated during 
and this is inconsistent wit 
anattenuated at all frequencies. 
This paper shows that if one s 
conductor, then all the modes, in 
ttenuated below a certain critica 


urface of the wave-guide is an imperfec\ 
cluding the zero-order mode, are heavily 
1 frequency. There are two modes 
whose amplitudes are much greater than those of all other modes, and 
hese both give an oscillatory response to a lightning flash at a great 
listance. There is no longer any mode of propagation which would give 
wn undistorted impulse at a great distance from the source. 


* Communicated by J. A. Ratcliffe. 
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The critical frequency depends most markedly on the height of the 
ionosphere, and to a smaller extent on the characteristics of the ionized 
medium. The changes in the measured curves of attenuation versus 
frequency from day to night, and during sudden ionospheric 
disturbances are consistent with the changes in the height of reflection 
of very long waves which are known to occur at these times. 

The relation between the treatment of propagation in terms of wave- 
guide modes, and in terms of successively reflected “‘ rays ” is discussed. 


§ 1. InrTRoDUCTION. 

Tue nature of the radio-atmospherics observed with a receiver at a given 
point depends on: (a) the nature of the sources; (b) the length and 
characteristics of the propagation paths. Some recent experiments 
(Gardner 1951, Bowe 1951) have been concerned with the characteristics. 
of the propagation path for different frequency components in 
the spectrum of the atmospherics. They have shown that frequencies 
below about 7-5 ke./s. are heavily attenuated. The amount of attenuation, 
and the exact frequency below which it occurs, depends upon whether 
it is day or night, and changes also during a sudden ionospheric 
disturbance. The object of this paper is to suggest a possible mechanism 
to explain these results. The mechanism suggested will also account 
for the oscillatory nature of the wave-form of atmospherics coming — 
from great distances. 
In §3, §4 and §5, it is assumed that the atmospheric travels in a 

“ wave-guide ” formed between a perfectly conducting earth and an 
imperfectly conducting, sharply bounded ionosphere. Previous authors 
have discussed the propagation of waves in this type of wave-guide, 
but their discussions have been inadequate for our present purpose for — 
various reasons. ‘Thus Hales (1948) has discussed the propagation of an 
atmospheric in a wave-guide with perfectly conducting walls. We shall 
show that, although this model at first seems attractive, the ‘‘ zero-order” 
mode is unattenuated at all frequencies, and hence the experimental. 
results cannot ,be satisfactorily explained. Eckersley (1932), in a 
mathematical paper designed primarily to illustrate the nature of his. 
““ phase-integral” method of calculation, arrived at expressions like those 
which we derive here, but he did not apply them to our problem. 
Most observers are agreed that the wave-forms of atmospherics can be 
divided roughly into two main types (Schonland et al. 1940, Laby 
et al. 1940, Rivault 1945). In one type the wave-form is short and 
irregular in shape, but the irregularities are repeated roughly at intervals 
which correspond to the repeated reflection of a pulse from the ionosphere. 
This type of atmospheric appears to be fairly well explained by the 
mechanism of successive reflection. The other type is oscillatory and the 
wave-form is much smoother than for the first type. It has been vaguely 
suggested (e.g. Schonland et al. 1940, p. 184) that it could be produced — 
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by the smooth merging together of the separate reflections of the first 
type, but no detailed explanation has been given of how this could occur 
and, indeed, the two types of wave-form are so different that it is difficult 
to see how one could be converted into the other. 

Hales (1948) came near to explaining the wave-form of an atmospheric 

from a distant source. He first discussed the explanation of the 
oscillatory type of atmospheric in terms of a wave-guide with perfectly 
conducting walls, and concluded that its wave-form could be explained 
with this model. We show in §6 that we do not agree with this 
explanation which again neglected the zero-order mode. Hales next: 
considered the case where one boundary had a finite conductivity, and 
showed that, with this model, he could explain the existence of the 
* slow tail” of the atmospheric, but he did not reconsider the production 
of the oscillatory type of atmospheric by this model. 

In § 6 it is shown how the oscillatory type of atmospheric would be 
produced by transmission in a wave-guide system of the type discussed 
in this paper, and in §7 it is indicated how the transition from the 

first to the second type of wave-form can come about. 


§2. THE NaTURE OF THE SOURCES. 


The results of Gardner and Bowe were analysed in such a way that the 
characteristics of the source were eliminated. In discussing the frequency 
dependence of the characteristics of the propagation path, therefore, 
it is unnecessary to consider the exact nature of the source of an 
atmospheric. The wave-form of a single atmospheric, however, does 
_ depend upon the disturbance producing it, and some assumption about 
the source must therefore be made. It will be assumed in what follows . 
that the atmospheric is the result of the instantaneous destruction of 
a short vertical electric dipole situated at the surface of the earth, so that 
the dipole moment of the source varies with time as shown in fig. 1 (a). 
The signal that would be received at a distance in free space from such 
a source would be a double impulse, as shown in fig. 1 (bo). In the 
continuous spectrum of this source the amplitude is proportional to the 
frequency. 


§ 3. Rays AND WAVE-GUIDE MODES. 


The propagation of an atmospheric from the initial lightning flash to 
the receiver may be discussed by either of two methods ; these will be 
called (a) the method of successive reflections, (b) the method of wave- 
guide modes. In (a) the received e.m.f. is regarded as a succession of 
impulses arriving at the receiver, the first after propagation direct from 
the source, and the succeeding ones after one, two, etc. reflections at the 
ionosphere. 

In method (b) the space between the surface of the earth and the 
ionosphere is treated as a wave-guide.. In this guide the propagation 
of certain modes is possible. The characteristics of each of these modes. 
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must be found, and also the amplitude excited in each mode by a given 
source. The received e.m.f. is then found by adding the contributions 
from all the modes. 

If the ionosphere is an imperfect reflector, method (a) is the simplest 
when the receiver is close to the source, say within a few hundred - 
kilometres, since then the amplitude of the direct wave and the first and 
possibly second reflected waves are large compared with the reflected 
waves of higher order, and these can therefore be neglected. At larger 


Fig. 1 (a). 
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Electric dipole moment 


Idealized lightning flash. Variation of dipole moment of source with time. 


Fig. 1 (0). 
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Signal received in free space from idealized lightning flash. 


distances the incidence of the waves on the ionosphere is more oblique, 
so that the reflection coefficient is higher, and consequently the received 
e.m.f. is composed of a large number of reflections of comparable 
amplitude. In this case method (b) is simpler, since it is normally found 
that there is a marked difference between the attenuation coefficients of 
the various wave-guide modes, and at large distances the least attenuated 
mode is much greater than all the others. At short distances, the ampli- 
tudes in a large number of modes are comparable, and calculation of the 
received e.m.f. by the “ mode ’* method is therefore complicated. 
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§4. HKarrn anp IoNOSPHERE CONSIDERED AS A WAVE-GUIDE. 


Let us first consider the over simplified case where the earth and the 
ionosphere are both perfect conductors. A wave-guide of rectangular 
cross-section has four surfaces. We consider the limiting case where two 
. of these surfaces are separated to an infinite distance, and we shall say 
that the guide has infinite horizontal width, but a finite vertical width, h. 
We consider only signals in which the electric vector is in the plane of 
propagation. , 

The elementary theory of wave-guides shows that waves can be 
propagated in a series of different modes. When both surfaces of the guide 
are perfect conductors, there is one mode in which the electric field is 
constant across the guide and perpendicular to its surfaces. This mode 
suffers no attenuation at any frequency. It will be called the mode of 
order zero. Higher-order modes are propagated without attenuation 
if the frequency is greater than the “ cut-off ” frequency, which is different 
for each mode, but for frequencies below the “ cut-off’ there is heavy 
attenuation. The “ cut-off’? frequency depends on the width of the 
guide, and the order of the mode. 

It is tempting to suggest that this “ cut-off’ accounts for the great 
attenuation observed on the lowest frequencies, but it can readily be seen 
that this simplified model is not in accord with observation. Thus, if 
we assume a distance of 50 km. between the two boundaries, the first- 
order mode would cut-off at a wavelength of 100 km., i.e. 3 ke./s., the 
second-order would cut-off at 6 ke./s., and so on. The emf. received 
at a great distance from a lightning flash, on a receiver tuned to 10 ke./s. 
would be comprised predominantly of four modes of orders 0, 1, 2, 3, 
while the e.mf. from the same atmospheric at 5 ke./s. would include 
two modes of orders 0 and 1. The resulting e.m-f. at either frequency 
would depend on the relative phases of the e.m.f.s in the various modes, 
but the ratio of the resultant amplitudes received at these two frequencies 
would clearly be of the order of 2: 1. The zero-order mode is propagated 
without attenuation, so that it is present at all frequencies, and this 
model therefore predicts that the signal should have appreciable amplitude, 
however low the frequency to which the receiver is tuned. But the 
observed amplitude at 10 ke./s. is very much greater than twice the 
amplitude at 5 ke./s., and at still lower frequencies the e.m.f.s from distant 
perceptible, showing that there can be no 
Moreover, on this model a curve showing 
hould have a series of 


atmospherics are almost im 
unattenuated zero-order mode. 
the variation of amplitude with frequency s 
approximately equal steps at, say, 3 ke./s., 6 ke./s., 9 ke./s., corresponding 
to the cut-off frequencies of the various modes. The curve plotted 
from the observations (Gardner 1959, fig. 4) shows a single steep increase 
at about 6 ke./s. There are no similar increases at 12 ke./s., 18 ke./s., ete. 

We can therefore reject the model using a wave-guide with perfectly 


reflecting surfaces for two reasons -—— 
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(a2) Observations show that at very low frequencies the signal 
amplitudes from distant atmospherics are much smaller than they 
would be if an unattenuated zero-order mode were present. 


(b) The observed variation of signal amplitude with frequency is of 
the wrong form. 


As the next approximation we therefore assume that the ionosphere 
is a homogeneous medium, whose plane lower boundary forms the upper 
surface of the wave-guide which has a plane, perfectly conducting earth 
as its lower surface. The effect of the earth’s magnetic field is neglected. 
The medium is assumed to contain N electrons per c.c., and these are 
assumed to make an average of v collisions per second. Let p be the 
angular frequency of a component of the signal. At the level where 
very long waves are reflected v is of the order 107 sec.-1 We shall, | 
therefore, assume that v>p. Then, if » is the refractive index of the 
medium, we have 


p?=1—)p,/p, 4 io 2 
where . p,=4nNe?/mv, 
e=charge of electron, 
m=mass of electron. 


We shall now investigate the modes which can be propagated in this 

guide. We consider only waves in which the magnetic vector is horizontal. 

Any one mode can be considered as composed of two component plane 

waves * equally inclined to the boundaries of the guide. Either of these 

may be considered as being continuously produced by the reflection of 
the other at a bounding surface. After a part of a wave-front has 

undergone two reflections, it is parallel to its original direction, and the 

phases and amplitudes of the original and of the twice reflected wave 

must agree. This means that only certain discrete values of the angle of 

inclination of the wave normal are possible. The excess path travelled 

by the wave-front of the twice reflected wave is 2h cost. The difference 

of phase of the original and of the twice reflected wave is therefore 
2(cos 7)h . p/e (this is, in general, complex), and the (complex) ratio of 
the amplitudes resulting from this excess path is exp {—j(2 cosi.h.p/c)}. 

There is also a change of phase and of amplitude due to the reflection 

of the component plane wave at the upper boundary, and if the (complex) 

reflection coefficient is R, then for a self-consistent mode, we shall have 


R . exp {—j(2 cos i. h pic)}=1=e™, 


where 7 is an integer and is the phase difference, measured in cycles, 
between the original and the twice reflected wave. 


* The two component plane waves which form a single mode are not, of course, 


to be confused with the direct and reflected rays of different order which are 
considered in the ray picture. 
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Now R is the Fresnel reflection coefficient of the wpper surface for angle 
of incidence i. It is given by 


R= {p? cos 1—4/(u?— sin? 1)}/{u? cos i++/(w2— sin?1)}, . . . (2) 
so that we obtain 
{u? cos t— 4/(u?— sin? 2)}/{u? cos t+-+/(u2— sin? z)} 
=exp {j(2 cost h p/lce+2mn)}, . . . (8) 
where h is the width of the wave-guide (height of the ionosphere), c is 
the velocity of electromagnetic waves in free space, and n is a positive 
or negative integer which is the “ order’ of the mode. 
This formula, in slightly different form, was given by Eckersley (19382, 
equation 3.21). From it the angle of incidence i can be found for each 


value of n, and in general 7 is complex. The attenuation coefficient of 
the mode is then easily found, and is equal to 


PSD ATR TC ieee, By see eames (2) 


(* # ”’ means “ imaginary part of ”’). 

It is still necessary to find the amplitude which is excited by a given 
source in each of the modes. For a vertical electric dipole, the amplitude 
radiated along a ray at an angle 7 to the vertical is proportional to sin 1. 
This means that in a mode in which one component wave has its wave 
normal at an angle i to the vertical, the amplitude excited by a vertical 
electric dipole is proportional to sin 2. 

The amplitude also depends upon the frequency spectrum of the 
source. For a source of the type described in fig. 1, the amplitude is 
proportional to the frequency. On general physical grounds, therefore, 
we should expect the amplitude excited by a vertical electric dipole in 
any one mode to be proportional to 


DLsi0 b, Geieee a eose ta oO) 


where 7 is given by equation (3). A more mathematical proof of this is 
given in Appendix A. 

If R=1 (perfect reflector), sin i has the same value for modes of order n 
and —n. The nth mode is then really a degenerate case, where the two 
modes of orders n and —” are superimposed. The degeneracy is removed 
when R41. The use of a negative order to describe a mode is merely 
a matter of convenience. All modes, whether of positive, negative, or 
‘zero order, have the same general form. It is, however, important to 
stress that, by inserting values of n of equal magnitude but opposite sign 
in equation (3), we derive the properties of different modes. 


§ 5. ATTENUATION CURVES. 
Our main reason for considering the finite conductivity of the upper 
boundary was to see whether it alters the attenuation of the zero-order 
mode. - The formule (3) and (4) have therefore been used to calculate 
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the attenuation coefficient as a function of frequency. Equation (3) does 
not give a convenient explicit expression for cos 7, and a method of 
successive approximation was therefore used. The results for the zero- 
order mode are shown in fig. 2 for h=50 km. The curves show the 
fractional reduction of signal amplitude for a path of 1000 km. in the 
wave-guide. (A similar curve has been given by Hales (1948), loc. cit., 
fig. 4, curve, 5, but it refers to a mode of non-zero order, and no curves. 
are given for the zero-order mode.) 

The parameter p, determines the nature of the upper boundary. When 
p,>, the upper boundary becomes a perfect conductor, and the zero- 
order mode should be unattenuated at all frequencies. It will be seen that, 
for large values of p,, there is very little attenuation at any frequency. 
For smaller values of p,, the medium is an imperfect conductor, and there- 
is then some attenuation which is most marked at frequencies between. 


Fig. 2. 


Amplitude reduction in 1000 km. path 


Frequency in ke/s. 


Fractional reduction in received e.m.f. for 1000 km. path in the zero-order mode 
. ae bl ts 
as a function of frequency. h=50 km. 


2 and 8 ke./s. The increase in the attenuation as the frequency decreases. 
from about 10 to about 5 ke./s. is very similar to that observed in 
the experiments mentioned in §1. The frequency near which this 
decrease of amplitude occurs is seen to be fairly independent of the 
characteristics of the upper surface. It does, however, depend markedly 
on the width of the guide (height of the ionosphere), as is illustrated. 
in fig. 3 in which the various curves all relate to the same value of p 
but different values of the height h. , 
When p, is zero, the upper medium is a “ dielectric ”’ with dielectric 
constant unity. The reflection coefficient is zero, so that all the energy 
incident on the upper boundary is propagated into the medium. If p, is 
very small, a large proportion of the energy is propagated into the medium 
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and all the modes are heavily attenuated, and are “leaky”. (For a full 
discussion of “leaky ’’ modes, see Booker and Walkinshaw (1946)). 
It would be necessary to add the effects of a very large number of modes to 
calculate the received signal in this case. However, the order of 
magnitude of N and v for the ionosphere are known, and these show 

that it is unlikely that we shall be concerned with very small values of 7,. 
' The form of the curves of figs. 2 and 3 depend only on the two. 
parameters p,andh. A small change of p, would not result in any marked 
change in the frequency where the attenuation curve has greatest slope. 
A change of 10 km. in f would have an appreciable effect on this frequency. 
It has been shown (Gardner 1950, Bowe 1951) that the curve relating: 
attenuation to frequency changes from night to day and during a Sudden 
Ionospheric Disturbance, in such a way that the cut-off frequency 


Fig. 3. 


Amplitude reduction in 1000 km. path 
Oo 
H 


re) 5 10 15 20 
Frequency in ke/s. 


Fractional reduction in received e.m.f. for 1000 km. os) in ieee zero-order 
mode, as a function of frequency. p,=6 x 10* sec. 


becomes greater. It is also known (Best, Ratcliffe and Wilkes 1936 ; 
Budden, Ratcliffe and Wilkes 1939 ; Bracewell and Straker 1949) that 
the height of reflection of very long waves falls in the ot ies 
night to day, and during a Sudden Ionospheric Disturbance. Itt se oe 
seems probable that the changes i Fs ee a curves are produce 
1 he changes in height of reflection. 

eee or ron Ride Oi Des itis necessary to find how many nee ie 
be used to calculate the received e.m.f. at a given disuates ay 
source. The results of calculations for the modes of so =e =. , 
0, +1, and for a distance of 1000 km., are shown in fig. 4 ne ; = 
fie —6x<104 sec.~! and h=—50 km. These modes oer 4 ots 
frequencies of 6 ke./s., 3 ke./s., 0 ke./s., 3 ke./s. respectively, 11 Dota gui 
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-surfaces were perfectly conducting. Fig. 4 shows that for frequencies 
up to about 16 ke./s. the signal amplitudes at 1000 km. from the source, 
in the modes of orders —2 and +1, are much smaller than in the modes 
of order 0 and —1. It is readily verified that the amplitudes in modes _ 
of still higher orders are much smaller still. Hence, for an explanation 
of the experimental results, it is sufficient, as a first approximation, to 
consider modes of orders 0 and —1 only. 


§ 6. THE WAVE-FORM OF AN ATMOSPHERIC. 


In the introduction we mentioned two methods of studying the 
-characteristics of the propagation path of an atmospheric. These 
were (a) a determination of the frequency response characteristic of the 
path (Gardner 1950, Bowe 1951); (b) the study of the wave-form of 


Amplitude reduction in 1000 km. path 


Frequency in ke/s. 


Fractional reduction in received e.m.f. for 1000 km. path in modes of various 
order (n). h=50km. p,=6 x10! sec.-} 


an atmospheric from a distant source. These two methods may at first 
appear to be unrelated, but in fact they depend on the same properties 
of the path. For, if the curve showing variation of amplitude with 
frequency is given for the path, the response of the system to an impulse 
can be calculated*. Hence, any model which is inadequate to explain 
the variation of amplitude with frequency is also inadequate to explain 
the oscillatory wave-form of an atmospheric from a distant source. 


* It would seem that we should also need to know how the phase varies with 
frequency, but in fact the amplitude and phase characteristics of any system 
are interrelated, and if the amplitude curve is given, it is possible, subject to 
certain restrictions which do not apply here, to find both the phase characteristic 
and the impulse response. See, for example, Murakami and Corrington (1948). 
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Hales (1948) has attempted to explain the oscillatory type of atmospheric 
wave-form by considering the distortion of an impulse as it travelled 
through a wave-guide with perfectly conducting walls, but he took 
no account of the zero-order mode. The zero-order mode must be present, 
and since it would be unattenuated at all frequencies, it would not produce 
any distortion, and so there should be a large undistorted impulse at the 
beginning of an atmospheric of the oscillatory type. This is not observed, 
and so we cannot accept Hales’ explanation without modification. We 
have now also shown that the observed variation of amplitude with 
frequency cannot be explained if the guide has perfectly conducting 
walls, since the zero-order mode is not attenuated. It is therefore 
necessary to assume that one surface of the wave-guide has finite 
conductivity. 

We shall now examine the wave-form which we should expect to 
receive after an atmospheric has travelled for a considerable distance 
(1000 km.) through such a guide. The calculation has been made for the 
modes of order 0 and —1, and for a source of the type described by 
fig. 1. For this purpose, we have to evaluate the integral 


I(t) = { p sin 4 . eHi,(kp sin 1) .j . dp 


(see Appendix B), where 7 is given by equation (3), p is the distance 
travelled, and k=p/c. It is sufficiently accurate to use the first term of 
the asymptotic expression for the Bessel function Hi, and with this 
approximation the integral reduces (omitting a constant factor) to 


| (p sin i)te77/4 exp {j(pt—kp sin t)} . dp. 


The integrals were evaluated numerically with the help of Beevers- 
Lipson strips (Beevers and Lipson 1936). The results for the modes of 
order 0 and —1 were found to be identical within about 10 per cent. 
This was not unexpected in view of the similarity of the two curves for 
these modes in fig. 4. Modes of higher orders have much smaller 
amplitudes (see fig. 4) and therefore are not included. Fig. 5 shows the 
result for either mode, and this is also the form of the e.m.f. which results 
when the two modes are added together. It is seen that the wave-form 
is oscillatory, and that the instantaneous frequency and amplitude 
both decrease with time. This behaviour is typical of many atmospheric 
wave-forms observed in practice. This calculation extends and completes 
that of Hales (1948), in which both walls of the wave-guide were assumed 
to be perfect conductors. aa 

Hales has suggested that the characteristics of the zero-order mode at 
very low frequencies in a wave-guide with imperfectly conducting elle 
may explain the “ slow tail’? which is often observed to follow the mam 
part of an atmospheric wave-form. He has shown that the observed 
delay between the beginning of an atmospheric and the peak of the 
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“slow tail’ is consistent with a reasonable value of the width of the 


guide. (Hales’ discussion of the main part of the atmospheric appears 


to refer only to modes of non-zero order, and his discussion of the “ slow 
tail” only to the zero-order mode.) 


400 —— 600 


Time in Microseconds ——> 


Received E.M.F. 


Received e.m.f. in mode of order 0 or —1 at 1000 km. from idealized lightning: 
flash. h=50km. p,=6X104sec.-1 


Fig. 6. 
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Received E.M.F. 


Received e.m.f. in plane wave resulting from the refection of impulsive plane- 
wave at the surface of an ionized medium. Angle of incidence 60°. 


§ 7. EQUIVALENCE OF “ Ray” PicTURE AND “ MopE” ProTuRE. 


The calculations of the received signal, treated either by the ray 
method or by the mode method, can both be made exact, and must 
therefore give the same result. It is, however, interesting to discuss the. 
interrelation of the two methods. It would seem at first sight that the 
ray method would result in the arrival at the receiver of a series of sharp. 
impulses resulting from varying numbers of reflections of the origina 
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impulse at the earth’s surface and at the ionosphere, and it is difficult 
to see how a series of sharp impulses can combine to give a smooth curve 
like that of fig. 5. It must be remembered, however, that the Fresnel 
reflection coefficient (equation 2) is a function of frequency, and a single 
reflection therefore results in distortion of the incident signal. Fig. 6 
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shows the form of this distortion for reflection at an angle of incidence 
of 60° from the plane surface of an ionized medium of the type already 
discussed, with p,—6 104 sec.-!_ Hence it is clear that the successive 
reflected waves are no longer sharp impulses, but have comparatively 
smooth wave-forms, which, when added together, would give a smooth 
curve, such as that of fig. 5. 
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If both surfaces of the wave-guide are perfect conductors, there is. 
no distortion of the signal at reflection, and the received e.m.f. must then 
be a succession of sharp impulses. But this succession of impulses must: 
also result from the addition of a large number of modes. The way in 
which this occurs is illustrated in fig. 7. Here the curves are plotted 
assuming that both guide surfaces are perfect conductors, that the source 
and receiver are both at the earth’s surface, and that h=50 km. The 
source is assumed to be an impulse of the type described by fig. 1. There 
is an analytic expression for the impulse response for one mode of order n 
for this case, at a distance from the source. If p>h and ct—p<p, this. 
can be shown (see Appendix B) to be approximately proportional to * 


{cos (7Xn/h)+(7Xn/h) sin (wXn/h)}/X8, . . . . (6) 


2 


where X?=c?t?—p 

In fig. 7 (a) the values of this function for the first three modes are plotted ; 
fig. 7(b) shows the result of adding together the first seven modes (n=0 
to 6 inclusive). Fig. 7(c) shows the succession of double impulses that. 
would be received on a ray picture. With only seven modes, the first. 
and second reflected impulses are already becoming accentuated, and it 
is clear that, by including a sufficient number of modes, the mode picture 
and the ray picture would give the same result. In this case the high- 
order modes are not attenuated, and therefore the mode method is not 
the best for calculating the received e.m.f. 

We have thus shown that, when the number of modes contributing to. 
the received e.m.f. is large, the resulting signal is a succession of impulses, 
so that the wave-forms are of the first type mentioned in the introduction. 
If the number of contributing modes is small, the received e.m-f. is. 
oscillatory and therefore of the second type. The transition from the 
first to the second type could occur as the number of contributing modes. 
decreased. Thus we might expect that atmospherics from nearby sources 
would be of the first type, while those from more distant sources would 
be of the second type. 


§ 8. PossiBLE EXTENSIONS OF THE THEORY. 


Since we are concerned with very long waves, it is probable that the 
ion density and collision frequency in the ionosphere vary greatly in 
a distance small compared with one wavelength, and therefore, in the 
theory of propagation of very long waves, many authors have treated 
the ionosphere.as a sharply bounded, homogeneous medium, as a first 
approximation. In discussing the propagation of very long radio waves 
to great distances, Rydbeck (1944) has used the method of Watson 
(1919), which is essentially a treatment in terms of wave-guide modes. 
Both these authors allow for the curvature of the earth’s surface. Rydbeck 


* A series of terms like (6) is not uniformly convergent but can be shown to- 
be “summable (C2) ” (see Whittaker and Watson 1935, Modern Analysis, 
para. 8.43). If the electric dipole at the source were destroyed in a small but: 
finite time, the resulting series would be uniformly convergent. 
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has discussed the cases where the ionosphere is assumed to be (a) a 
sharply bounded, homogeneous medium, (b) a medium in which the ion 
density has a maximum and varies parabolically with height, (c) a medium: 
in which the ion density varies as the square of the height Boots its lower 
edge. He concludes that, for the modes of lowest order, there is not much 
difference in the behaviour for the three cases. Rydbeck gives numerical 
results for only a few single frequencies. His formule could be used 
for calculating curves showing attenuation as a function of frequency, 
but the calculation would be laborious because of the complexity of the. 
formule. G 

No allowance has been made for the effect of the earth’s magnetic 
field. If this effect is included, the ionized medium becomes doubly 
refracting. It is still possible to consider propagation in terms of wave- 
guide modes. When a plane wave is incident from below on the boundary’ 
of the medium, the state of polarization of the reflected wave is in general 
different from that of the incident wave. For a given mode, it is necessary 
to ensure that a component wave, after two reflections at the wave-guide. 
surfaces, has not only the same amplitude and phase as the original 
component wave, but also the same polarization. The result is that the 
number of possible modes of a given order is doubled, and the electric 
field of the wave is no longer wholly in the plane of propagation. 
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APPENDIX A. 


Wave-GuiIpE MopES EXCITED BY A VERTICAL DIPOLE. 


In this appendix, we deduce an expression for the horizontal magnetic 
field produced in a wave-guide of the type discussed in the paper, by a 
vertical electric dipole of constant angular frequency p at the surface of 
the earth. In Appendix B, we shall use the result to deduce the field 
when the source is a lightning flash of the type shown in fig. 1, and when 
both surfaces of the wave-guide are perfect conductors. 

Suppose the source is a vertical electric dipole of moment P=exp (jpt). 
‘Then we shall show that the horizontal magnetic field %, at a horizontal 
-distance p from the source and at height z, is proportional to 


+ 0 
x p®?.oc,Hi,(k . op) . cos {(j¢,—27n)z/2h} exp (jpt). . . (7) 

Each term of this series is the field in one mode. o, is the value of 
:sin?, for the nth mode when n is zero or positive, and the complex 
conjugate of sini when 7 is negative. We have written R=e~¢ and q, is 
the value of g when i has the value appropriate to the mth mode. The 
angle i is given by equation (3). Hi is Jeffreys’ notation for a Bessel 
function of the third kind (Jeffreys and Jeffreys 1946, Mathematical 
Physics (Cambridge), p. 544). 

The expression (7) can be predicted on general physical grounds. 
‘Thus the factor p?c, when multiplied by 1/p (which is proportional to 
the amplitude of the frequency p in the dipole moment of the ideal lightning 
flash) gives psini which is simply the expression (5) of §4. The 
cosine term is a “standing wave” term which arises through the 
combination of the upgoing and downgoing component waves of each 
‘mode. The Bessel function represents a wave propagated outwards 
from the source with cylindrical symmetry. For large values of p, it is 
proportional to {exp (—jko,,p)}/1/(ko,,p). 

The mathematical derivation of (7) is given below in outline. The 
method used is similar to that developed by Watson (1919). 

We use spherical polar coordinates r, 0, 4, with the dipole source at 
the origin. We also use z for distance measured vertically upwards, 
and p=1/(x*+y") for the horizontal distance measured radially from the 
origin. Then the field of the dipole in free space is given by the Hertz 
vector [T)=exp (jpt) . {exp(—jkr)}/r. This may be written 


; /2-+4 0 a 
IT)=(jk/27) . exp (jpt) 5 ; dd ol (S/C) dS . exp {—jk(a cos ¢ sin 6 
2—joo —-o 


+y sin ¢ sin 6+-2z cos 6)}, 
-where S=sin 6, C=cos 0. 

The integrand represents a plane wave whose wave normal has the 
direction (0, ¢). The radiated field is thus expressed as the sum of an 
infinite number of such waves. Note that it is necessary to include 
complex values of ¢ and @. 
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Each plane wave undergoes successive reflections at the earth and the 
ionosphere. To get the complete resultant Hertz vector, IJ, we must add 
to the original Hertz vector the expressions for the various reflected 
waves. It is permissible first to perform the integration with respect 


to ¢, since the reflection coefficient, R, of the ionosphere is independent: 
of 6. Then we have 


IIp=}jk exp (jpt) . | Hi,(kSp) . exp (—jkOz) . (S/O) . dS. 


The phases of the successive reflected waves are easily seen to contain’ 
factors exp (+j .2hkC), exp (+7 .4hkC), etc., since they originate at: 
image sources at heights +2h, +4h, etc. Hence the complete field is 
' given by 


SIC. dS Hiy(kSp) 


co 
=O: 


IT=jk exp (Jp!) | 


x[ y R’ . exp (—j . 2vhkC)-+ 2 R” exp (+p) | . 
vy=0 y=1 

' 
By suitable manipulation of the contours, this series may be transformed. 
into another series, corresponding to the “residue” series of Watson. 


The result is 


= (jmnj2h) exp (jpt)| | E _Hig( 8 qp) 005 ((jtq—2rm 2h} 


4S Hi,(kS,,p) - cos {ign 202/283 |. Se cee) 
n=0 A“ 


Here we have written R=e-¢. An asterisk * denotes a complex conjugate. 
The terms of the last series are the various wave-guide modes. The 
values S,, are the poles of the integrand, and are given by 


Dk, =n 27M, - 2 2 ee et (9) 


where C2=1—S?. ; wile . 
Equation (8) is derived by an application of the residue theorem 
of complex variable theory. We note at this stage how the full 


mathematical treatment gives two results :— 
(a) the characteristics of each mode; these arise from equation (9), 


which leads directly to equation (3) and gives the positions of 
the poles of the integrand in the complex S-plane ; 


(b) the amplitudes excited in each mode by the source ; these are the 
terms of equation (8) and are the residues of the integrand at 
the poles. pe’, 


c 
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If the ionosphere is a perfect reflector, g=0, and the §,’s are all real 
or purely imaginary. Then (8) becomes 


IT=(jn/2h) . exp (jpt) | His(kp)-+ 2 2 cos (rnz2/h) . Hig(8,0) | .. (10) 


n=1,2.. 


The summation includes all values of x for which §,, is real. Terms for 
which S,, is imaginary refer to heavily attenuated (evanescent) modes, and 
may be neglected when p is large. 

The horizontal magnetic field W of the signal is derived from the Hertz 
vector thus 


H =G7II/dt . Op. 


Then for the mode of order n, we have, from (8) (imperfectly reflecting 
ionosphere), 


8, : Ss ; 4 n>0 
Ea — ae (wp?/2hc) Hi, {kp ‘oe cos {iq¢,—27n)z/2h} exp (jpt) .... a . 
‘This leads to equation (7). 
For a perfectly reflecting ionosphere, we have, from (10), 


KH ,=(7pS,,/he) . Hi,(kS,p) . cos (7nz/h) exp (ppt). 


APPEND TX —B: 


_ IwpuLse RESPONSE OF ONE MODE. 


Equations (8) and (10) give the Hertz vector of the signal when the 
source is a vertical electric dipole at the earth’s surface, oscillating with 
constant angular frequency, p, and constant amplitude. 

If the source is a lightning flash as illustrated in fig. 1, the dipole 
moment of the source (fig. 1(a)) can be expressed as the Fourier integral 


PQ)=|-_ {exp (jpi)}ip.dp. 


Each component frequency in the Fourier spectrum initiates a series of 
modes such as that of equation (8) or (10). 

We consider only one mode, the nth, for the case where both earth 
and ionosphere are perfect conductors (equation 10). We assume that 
the receiver is at the earth’s surface, so that z=0. We use the value 
of 8, given by (9) when q=0. The Fourier component of angular 
frequency p then gives, for the nth mode only, 


I1=(jn]2h) exp (jpt) . 2Hig{kp-/(1—n?n?/B12)}. 
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‘The complete impulse response for the nth mode is therefore 


IT ,=(x/h) jee {exp (jpt) }/jp «dp . Hig {kp./(1—1?n?/h?k?)}, 


-Hence 
II, /Ot—=(n/h) | : exp (jpt) . jdp . Hig {s/(p?—n2n2e?/h2) . ple}. 


‘This integral can be reduced to a standard form (see, for example, 
McLachlan and Humbert 1941, Mémorial des Sciences Mathématiques, 
Fascicule 100, p. 34). Its value is 


OlT,,,/Ot=(ae/hX)[ cos (7nX/h)— (2/77) sin (wn X/h) log, {(ct-+X)/p}], 
-where X?=c?t?— p?. 

‘To find the horizontal magnetic field which the nth mode contributes 
to the received signal, this expression must be differentiated again with 
respect to p. We are interested in values of X which are small compared 
with ct and p. We also assume p>h. Then it is readily shown that the 
log. term is negligible, and we have 


KH =P, [dtd p= (mep/hX*) {cos (mmX/h)-+(mnX/h) sin (nn X/h)}. 
‘This leads to equation (6). 
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ABSTRACT, 


Measurements of the attenuation of a beam of high energy neutrons 
in carbon and polythene have been made in good geometry. The total 
cross-sections of hydrogen and carbon have been found to be 
0:0464-+0:0012 x 10-24 em.? and 0-330-+.0-003 x 10-*4 em.” respectively, 
at an effective neutron energy of 156+3 MeV. 


INTRODUCTION. 


THE present work has been the first of a series of experiments on high 
energy neutrons produced by the 110-inch frequency-modulated cyclotron 
(Pickavance, Adams and Snowden (1950)) at the Atomic Energy Research 
Establishment, Harwell. In the course of the measurement of the neutron— 
proton cross-section, the total cross-section of carbon was also measured 
and this result is included. The experiments have already been reported 
briefly (Taylor, Pickavance, Cassels and Randle (1950)), and the purpose: 
of the present paper is to describe the work in more detail. 

The total neutron cross-section of hydrogen has been measured at. 
40 MeV. and 90 MeV. by Hadley, Kelly, Leith, Segré, Wiegand and 
York (1949), at 95 MeV. by De Juren and Knable (1950), and at 260 MeV. 
by Kelly, Leith, Segré and Wiegand (1950); recent measurements by 
Moyer and de Juren (1950) at energies between 90 MeV. and 270 MeV. 
were published almost simultaneously with our preliminary report. The. 
experimental results so far obtained are in mutual agreement within 
experimental error. 

The method of observation in the work described here was similar to. 
that used at Berkeley by Hadley et al. (1949) and by Kelly et al. (1950). 
Measurements of the attenuation of a neutron beam in carbon and poly- 
thene were made in good geometry. The neutron detector was a triple- 
coincidence counter telescope, into which protons were scattered from a 
polythene disc placed in the neutron beam. 


*Communicated by the Authors. 
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APPARATUS, 


Collamation of the Neutron Beam. 
A beryllium target, 8-9 mm. thick, was mounted in the cyclotron vacuum 
chamber near the edge of the dee at a radius of 484 inches, where it 
intercepted a beam of about 1 microampere mean current of 171 MeV. 
protons. Neutrons emerging near the forward direction (fig. 1) passed 
through a thin aluminium window in the vacuum chamber wall, and were 
then collimated by a series of 5 cm. diameter holes in concrete blocks and 
_ shielding walls. The shielding, which was divided into three parts, had 
a total thickness of 4:3 metres. The neutron attenuators were placed 
between the second and third walls and the neutron detector was in a long 
tunnel beyond the third wall. The collimating holes, the polythene 
radiator and the neutron attenuator holder were aligned with the target 
with the help of a theodolite. 


Fig. 1. 
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The proton radiator (R, fig. 1) was a polythene disk 19 mm. bea ek 
6-3 cm. in diameter, placed at a distance of ut metres from ai are : 
Protons ejected from the dise at an angle of 10° to the jopeen irec ae 
were detected by a telescope consisting of three proportiona clea sit 
€,, Cz, Cz; separated, by carbon absorbers A, Ay. All ei ped sea : 
clear of the main neutron beam. Each counter was aye rica an i p : 
5 cm. in diameter and 5 cm. long, with a 46 mg./cm.” copper ie “ a) 
each end and a 0-02 cm. tungsten wire stretched across the diameter. 
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The gas filling was a mixture of argon (96 per cent) and carbon dioxide 
(4 per cent) at a total pressure of 52 cm. of Hg and the operating voltage 
was about 2 kilovolts. 

Each of the counters was connected to a linear amplifier channel (fig. 2) 
with a high frequency response up to 1-5 Mc./s., consisting of a head 
amplifier feeding through about 50 yards of correctly terminated coaxial 
cable to a main amplifier in a remote counting room. A differentiation 
time constant of 0-6 microseconds was used and all coupling time constants. 
in the later stages of the amplifier were kept short to reduce the dead 
time following a large background pulse. With the actual neutron 
background encountered there was no difficulty from dead time effects. 

The output pulses of the main amplifiers were mixed in a standard 
triple coincidence unit Type 1036A. The “ tuning ” of this coincidence 
unit was carried out in the usual way, by adjusting gate-open times and 
delay times in the individual channels until all real coincidences were 


Fig. 2. 
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counted ; the gate-open times finally used were 0-5 microseconds. The: 
amplifier and discriminator were adjusted to register a count on the arrival 
of 10,000 electron charges in a pulse. Standard scaling units type 1009A 
were used to register the mixed pulses and the counts in each channel. 

A second triple coincidence unit enabled the accidental coincidences. 
to be measured during the normal counting runs. Double coincidences 
between channels 1 and 2 were fed into the final mixer stage of the second 
coincidence unit, and mixed with the single pulses from channel 3. The 
pulses from channel 3 were delayed by 24 microseconds in a coiled inner 
coaxial line, so that no genuine coincidences could be recorded. There 
was a considerable counting rate of genuine double coincidences caused 
by recoil protons which had sufficient energy to pass through the first 
two counters, but were brought to rest in the absorber A, between the 
second and third counters. Thus the background consisted of accidental 
coincidences between the mixed channel 1 and 2 pulses and the random 
counts of channel 3 caused by relatively slow neutrons emerging from 
the shielding walls. ; ; 


——e 
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All three counter voltages were adjusted to be on a plateau of triple 
coincidence counting rate against voltage in any arrangement used, thus 
ensuring that the fastest protons were being counted and that any aight 
drifts of voltage had a negligible effect on the counting rate. 


The Neutron Attenuators. 


A revolving drum (D, fig. 1), carrying four tubes 1-6 metres long and 
6-3 cm. inside diameter, was slung between the second and third shielding 
walls, about 8 metres in front of the neutron detector. A carbon 
attenuator, 120-23 gm./cm.?, was mounted in one tube and a polythene 
((CH,),,) attenuator containing 136-93 gm./em.? of polythene in another. 
One tube was left empty, and a copper bar, sufficiently long to remove 
virtually all the neutrons from the beam, was fixed in a fourth position. 
The drum could be driven by remote control from the counting room, 
so that measurements in the various positions could be made in rapid 
succession. Impurities in the carbon and polythene were such as to 
have a negligible effect on the attenuation. A second pair of attenuators, 
containing 70-76 gm./cm.? of carbon and 82-41 gm./cm.’? of polythene 
respectively, was also used. : 


THe ErrectivE NEUTRON HNERGY. 


The triple coincidence telescope detected neutrons in two distinct ways, 
since protons were ejected from both the hydrogen and the carbon content 
of the polythene radiator. A test with a graphite radiator showed that 
during the total cross-section measurements 93 per cent of the counts arose 
from hydrogen recoils and 7 per cent from nuclear reactions in carbon. 
The hydrogen recoils were elastic so that neutron energies could be 
inferred directly from the energies of the protons actually detected. On 
the other hand, only a lower limit could be placed on the energy of the 
neutrons detected through nuclear reactions in the carbon. This lower . 
limit was fortunately rather high, because of the ?C(, p)B threshold of 
13-2 MeV. 

- Throughout the total cross-section measurements the absorbers A, 
and A, were kept equal to 11:26 and 3-35 gm./cm.? respectively. The 
carbon equivalent of all the copper windows of the counters, together 
with the air separating them, was 0-21 gm. jom.2. A further 1-00 gm./cm.? 
of carbon equivalent was added because the protons, on the average, 
traversed half the thickness of the polythene radiator. With the help 
of the range-energy relation in carbon (Taylor 1950) the minimum mean 
proton erergy detected was found to be 145 MeV. Since the telescope 
was set an at angle of 10° to the neutron beam the minimum neutron 
energy detected was, on the average, 149 MeV. Now, the energy of the 
primary protons in the cyclotron was 171 MeV., so that there were no 
neutrons of energy greater than 169 MeV. The effective energy of th 
experiment was thus quite closely defined from the outset. 
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If the total cross-sections measured are assumed not to vary rapidly 
with energy, then the effective energy can be fixed still more accurately. 
First, the effective energy EE, of the neutrons detected through hydrogen 
recoils was found by means of the following subsidiary experiment. The 
absorber A, in the telescope was increased in steps to a maximum of 
9-15 gm./em.? and the triple coincidence rate determined at each setting, 


Fig. 3. 
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first with a polythene radiator and then with a carbon radiator. The 
integral spectrum of the hydrogen recoils (fig. 3) was then found from the 
difference between these two sets of results. Small corrections of the 
order of 2 per cent had to be applied to allow for the increased nuclear 
absorption of the protons in the extra absorber. After translating 


hydrogen recoil energies into neutron energies, EH, was found to be 
155+3 MeV. 
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The effective energy E%, of the neutrons detected through nuclear 
reactions in carbon could not be found experimentally. This did not 
matter greatly because only 7 per cent of the protons came from this 
source. = 
The high n,; p threshold in carbon restricted the corresponding neutron 
energies to the comparatively narrow range from 162 to 169 MeV. The 
value of ES, was guessed to be 165+2 MeV. 

Giving E#, 93/7 times the weight of E%,, the overall effective energy 
of the experiment was found to be 156+3 MeV. This value has been 
slightly revised from that given in our earlier report. 

The estimated error of the effective energy includes some allowance 
for possible inaccuracies in the range-energy relation. The results 
displayed in fig. 3 show clearly that this source of error was small. The 
maximum neutron energy was 169 MeV., so that the maximum possible 
- hydrogen recoil energy was 164 MeV. The curve of fig. 3 should therefore 
reach the energy axis at 169 MeV., remembering that haif the radiator 
thickness has been added in arriving at the energy scale. In fact, the 
experimental intersection is at 169-2 MeV. 


THE Bram MONITOR, 

The neutron output of the cyclotron was monitored by counting the 
pulses in an ionization chamber placed in a standard position behind the 
second shielding wall (fig. 1), where it was not affected by changing the 
neutron attenuator. Tests showed that the monitor counting rate was 
proportional to the cyclotron beam current. 


EXPERIMENTAL PROCEDURE. 


A complete cycle consisted of three successive measurements of the 
ratio of triple coincidence counts to the monitor counts with a polythene 
attenuator, a carbon attenuator and finally with no attenuator. The 
duration of counting was so arranged that about 300 triple coincidences 
were recorded in each run, and with a proton beam of 1 microampere 


mean current a complete cycle could be completed in about half an hour. 


Occasionally, between cycles, a blank run was performed with the long 
copper attenuator. In all, 19 cycles of readings were taken with the 
longer carbon and. polythene attenuators, and 6 cycles with the shorter 
ones as a test for ‘“ hardening ” and other effects which might have 
introduced errors. 
Accidental triple coincidences represented a background which had to 
be subtracted from the counts observed during the experiment ; they were 
measured by the method described above. Although the single channel 
counting rates were all far greater than the real triple coincidence counting 
rate, the background was only about 1 per cent with an. attenuator in 
position. The background represented by the counting rate bee the 
long copper attenuator was negligibly small, showing that there were no 


serious leaks in the shielding. 
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RESULTS. 


The measure of the incident neutron beam upon the counter telescope 
has been taken as the number of real triple coincidences divided by the 
monitor counts. “ 

Calling the ratio Ty, T,, T, where the subscripts refer to the beam 
unattenuated, and attenuated by carbon and polythene respectively, 
then : 


In-T,/T.=(No); 
and In T,/T,=(Ne)p, 
where N is the number of atoms per cm.? of attenuator and o their cross- 
section. Hence 
In T,/T,=(No),,—(No),. 
If (4No), represents the excess of carbon in the polythene over that. 
in the carbon attenuator, then : 


(No)y=In T,/T,,—(ANo),, 


where the subscript H refers to the hydrogen content of the polythene. — 
The results are presented as the ratios of T,/T, and T,/T,, for each cycle 

of the experiment. The mean value of these ratios then is used and the 

standard deviation quoted is that determined from the deviations from 

this mean. The values of the standard deviations so found differ little 

from the values expected from the total number of counts in the whole 

experiment. 


TABLE I. 


Short Attenuator. 


Cycle L./T, Ty/T. 
A 1-360 3-780 
B 1-436 2-981 
0 1-320 2-919 
D 1-413 3-598 
E 1-312 3-373 
F 1-412 2-890 


From the results in Table I, the values of T, Lee atid k /T,, are 
1:376+0-025 and 3-257-+-0-182 respectively ; the Baa Pile re 
expected from the total number of counts are 0-044 and 0-104 respectively. 


With these values the total cross-sections o,, oy for carbon and hydrogen. 
are : 


o,=0-333+0-016 x 10-74 em.? 
o=0-0454-+0-0025 x 10-4 em.? 


of Hydrogen and Carbon for H igh Energy Neutrons 27 


From the results in Table II. the m 
: ean values of T,/T, and T,/T, 
are 1-649+-0-0242 and 7:310-+-0-127 respectively, the peer Hatiatious 
expected from the total number of counts being 0-0271 and 0-098. 
respectively. The total cross-sections with these values are : 


6:=0°330-40-003 x 10-24 em.2 
oy—=0-0467+0-0014 x 10-*4 em.? 

The weighted means of these results are : 
o,=90:330+0-003 x 10-24 cm.? 


oy—0:0464+.0-0012 x 10-24 cm.? 


TABLE II. 


Long Attenuator. 


Cycle tyl p T/T. 
A 1-714 7-810 
B 1-596 7:364. 
C 1-913 7-139 
D 1-779 7-682 
K 1-402 8-406 
F 1-845 6:°747 
G 1-547 6-915 
H 1-625 7-674 
I 1-750 6-765 
J 1-640 6-807 
Ik 1-622 6-724 
L 1-672 7-235 
M 1-675 6-808 
N 1-600 7-431 
O 1-629 6-866 
P 1-558 7-812 
R 1-650 7:955 
g 1-493 7-004 
ap 1-683 
U iheralss 

CONCLUSION. 


This value of the hydrogen cross-section is in line with the results of 
experiments at Berkeley (Hadley et al. 1949, Kelly et al. 1949, De Juren 
and Moyer 1950, De Juren and Knable 1950). However, it is too low to 
be compatible with the theoretical total cross-sections calculated with. 
the interactions suggested by present meson theories (Burhop and Yadav 


1949, Christian and Hart 1950). 
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SUMMARY. 


Schlieren photographs have been taken to show the reflection of oblique: 
shock-waves of different strengths at the turbulent boundary-layer 
on the wall of a supersonic wind-tunnel. A weak shock-wave is reflected 
as a compression followed by an expansion; this is in agreement with 
Lighthill’s theory in which the boundary-layer is replaced by a non-turbu- 
lent inviscid layer with the same velocity distribution as a turbulent 
boundary-layer. For stronger incident shock-waves the reflected wave: 
moves upstream of the incident one to a point where the boundary-layer 
starts to thicken. Thus the incident and reflected shock-waves cross. 
outside the boundary-layer. Near the point where the incident shock- 
wave meets the boundary-layer, the rate of growth of the boundary-layer 
suddenly decreases and an expansion occurs. For very strong shock-waves 
a Mach type of interaction occurs, with a short normal shock-wave 
between two intersections. 

The supersonic wind-tunnel and optical system used for these 
experiments are briefly described. 


§1. INTRODUCTION. 


THE experiments were made in a small intermittent supersonic wind- 
tunnel at a Mach number of 1-965, using a wedge to produce an oblique: 
shock-wave. The strength of the shock-wave was varied by changing 
the angle of incidence of the wedge. The interaction between the: 
shock-wave and the turbulent boundary-layer on the tunnel wall was. 
investigated by schlieren photography. 


§ 2. Novation. 


p,=Sstatic pressure before incident shock-wave. 
p,=static pressure after incident shock-wave. 


Y= PoP =strength of incident shock-wave. 


2 . ° 
: $=stream deflection produced by incident shock-wave. 


* Communicated by the Authors. — 
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%—theoretical angle between undisturbed stream and incident 
shock-wave, calculated from a given value of 6. 

x,—=angle between undisturbed stream and incident shock- 
wave, measured at upstream edge of “fan” in photograph. 

x,—angle between undisturbed stream and incident shock-wave, 
measured at downstream edge of “ fan ” in photograph. 


§ 3. DescRIpTION OF WIND-TUNNEL. 


The arrangement of the wind-tunnel is shown in fig. 1. The large 
-reservoir M consists of three modified Lancashire boilers having a total 
capacity of about 130 cum. This reservoir is connected to the other 
parts of the wind-tunnel by the quick-operating valve L. Before starting 
the wind-tunnel the reservoir is evacuated to an absolute pressure of 
about 0-1 atmosphere. The valve L is then opened and air enters the 
intake A from the atmosphere, is dried by the silica gel C and passes 


Intermittent supersonic wind tunnel. 


through the contraction E to the working section F. The working 
section is rectangular, with a height of 12-7 cm. and a width of 10-2 em., 
the shaped nozzle blocks being at the top and bottom, with parallel 
glass plates at the sides. Downstream of the working section there is 
a gradually expanding passage G, in which the supports for the models 
are fixed. This is followed by a transition to a circular section H and a 
subsonic diffuser K leading to the valve L. 

In the drier the silica gel bed C is 19 em. thick and has a frontal area 
of 2:65 sq. m. With a Mach number of 1-965 in the working section 
(area 129 sq. cm.) the mean velocity of the air passing through the silica 
gel is only 0-58 m. per sec. This small velocity makes it possible to use a 
cotton-wool filter B to clean the incoming air and so prevent contamination 
of the silica gel by dirt from the atmosphere. A wire gauze D is placed 
after the silica gel to reduce the turbulence caused by the structure 

.supporting the silica gel compartment... The total pressure drop. caused 
by the cotton-wool, silica gel, and wire gauze is about 15 em. of water. 
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A small centrifugal fan and a 12 Kw. electric heater are provided for 
regenerating the silica gel as required by blowing hot air through it. 
The maximum duration of each run. of the wind-tunnel is about 
45 seconds. For many purposes a much shorter time is sufficient, and 
several runs can then be made with one evacuation of the reservoir. 
At a Mach number of 1-965 the pressure in the reservoir when the normal 
shock-wave begins to move upstream into the working section is 
0-665 times the stagnation pressure upstream of the throat. The diffuser 
efficiency may be defined as the work done in isentropic compression, 
between the condition before and the pressure after the diffuser, divided 
by the decrease of kinetic energy in the diffuser. Using this definition 
the efficiency is 0-75, a value consistent with the results obtained by 
Neumann and Lustwerk (1949) for a similar type of diffuser. | 

The distribution of static pressure along the centre line of one side of 
the working section was measured, using a brass side-plate fitted with 
pressure holes instead of the usual glass window. From these static 
pressures the corresponding Mach numbers were calculated, assuming 
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Distribution of Mach No. along centre line. 


—— 


that the stagnation pressure was equal to the static pressure in the drier 
immediately after the gauze D. (Measurements with a pitot tube at 
high subsonic speeds have shown that this assumption 1s correct, outside 
the boundary-layer, within 0-2 per cent of the dynamic pressure. ) The 
distribution of Mach number obtained in this way is shown in fig. 2; 
for the purpose of the present experiments the mean Mach number is 
assumed to be 1-965. 


§ 4. Humipirty MEASUREMENTS. 


The humidity of the air passing through the wind-tunnel cannot be 
‘measured in the limited time available while the tunnel is running, and 
an indirect method is therefore necessary. .'The tunnel is run about three 
times, starting with an absolute pressure of about: 0-1 atmosphere in, the 
reservoir and continuing each run until the pressure im the reservoir is 
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atmospheric. When this has been done nearly all the air in the reservoir 
has passed through the drier, and a sample of air from the reservoir is 
therefore representative of the air passing through the wind-tunnel. 
A compressor is then used to pump air from the atmosphere through the 
silica gel drier into the reservoir, so that the pressure is slightly greater 
than atmospheric. A steady stream of air from the reservoir is then 
passed through a frost-point hygrometer of the type described by Brewer, 
Cwilong and Dobson (1948). 

When the silica gel has just been regenerated the water vapour content 
of the dried air passing through the wind-tunnel is always less than 
1 part in 5000 by weight, and is sometimes very much lower. Normally, 
the silica gel is regenerated when the water content of the dried air has 
risen to about 1 part in 2500 by weight. ' 


§5. OPTICAL SYSTEM. 


The optical apparatus used for the schlieren photographs was a 
conventional two-mirror system as shown in fig. 3. The two spherical 


Schlieren, optical system. 


concave mirrors M, and M, have a focal length of 152 cm. and a diameter 
of 12-7 em. lL, is a condenser lens and L, is a long-focus achromatic 
object lens. For setting up and adjusting the apparatus, and for visual 
observation of the schlieren image on a screen, a “‘ compact source ”’ 
mercury vapour lamp is used. The very high brightness of this lamp 
gives good illumination on the screen, even when the second knife-edge 
is very close to the position giving complete cut-off. This facilitates 
exact adjustment of the knife-edges and mirrors. 


§ 6. PHOTOGRAPHS OF OBLIQUE SHOCK-WAVES. 

The wedge used to produce the shock-waves was 3-8 cm. long in the 
direction of the stream, had an apex angle of 10-1°, and was made of 
hardened high-carbon steel. It was ground and lapped to give a constant 
angle with as sharp an edge as possible ; the thickness of the leading edge 
was measured and found to be between 2:5 and 12u. The wedge extended 
across the whole width of the working section and was supported by two 
arms on the downstream side. 


BARDSLEY & MAIR Phil. Mag. Ser. 7. Vol. 42. Pl. I. 
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The schlieren system was arranged with horizontal knife-edges so that 
increasing air density in an upward direction corresponded to a light 
region in the photograph. 

Preliminary visual observations showed that the flow was steady, so 
that a relatively long exposure (=th second) could be used for the 
photographs. The mercury vapour lamp gave considerably more light 
than was necessary for photography, and a slow “ process” film was 
used, with a yellow filter to reduce its effective speed. This arrangement 
was convenient, as it enabled a single light source to be used for visual 
observation and photography. Figs. 4 to 9 (Pls. I. to IV.) are 
photographs of the flow for positions of the wedge giving deflection 
angles (8) of 2-9°, 5-0°, 7-3°, 9-8°, 12-1° and 15-2°. In each case the flow 
is from left to right. Table I. gives the strengths (y) of the incident 
shock-waves, calculated from the values of 5. 


§7. Discussion oF PHOTOGRAPHS. 


The arms suporting the wedge continue downstream beyond the edge 
of the field of view. The position of the back of the wedge is revealed 
by the two large expansion regions. The joints in the supporting arms 
produce two shock-waves at the right-hand side of each photograph. 
In addition, several wavelets from the liners are visible and a scratch on 
one of the glass walls may be seen in the bottom left-hand corner of each 
photograph. In fig. 7 the vertical black line across the middle of the 
photograph is a reference mark to show the position of the wedge in the 
tunnel. 

The thickness of the light region representing the boundary-layer on 
the tunnel wall is about 4 mm. at a point 24 cm. downstream of the 
throat. This indicates that the boundary-layer was turbulent, as would 
be expected at such a high Reynolds number (1-34 x 10° per cm. in the 
working section). In all cases the tunnel wall was flat at the point of 
incidence of the shock-wave. mM 

In each photograph the thickness of the line representing the incident 
shock-wave increases with distance from the leading edge of the wedge. 
This is probably because the flow is not two-dimensional, so that the 
- shock-wave is not plane over the whole section of the tunnel. The 
reduced velocity in the outer part of the turbulent boundary-layer on 
each glass wall increases the shock-wave angle ii that region and hence 
explains this fan-shaped appearance. Thus the downstream side of the 
“fan” represents the true position of the shock-wave outside the 
boundary-layer. This explanation is supported by fig. 4, showing the 
weakest incident shock-wave, in which the fan-shaped region is very 
faint except at the downstream edge. Similarly, in figs. 4 to 7, the pamie 
effect can be seen in the fan-shaped region representing the shock-wave 

ace. 
bernie ae deflections and shock-wave angles are considered. 
For a given stream deflection angle (5) the theoretical value of the 
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shock-wave angle (%)) has been calculated and compared with the 
measured angles at the upstream and downstream edges of the fan 
(x, and a). It can be seen that % agrees much better with «, than 
with «,, again showing that the downstream edge of the fan corresponds 
to the position of the shock-wave outside the boundary-layer. Part of 
the difference between the values of %) and «, may be explained by the 
growth of the boundary-layer on the wedge, causing an increase in the 
effective value of 6. 

In the light fan-shaped region representing the incident shock-wave in 
figs. 4 to 7, there is a dark line near the downstream edge. The reason 
for this is not understood but further work is being done to investigate 
the flow at the sharp leading edge of a wedge. 


TABLE I, 
Fig. no. 8 x Xo ay Xo 
4 2-9° 0-17 33-0° — 33-6° 
5 5-0° 0-31 34-9 37-4 35°6~ 
6 ives = 0-48 cba 39-5" 37-8° 
7 9-8” 0-68 39-8 42-1 40-1° 
8 PZ}? 0-90 42-6° 45-9° 43-4° 
9 15-2° 1-20 46-6° 49-4° 47-12 


Figs. 4 and 5 show that a weak shock-wave is reflected as a compression 
followed by an expansion, 7. e. as a pressure ridge. The dark region after 
the expansion indicates that there is a further compression there, with 
a smaller density gradient than in the first compression. The expansion 
appears to get weaker with increasing distance from the boundary ; this 
would be expected as a result of interaction with the compression waves. 
At a sufficient distance from the boundary the expansion probably 
disappears completely, and the reflection of a weak shock-wave is then 
simply another shock-wave. These conclusions are in agreement with 
Lighthill’s theory (1950) in which the boundary-layer is replaced by a non- 
turbulent inviscid layer with the same velocity distribution as a turbulent 
boundary-layer. 

The above description applies to incident shock-waves with strengths 
(xX) up to about 0-4. For these strengths there is no appreciable 
thickening of the boundary-layer upstream of the shock-wave and no 
separation. When the strength of the shock-wave is more than about 0-4 
however (figs. 6 to 9, Pls. IT. & IIT.), the boundary-layer begins to thioken 
and separate on the upstream side of the incident wave and the reflected 
shock-wave moves up to this point. Thus the incident and reflected 
shock-waves cross at a point outside the boundary-layer. As the 
incident shock-wave becomes stronger the foot of the reflected wave 
moves further upstream and the point of intersection moves further away 
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from the boundary. An expansion still occurs after the reflected shock- 
wave, but this remains near the point where the incident shock-wave 
meets the boundary-layer. In this region the rate of growth of the 
boundary-layer suddenly decreases, and there is a corresponding 
deflection of the stream due to the expansion. In figs. 6 and 7 there 
is a dark region after the expansion, indicating a further gradual 
compression similar to that occurring with the weaker shock-waves. 
The concave shape of the outer edge of the boundary-layer in this region 
can be seen clearly in fig. 7. The boundary-layer must always have a 
concave shape where there is a compression region, but it is not easily 
seen in the other photographs. The second compression would be 
expected to occur also for the stronger shock-waves in figs. 8 and 9 but 
in these photographs it is almost outside the field of view. (The dark 
region representing the compression after the expansion in figs. 6 and 7 
should not be confused with the dark expansion region caused by the back 
of the wedge.) 
Fig. 11. 


DIRECTION OF FLOW D 


Reflection of oblique shock-wave of medium strength (0:5<X<0-9). 


The dark region near the tunnel wall in figs. 7 to 9 shows that the 
boundary-layer has separated near the foot of the reflected shock-wave. 
This dark region is just visible in fig. 6 (x=0-48), suggesting that separation 
of the boundary-layer occurs if x is more than about 0-45. Figs. 6 and 7 
show a re-attachment of the boundary-layer further downstream. This 
probably occurs in all cases, but in figs. 8 and 9 the point of re-attachment 
is i field of view. RY! 

i. ney eee showing the typical configuration for Pn 
shock-waves of strength (x) between 0-5 and 0-9. AB is ee a ont 
shock-wave entering the boundary-layer at B. The eine ary- ahs 
separates and begins to thicken at C, and the reflected ure Tae b 
passes through this point. At a short distance Sd ge : eae ere 
is an expansion region E where the rate of growth of vas pe eke 
decreases rapidly through zero. After the expansion the ‘iS er = ge Oo 
the boundary-layer becomes concave and this causes a gradual compression 
in the region F, 
Fe bz 
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In fig. 9 the strength of the incident shock-wave is so great that its 
interaction with the reflected wave is of Mach type, with a short normal 
shock-wave between two intersections. Experiments by Fage and 
Sargent (1947) have shown that a strong normal shock-wave has a 
bifurcated foot where it enters a turbulent boundary-layer. It is 
suggested that the bifurcation at the foot of the Mach shock-wave in 
fig. 9 is similar to that of any normal shock-wave interacting with a 
turbulent boundary-layer at this Mach number, since the conditions 
at the upper end of the Mach shock-wave should not have much effect 
on the bifurcation at the lower end. For example, a similar configuration 
was observed with the wedge set at an incidence so high that the tunnel 
was choked and the normal starting shock-wave remained in the working 
section (fig. 10). (For this photograph a short-duration flash was used, 
as the flow was not quite steady.) The bifurcation of the normal shock- 
wave in this case is similar to that shown at the foot of the Mach shock-wave 
in fig. 9. 

This means that for any oblique shock-wave with x >1*2 the bifurcated 
foot of the Mach shock-wave and the separating boundary-layer would 
be generally similar to figs. 9 and 10. 

[t is perhaps of interest to point out the vortex sheets that arise at the — 
intersections of the shock-waves in figs. 8, 9 and 10. 

At lower Mach numbers than in these experiments (7. e. M<1-965) 
Mach reflection would occur for weaker incident shock-waves. Thus the 
upper limit of shock-wave strength giving a cross formation similar to 
fig. 11 (without a Mach shock-wave) would decrease with Mach number. 
At the same time the lower limit of shock-wave strength giving this 
configuration would probably not change much with Mach number, since 
it depends on the separation or thickening of the boundary-layer caused 
by the shock-wave. Thus the range of shock-wave strengths giving a 
cross formation, without a Mach shock-wave, would become more 
restricted with decreasing Mach number. This means that for sufficiently 
low Mach numbers this formation would not occur for a shock-wave of 
any strength. Weak shock-waves would give reflections like figs. 4 and 5 
and stronger ones would give Mach reflection, either with or without a 
bifureated foot. This probably explains why Fage and Sargent (1947) 
did not find configurations like figs. 6, 7 and 8. 
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. SUMMARY. 

A theoretical formula is deduced for the broadening and shift of 
spectral lines due to collisions ; this formula is valid in the microwave 
region, for gases at low pressures. It is assumed that the partial collision 
amplitudes f for the bimolecular collision problem are known, and using 
these, the equation of motion of the two-molecule density matrix is 
found. A formula for the line broadening and shift follows from this. 

Using the equation of motion of the density matrix, it is also shown 
that a Boltzmann distribution is the only equilibrium one. 


§1. INTRODUCTION. 


THERE are three causes of the broadening of the spectral lines emitted 
by a gaseous radiator, namely, radiation damping, Doppler broadening, 
and effects due to intermolecular interactions. These effects are very 
complicated, and it is only at low pressures, where bimolecular collisions 
are the dominant interaction mechanism, that a detailed theoretical 
treatment is feasible. The low-pressure effect is called collision 
broadening. 

Reducing the gas pressure naturally reduces the pressure broadening, 
and in the optical region, this effect is obscured by Doppler broadening. 
However, Doppler broadening is proportional to the radiated frequency a, 
while collision broadening is independent of w; hence, at sufficiently 
low frequencies, collision broadening must predominate. It turns out 
that this is the case in the microwave region (A~ 1 cm.). 

Since the war, a number of accurate experiments have been done on 
collision broadening in the microwave region—notably those of Bleaney 
and Penrose (1947) on the inversion spectrum of NH,—and this has 
stimulated theoretical work on the subject. 

Karplus and Schwinger (1948) have taken the important step of 
transcribing the classical theories of Van Vleck and Weisskopf (1945) 
and Frohlich (1946) into quantum form ; they are then able to account 
for saturation. They show that the correct tool for this type of calculation 
is von Neumann’s density matrix; its equation of motion may be 
written 
Oe 
Ere a Og gaia ee (1.1) 


a (HU Un? 
aeons as 
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where H?® is the internal Hamiltonian of one molecule, and 0,/dt denotes 
the rate of change of an operator due to collisions. Their assumption 1s 
equivalent to setting 


¢— ——o(U—U,), . . +--+. > (1.2) 


where U, is the equilibrium statistical matrix (we use the terms 
“ density ” or “ statistical ’? matrix indifferently). The original collision- 
broadening theory of Lorentz was equivalent to setting 


U,=const. exp (—H/kT). 


Van Vleck and Weisskopf (1945) pointed out that in the microwave 
region, the oscillating electric field applied to the gas (such measurements 
are always made in absorption) was essentially “ slow,” 7.e., that the 
collisions would attempt to maintain an instantaneous Boltzmann 
distribution. 


U,.= const. exp {—(H°—u..F cos vt)/kT } 


v being the frequency of the applied field. Owing to the rapid variation — 
of this field, the actual statistical matrix will be different from U,. This 
assumption gives a line shape in agreement with experiment. 

However, this phenomenological approach does not give any theoretical 
connection between the line width o and the intermolecular potential V(r). 
A theory due to Anderson (1949) has had considerable success in this 
direction, but it contains a number of assumptions which can be avoided. 
The theory presented below is restricted only by fundamental physical 
assumptions, namely that the gas pressure is small, and the applied field 
slow, in the sense explained above : this is believed to be true in the 
microwave region. 


§ 2. CALCULATION OF 0,U/dt. 


Karplus and Schwinger’s equation (1.2) is certainly too restrictive, 
but for a U not far removed from equilibrium, the equation 


U 
Panioeatn =2'L(ab, cd) (e|U—U,|d), 9 te 
x at ed 2 
must hold, since the definition of equilibrium is 
0.U, —() 
a 


The. object of this section is to calculate the coefficients L(ab, cd) in terms 
of known quantities. 

We have explained that we are taking bimolecular collisions to be the 
only mechanism producing broadening, and we therefore want to 
investigate the effect of bimolecular collisions on §, the statistical matrix 
of an individual molecule (U is an averaged statistical matrix). We shall 
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have to use a generalized statistical matrix D describing the external 
as well as the internal, motions of the molecule. The actanrial poenditres 
are continuous, and hence, in a coordinate representation, D will be a 
matrix in the internal coordinates, with each element a function of two 


external vector variables r and r’ (such a representation was first discussed 
by Born and Green, 1947). 


Sis very simply related to D, being its external trace ; 

Def. f 40th rds. ti (Dieter). dr... (2:2) 
the integral being taken over unit volume. As a simple example, the 
matrix 

D=S exp [tk. (r—r’)], 
describes a plane wave in a specified internal condition S. 
To discuss bimolecular collisions, we must work in a 2-molecule Hilbert 
space, which is the direct product of two one-molecule spaces. 
=H KH?” 


As long as there is no correlation between molecules about to collide, 
which is certainly true in dilute gases, D, the density matrix in #9, is 
given before collision, by 


aa 1 2 
D,=D® .D®, 


where the symbol a denotes the Kronecker product of two matrices. 
D, has to describe the external motion of two molecules, and will 
therefore be a function of four vector variables, D,(r,, f2 ; fj, ¥3)- 
S,, the two-molecule internal density matrix, is obtained from D, by 
a process analogous to (2.2) 


S.—=1fD3(t;, ry 3 ri, rp) . dr, . ie ee . . . (2.2 () 
S,, the one-molecule internal density matrix, is then found from 
SM=Tr,(8,). ee eect (22,0) 


Tr, denoting the operation of taking the trace over the internal 
coordinates of molecule 2. 

We shall suppose that we have, as a basis for #, a complete normal 
orthogonal sequence (¢.n.0.s.) {ly)} of eigenvectors of H®, the internal 
Hamiltonian of one molecule, corresponding to energies {E,,}. A suitable 
basis for 7, is then the c.n.o.s. 


{Jay}= {ly 125 


whose energies are {E,}, E,=E,+H, 
The Hamiltonian of a system of two interacting molecules is 


Ligkae ty 
H,=H(1)+H(2)+ ay (Pit Pa TY: 
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where H°(1) and H%(2) are the internal Hamiltonians of the molecules (1) 
and (2), p, and p, their momenta, and V the interaction potential between 
them. Since we are assuming the electric field to be slow, we do not 
need to. consider it while discussing collisions. We set 


H°(1)+-H°(2)= Hb, 
the two-molecule internal Hamiltonian, which satisfies 
H?|a)=E, |). 
We also set 
r=r,—r, , R=}(r,+r,), 


and similarly with rj, r5; then 


2 
H,=HI— = (v44VQ+Vie sD 
The equation of motion of D, is 
th Be =H,D,—D,H,, 


and with H, in the form (2.3), this can be separated. We write 


H,=H,— ee V2, where Hi;=H?— - V?+ V(r) 
and 
. DireR rR =f, £3.01 KR. RR): 
p and J must then satisfy 


30 
in = Hip pH, asta a 8. 5 
vets A fre 5s p 

and ha = re Ve)d(R, R'), > oS. eee b) 


p is therefore a matrix in the internal variables, but J is a function of 
R and R’ only. A particular solution of (2.4 6) is 


J(R, R’)=exp {ix.(R—R’)}, 


representing uniform motion of the centre of gravity of the two-molecule 
system. A more general solution of this type is 


J(R, R’)=J f(x) exp {ix.(R—R’). d3x, 


where f(x) is a normalized function giving the distribution of velocities 
of the centres of gravity of the two-molecule systems. At equilibrium, 
f(x) will be a Maxwellian distribution, since the oscillating electric field 
will not influence the external motion of the molecules. 
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We cannot solve (2.4 a) in the same way, since as we shall explain 
shortly, V(r) must be taken to be time-dependent. S, is then given by 


So=Jfe(r, r)J(R, R). dr d'R, 
(since d®r, . d’r,—d*r . d*R) 


=fJo(r, r) f(x). dr. d?x=Jo(r, r). dr. 
The coordinate R does not, therefore, enter into the problem, and we 
have only to find a suitable solution of (2.4 a). 

Now if we were to solve the ordinary static collision problem, we 
should only find the effect of one collision on a molecule. To attempt 
to find the effect of repeated collisions using a time-independent V would 
mean solving a many-body dynamical problem. 

This can be avoided by solving a two-molecule time-dependent 
problem (it is our assumption that binary collisions are the only important 
interaction mechanism that makes this simplification possible). We 
take the potential to be given by 


V(r, t)=V(r) if T—7r<t<T (T>0), 
=0 otherwise, 
V(r) being the given time-independent interaction potential. The initial 
condition is that at time t=0 


p(t, r’)—=[S.g(k) exp {ik.(r—2')}. dk, 


corresponding to a superposition of plane waves ; g(k) is a velocity 
distribution function which we shall later be able to show to be 
~ Maxwellian. We shall show how the solution of this problem can be used 
when we have obtained it. 

We shall assume that the asymptotic solution of the time-independent 
two-molecule collision problem is known in wavefunction form. In 
detail, we shall know that if the incident wave is given by 

(a, r|1)= [exp (tk. 1r)lam> - + 2+ + + (2.5 a) 


(representing a wave with momentum hk incident on a scatterer in 
state |m)), then the asymptotic scattered wave will be given by 


(a, x[8)— =P Mae!) (| fe, ua) |), © + + (2.58) 
where | 


12, =k? + 5 (Ep—E,). pte 8, WER Be (26) 


The coefficients f all being known. The complete asymptotic solution 
of the collision problem is then 


[C= [lee es ee ee meee 
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All states occurring in formula (2.5) are two-molecule states. _ The 
ket|C) is usually thought of as being given by Schrédinger’s equation— 
a differential equation—but it can equally well be regarded as being 
determined by the integral equation 


k 


Cinch ae 5 [SRP Mina?) a] V(s)|b)<0, 810). Ps, _ (2.7) 
b 


a 
 Agh? o 


where o=r—s. 


The determination of the f’s in the case of inelastic scattering is always 
very difficult, and it is made doubly so in these collision-broadening 
problems by the fact that we are usually dealing with non-central 
potentials. The Born approximation is not. satisfactory, since the 
thermal motions with which we are concerned are quasi-classical (see 
Anderson, 1949). This problem is quite distinct from the statistical 
problem now under discussion, but it is hoped to take it up in a later 
paper. 

The scattered wave (2.56) has the form of a spherical wave in free 
space ; this form can therefore only be valid at large distances from the 
scattering centre, where the interaction potential is inappreciable (this 
is what we mean when we call it an asymptotic solution). A solution 
of this kind will only tell us about that part of the line broadening due to 
‘permanent ’’ changes in the molecules produced by collisions ; this 
type is referred to as “impact” broadening, and will plainly depend 
on the intermolecular velocity. 

The other type, called ‘‘ statistical’ broadening, is due to transient — 
changes in the molecules, which disappear as they move apart. Statistical 
broadening may be pictured as being due to changes in the energy levels . 
of molecule (1) due to the influence of molecule (2). This will alter the 
radiated frequency, and the aggregate effect is to broaden the lines of 
the spectrum. These two types of broadening have been carefully 
discussed by Margenau and Watson (1936). 

The actual broadening will be due to a mixture of the two types, 
but it is generally thought that the collision broadening of microwave 
lines is mainly of the impact type : this is because collisions only last for 
a fraction of a cycle of the applied electric field. Recently, Howard and 
Smith (1950) have measured the temperature dependence of a prominent 
line in the inversion spectrum of NH, ; their measurements support this 
view. 

We must now consider the solution of (2.4). It may be written in 
the form 


p(t)=Ri(t)p(O)R*(), 


where the unitary operator R satisfies 


SUA a I kc) 
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Because of the form of p(o), the equation for p(t) is equivalent to 


(a, | p(t) [b, s)=2 |o(k) Pk. (a[ f(x, #)[¢)<c|S,(O) |d) (4 | o*(s, #)[6), 


(2.9) 
where ¢ is defined by 
{a| d(x, t)[6)=Ja, r| RM, s) exp (ik. s) . ds. 
We shall now show that ¢ can be simply expressed in terms of the /’s. 
We write 
Hj{=H)+ Vir, t), 
where H,, the “ non-interaction ” part of H,’, is defined by 


1 
H, = H$+ mee 


Later, we shall need the symbol H, to denote H°—yw.F cos vt: this sort 
of duplication is inevitable if we wish to use “H” for all the 
Hamiltonians occurring in the problem. 

(2.8) is equivalent to the integral equation 


1 ft 
R= oxp(Hyt/it) + = { __ exp (Hy(t—1 60} REP) ae 


From this, we deduce that ¢ satisfies 


h?k 
{a|d¢(r)|m)= [exp a ir)] 103, 6xD sae: 


74h KG rjexp {H)(t—t’)/ih} |b, s) d’s 
+32 7, 


x <b] Vis) |e)<c|4(s) |) dt. 
The interpretation of the representative — 
(a, rjexp {H,(t—t’)/ih|5, s)} 
is not immediately obvious, but, by transforming from a momentum. 
representation, we find that it is equal to 


in(e) = (&e-=™) 
(Qnh) (Z) XP 1" \ aah? a 


where x—=(t—t’)/m# and c=r—s. Now for t<T—7, ¢|m) varies with 
t according to 
eta a h?k?/myt ; 
ah 

i it i _vector corresponding to a state of energy 
ea) ae ie pate we are cece is PR e: so 
that its energy will be the same for t>T—r. Hence, if we define # by 

(Ey thik" +hk*, [mye 
(al b(e, t) [m)—=(a" v(e,t)|m) exp 
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will not have any oscillatory dependence on ¢. It will satisfy the 
integral equation ~ 


{a] o(r)|m)=5 exp (tk. r+ Facmp® | lis) [as 


x f° (Z)" fox i (Ga +A Hhee) } M1 8) Lm) 


where 1 ree and 2%.= 


The a-integral may be evaluated by the method of stationary phases 
(see Jeffreys and Jeffreys, 1946, p. 474). The asymptotic form of the 
equation turns out to be 


{a | (1) [m= 9), n XP (ik. rn) ce Qin) | Cxn ie 


x a] V(s)|5)<b | H(s)|m) as, 


2k, a= T) 27k (t -T-+7) 


where Q(r iss 1¢ = 


SIS 


=0 otherwise 
Hence, by comparison with (2.7), we infer that 


exp (2 


(a| P(r) |” )=Sam - Oxp (ik ° r)+Q(r (ees (a| fk, k,,.)|™). 
Hence in this problem, the “ scattered wave ” (represented by the second 
term) is the same as in the time-independent problem inside a shell of 
radii 

2hKkimg(t—T) 2hk,.g(t—T+7) 

m m 

and zero outside this shell. This shell is just the region which is 
classically accessible to the molecules, so that this equation is what we 
would expect on classical grounds. Using it, we can calculate p(t), and 
hence §,(t) by integration over the external coordinates. We find, 
asymptotically, that 


(a|S4(t) |b) exp (—iergaf){(a|8y(0) [0)+72K (ab, ed)<e]8,4(0) | a}. 


(2.10) 
where fiiw,,—EK,—E,, and 


K (ab, ed) — |g (ke) dh | | Calf (Kk, keam)|c)<b| f(s, b,,m)|d) dn 


ah ———_________ 
+ Fp, (aU, Le)Sua~ TFT HITT} ud plea 
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n is the unit vector in the direction of r, and dn denotes integration 
over solid angle, 7.¢., over all directions of scattering ; g(k) is the 
velocity distribution function. There are energy restrictions on the states 
occurring in the various terms of K(ab, cd) ; these are 


W_,p—,, for the first term 
ac 9 i second ,, 
cf LAbEGOl pe 
The first term represents the true scattered wave, while the second and 
third terms represent ‘shadow scattering” or loss from the incident 
wave due to the scattering process (we shall later justify this identification 
by classical analogy). 

In the absence of the interaction V, the time development of the 
density matrix is given by 

{a | Z,(t) | b )=exp (—tw,, pt) (a | S_(O) | b» 

(the symbol Z will always be used in this sense in future). Equation (2.10) 
can then be written 


{a|S8,(t) [b= a] Z(t) |b)+72 K(ab, cd) c|Z,(t)|d), . (2.10 a) 


Wpq=0 


which represents the effect of the interaction as a term proportional to 7, 
the time during which it is effective. 

In future, the symbols [a),|b),|¢)>, and |d) will be reserved for 
two-molecule states, and all other symbols will refer to one-molecule 
states. The quantity we are really interested in is not S,, but 8, : now 
S,=—S,S®, or, in representative form 


(a|S215)=(g 180 |2><v |S 12). | 
Taking the trace of (2.10 a) over the coordinates of molecule (2), we get 


(g|SPH | A= <g| ZV | &> 
tr ZX K(gy, hy; tj, wx) |Z O| wg |Z |). 


j, WX, Y 


(2.12) 


From a quantum-mechanical view-point, the molecule (1) must be 
regarded as interacting with all the (N—1) other molecules in unit 
volume. Thus equation (2.12) is incomplete, and should actually read 


g [SPO [y= glZPO1h> 
+7 SS Kiygy, hy; yj, we) [ZO [w) Gj ZP0 | 2). 
r=2 1, wx,y 
The matrix U of §1 is defined by 
N 
NU= 2 Sf”, 
f= 
and similarly we may define ‘ 
NW= 2 Z?. 


r= 
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Since N is very large, (2.12) can be written, with negligible error, 


Cg | SPH) [> = Kg [ZH [2 
4+Nr 2 K(gy, hy; i, wx) | Z(H) [w) (i |W) [@). 


ij, Wx, Y 


The same equation holds for any other molecule in the gas, 7. ¢. 


(g[SP(t) [b> = g[ZPO)[ 2) 
+Nr SF K(gy, hy; ij, wx) [Z| w)<j| Wz). 


tj, WL, Y 


Summing this over all 7, we have finally 


(g|U() [> = g | WO|h> 
+Nr 2 K(gy, hy; ij, wx) i] W()[w><g]W)|2), 


tj, Wx, Y 


and in particular 


(g| U(T)[h) = (g| W(T) |b) 
+Nr 2X K(gy, hy; vy, wx) | W(T)|w> gj] W(L)|x). 
pina e Pe ler ts) 


This formula is only true as long as the interval (T—7, T) is short, since 
it does not allow for the possibility of a molecule making more than 
one collision. But we are interested in the behaviour of U(t) over 
relatively long intervals of time ; since we are assuming the applied field 
to be slow even one period is, for our purpose, a long time. 

The second term on the right-hand side of (2.13) represents the 
redistribution of molecules among the various states by the collision 
process ; this rearrangement is going on continually. In the problem 
we have, as it were, “ frozen” the gas in the distribution U(O) up to 
time t=T—vr, and have then found the effect of single collisions on this 
distribution. 

Now up to the time ‘=T—v7z, when the interaction is switched on, 
U(t)= W(t), so that we can subtract U(T—7)=W(T—7) from both sides 
of (2.13), to give 


<g| U(T)—U(T— 7) |h) =<g| W(T)—W(T—7)|h) 
+Nr 2 K(gy, hy; ij, wx) <i] W(T)[w)<j|W(L) |x) 


4), WX, Y 
This shows clearly that (2.13) represents an infinitesimal portion of the 
redistribution process, and must be regarded as a differential equation 


dU 
xy a= <a h>+N a K(gy, hy; VY), WH )i| Ww) gj] Wa), 


ij, wx, y 


oW 
“Ot 
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since it is only true for every short time-intervals. Now from the definition 


of V 


ow. I HOW 0 
so that 
wy CUP ol 
C9 |r |b > == HW WHA) 


+N 2 K(gy, hy; tj, wx) <i] W]w) j[ We). 
49, WH, Y 
Since the interval 7 is so small, the difference between U and V on the 
right-hand side of this equation may be neglected, and it may finally 
be written 
0U 


. 1 
=-\|h >= 37, <9 | H°U—UB®|h) 


9 Ot 


+N 2 Ky, hy; tj, wx) |U[w)<g[U[x). 
ij, WX, Y : 

This interpretation of (2.13) as a differential equation is most easily 
grasped by thinking of the redistribution process as occurring in 
successive “‘ rounds.” If the initial density matrix is V—Up, this will 
be reduced to U, by the first round of collisions. For the second round, 
we must then write V—U,, and obtain as a solution U=U,, and so on. 
Plainly, these rounds must be of infinitesimal length since, in a finite 
interval, there is always a finite possibility of a molecule making more 
than one collision. -By comparison with (1.1), we see that 

we he-N ZS K(gy, hy; ij, wx) i]U]w)<g[U|z>.  . . (2.14) 


tj, Wx, 


<9 


The equilibrium statistical matrix U, is defined by 0,U,/0t=0, 7. €. 
ES K(gy, hy; ij, wx) [U..w)<j|U.|%=0. 
aj, Wx, Y 
We are interested in the value of 0,U/dt when U is fairly near equilibrium, 
2. €&., when 
Y=U-U, 

is small compared to U,. Substituting this into (2.14) and neglecting 
the term of order Y?, we get 


a,U 
CIF 


h>=N  K(gy, hy; 1j, wa) [i |U,[w)<g][¥ |x) 


1j, WX, Y 
+] ¥[w)<j]|U.1*>] 

—N 5 [K(gy, hy; ij, we) +KGgy, hy; jt, x) 
ij, WX, Y 


x G[U.[w)>G1Y[*): 
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on rearranging the dummy indices on the second term. Comparing 
this last formula with (2.1), we see that 


L(gh, wx)=N 2 [K(gy, hy; iw, je) +K(gy, hy; wi, xj)]x @|U.[ 7): 
i,5,Y ¢ 
A i 


the specific formula for the coefficients L which we wanted to find. 


§3. THe Derivation oF U,. 
We should expect U, to be a Boltzmann distribution 


exp ( a % K,/kT) 


Yexp (—HJeT) voy 


G[U.|j)= 


but we have to show that this matrix satisfies the equation 0,U,/dt—0. 
We shall do this by showing that the solution of 


EK (ab, cd)Xc|Sy|d)=0 
cd 


st vie _ exp (—E,/kT) 
1S RA LD A Sarees (—E,/kT) ed> 
- 
for the validity of (3.1) then follows by taking the trace of this equation 
with respect to the coordinates of molecule (2). 

To justify this last equation, we need to prove a reciprocity theorem, 
analogous to the Lorentz theorem of electrodynamics, which states 


<b | f (ka, —k,,B)|a) = <a lf (kp _B,—ka) |b). . . | . (3.2) 


This is easy enough to prove for elastic scattering, but for inelastic 
scattering it becomes somewhat complicated, owing to the large number 
of coordinates involved. In order to keep the writing to a minimum, 
we shall only prove the theorem for the simplest possible case, that of 
the scattering of electrons by a hydrogen atom. The generalization of 
this to more complex cases is obvious enough, if rather tedious. 

We take the proton as the centre of coordinates (assuming it to be so 
massive compared to the electrons that it does not move), and let the 
coordinates of the bound and free electrons be r, and r, respectively. 
We expand the overall wavefunction #4(t,, >, ¢) in terms of the 
eigenfunctions {¢,,(r,)} of the scatterer 


es r,, t) Bale, t) 


We also take another solution 7,(t) of the equation of motion, and 
expand it in the form 


Po(— =A bil r,, 1)G,,(t2, —t) 
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[since ¢*(t) and 4,,(—t) have the same dependence on ¢]. We can thus 
show that F and G must satisfy the equation. 


PEPAObn(t) {iF 5 [E(—OF (0) © 1G, VR, —, v4,]} =0, 


the argument being much the same as that used to deduce the wave- 
mechanical current vector. We now integrate this over d°r,, and since 
the {¢,,} form a ¢.n.o.s., we get 


h2 


42 
Bfin 54 1Gu(—OF,()1— =—[ 


es ot Gr ya ve vCo) =0, 


If 4, and %, correspond to states of the same energy, the time-variations 
of F,, and G,, are equal and opposite, so that the first term is zero ; also 
the time factors cancel in the second term, giving 

2(G, VF, —F,, V*G,)=0 

n 


or, by Green’s theorem 
Z| (evG 2G Vi erve= 0 


S being any closed surface. This equation is true for any inelastic 
scattering process. 
In agreement with (2.5), we take F,, and G,, to be given asymptotically 
by 
F,(r)=[exp (ihr. 0) 8,4 PO”) (nf (hee, nat) |) 
exp (tk, yp?) 


/ 


G,,(r)=[exp (th pM- B)J8,,+ (2 boaBs Frat) |6)- 

These refer to states of the same energy ; F represents a wave incident 
in direction a on the scatterer in state |@), while G represents the scatterer 
in state |) and the incident wave coming in direction 8. Substituting 
these into (3.3), and performing the integrations, the theorem (3.2) 
follows at once. 

Another result which we need follows from the fact that the operator R 

introduced in equation (2.8) is unitary. This gives us the equation 


B | dale) [oy dx Cel $(2)]2)—Pun 
and using our equation for ¢, this reduces to 


2g | CFU Fees) (CL, 6) A= — Aeron {Calf 1D): | 


(3.2) and (3.4) represent basic 
ion problem. (3.2) tells us that 
is the same as the cross-section 


(subject to the restriction Wp =): ’ 
restrictions on any set of f’s for a collis 
the cross-section for any given process 


SER. 7, VOL. 42, NO. 324.—JAN. 1951 A 
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for the inverse process, while (3.4) is a consequence of the conservation 
of mass ; any particles scattered must be compensated by losses from 
the incident “‘ beam.” 

By repeated use of the reciprocity theorem, we can show that (3.4) is 
equivalent to 


> | Keo (By) C2. | (a|f(k, beam) ]e) <b] Fk, &,,m) |e). dn 


= —4nIm{{g(k) d®#k(a|f(k,k)|b>}, . . . . . (34a) 


g(k) being the velocity distribution function. 
Our object is to find a value of S, which makes 


(a]C|b)X K(ab, ed) <c| S| d) 
ed 


zero, © being proportional to the two-molecule scattered wave. Now 
we know, on general statistical-mechanical grounds, that S, must be 
diagonal in an energy representation, so that we need only consider the 
diagonal elements K(aa,cc) of the array K(cd,ab) (this can be proved 
directly, using the reciprocity theorem). We then have 


{a|C]a)= XK(aa, cc) {e|S,|e) 


m 


=2{| 2hkna (a | f(k, Kal) |e) | 2dn. 2k 


x {9(k)<c] S.| ¢)—g(k,,) <a] Sg] a>}. . . (8.5) 


from (3.4 a). (3.5) can be compared with the corresponding classical 
equation, for (c| S,| ¢) can be identified with p,, the fraction of molecules 
in the state |c), while 


k 
ea} a] f(le, knat) |e) 


is o,,(m), the differential cross-section for a collision in which the scatterer 
goes from state |c) to state |a@), the scattered particle going off in 
direction n. Remembering that 2hk/m is just v, the relative molecular 
velocity before collision, (3.5) may be written 


(alcja)= 2 | | V . Gea(M) {9(I)p,—9(K,,) py 2k. dn. 


This is identical with the classical expression obtained by direct physica] 
arguments, which confirms the accuracy of our analysis. The first term 
represents the increase in the population of state |a@) due to molecules 
being knocked into this state from other states; the second term 
represents losses in this population due to molecules being knocked out . 
of this state. This argument identifies the terms of (2.11). 
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Now an argument, due originally to Boltzmann shows that <a|C] a) 
is Zero if, and only if, 


h?2 \3/2 
g(k)= (<a) exp (—h?k2/mkT), 
and S.=A exp (—H9/kT) 


(A being a normalizing constant) ; ¢. e., if g(k) is a Maxwellian velocity 
distribution, and 8, is a Boltzmann density matrix. Hence we have 
shown that 
¢|U.] 3) =| ppl 4)9,, 
ore : .._ exp (—E,/kT) 
{t| pal) 2 exp (—E,/kT)’ 


(py, thus being properly normalized). The formula (2.15) for L(gh, wa) 
thus becomes 


Ligh, wt)=N2X[Kigy, hy; iw, 1x)-+K(gy, hy; wi, xt) <t|pp]?). (3.6) 
ay 


§4. DEDUCTION oF A CLOSED FORMULA FOR o. 


In future, we shall use the symbol o,, to denote the ‘“ complex 
breadth’ of the spectral line corresponding to the transition g—h. 
By this, we mean that o,,,is a complex number whose real and imaginary 
parts give the line breadth and shift respectively. We want to obtain 
an expression for o in terms of the L(gh, wz). 

The equation of motion of S, including the effect of the applied electric 
field, is 

ees DL (gh, we) (w|S—S 


where 
H,=H°—up. F cos vé. 


We solve this by making the weak-field approximation, putting 
S=S,+D.F 


- and ignoring terms in F?. In the microwave region, all energy differences 
are. <kT, so that 
H, u-F 


apt rT cos vt. 


S,~1 


This gives us the equation of motion 


oD 
“Ot 


p> = 250 gy |) — ing, g|D|)+ Z(Lgh, we) w|D |x). 


ie 
,I KT 
(4.1) 


52 D. C. M. Leslie on the 


Now rotational degeneracy must always be present. This is a con- 
sequence of the fact that, in an infinitesimal applied field such as we are 
considering here, all permitted orientations of the molecule have the same 
energy. We shall assume that there is no other form of degeneracy 
present ; this is usually so, and when it is not the theory given below can 
easily be extended. 

Hence any state must be labelled 


[9,7 
where m, is the magnetic quantum number ; E, is, of course, independent 
of m,. We shall assume that there is no accidental degeneracy, that is 
that w,;, cannot equal w,, unless |j)=|g) and |k)=|k). It seems 
likely that this assumption is correct for all actual atomic systems. 

Suppose now that the applied field is in the z-direction ; we shall then 

want to know D., and taking account of the orientational degeneracy, 
(4.1) may be written 

v sin vt 


Cee 
<q, m, | (5 +iey,) D;|h, m3)= aire <g, m,|u-|h, me) 


+ 2 Lig, m1; h, me: w, m3; 2, M4) Cw, my | D,| 2, m,). 
W, 7435 ©, M4 
Now since D. F is a scalar, D must transform like a vector ; that is to 
say that D and uw transform in the same way under a rotation, so that 


(9; My | ib} |r, Ms )—=<(g|A XCE my [u. |, Ms), 


where the representatives <g|A|h) are independent of the indices m. 
Then A is given by 


ae v sin vt 
<g, m,|u.|h, m2) {¢g | (5; iy.) rA|hy— ET \ 
= (wlA|x) LY Lig, my; h, mg: w, mg; x, m4) (w, ms [mz |, m4). 


W, 2% 3, 724 
(4.2) 


The general solution of this equation is extremely complicated. 
However, in the case where the breadth of the lines is much less than their 
separation, we can derive a simple formula for o. This is obviously the 
case of practical importance, since when the lines begin to overlap 
seriously, it will not be possible to measure their breadths. j 

The case not covered by this assumption is that of a very wide variation 
in the spacing of the lines. Thus, if two spectral lines were very close 
together, they might be merged while the spectrum was apparently well 
resolved. However, this line would probably have an anomalous shape, 
and if such an effect were found experimentally, it would be quite easy to 
check the above explanation. 

In mathematical form, this assumption is equivalent to 


Oi mn > L(gh, we), 


) 


Ay eed ie : ‘ =a 
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ae set of indices. We now suppose that v, the frequency of the applied 
a: e aaa. gy, in the sense that w,,—v—=A is of the-same order of 
agnitude as L. It is quite easy t6 show that, under these assumptions, 


“ja k» is of magnit BS: 
| gnitude EL aed (3 bq, le), 


while 

{g|A|h) is of magnitud Hovis 

|A|h» agnitude 745 - 

Hence <g|A|h) is much larger than the other representatives of 
A, and therefore in the sum on the right-hand side of (4.2), only 
the terms containing (g|A|/) need be retained. [This assumes that 
L(g, m,; h, Mea: J, Ms; h, m,) is not much smaller than the other coefficients 
ive but this would give rise to a very narrow line, which would plainly 
be anomalous]. (4.2) then becomes a simple equation for <g|A|/) 


Om. 
1 (5 + iy.) 9, mM, | w2 | A, Ms) 


eee Ly, my > h, Me: J, Ms; h, M4) <9; Ms [u.[4, m,)} 


M3, 14 


sin vt 
x dg [AW) [b= a gm jusliim,). 2 hs. (4.3) 


Since we know that <g[A|/) is independent of the indices m, we clear 
them out, and at the same time convert (4.3) from a vector to a scalar 
equation, by multiplying both sides by 


{h, Mz [te 9; M1 >= (9, My |u|, Ms; 
and summing over m, and my». This gives us 
igs | v sin vt 
Gi +i gn—o9n) (g|A|2)= aR > o oo Oo (4.4) 


where 


ZS L(g, mj h, Ma: 9, Mg; h, Ms)G, Ms | uz |b, ma Kh, Me [u21 9,1) 
mim2mn3ms 
XS |<g.my [H2| 2 ms) |? : 


myme2 


(4.5) 


the calculations of Karplus and Schwinger, we 
ds to a line of Van Vleck-Weisskopf shape, with a 

An equation substantially eqivalent to (4.5) has 
(1949), by a different method. 


By comparison with 
see that this correspon 
complex breadth o,,,. 
already been deduced by Anderson 
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We have equations giving L(gh, wx) in terms of K(ab, cd) and K(a6, ed) 
in terms of the f’s, so that we can obtain a formula for o,, in terms of the 
fs. The details of the substitution are quite trivial, though rather tedious ; 
the result is 

oa NI @on(keg(k) de VAL ° 4 Aiea 
of ~ IG: mM, | Hz | 2, Ms) |?’ 


where g(k) is the Maxwellian velocity distribution, and 


2h ; ‘ 
gh 2 <i | py] t><g, my | ue | A, m1 <h, ms [He | 9, Ms) 


M Li, mg; y. my 3 my mM 
x J{Ch, ms; ¥, My [Ff (Ik, bin) [h, 2,5 4, m;) 
XG, Mg; Y, My | f (Kk, kiym) |g, m4; %, m;) 
+ Ch, ms; y, my |f(, kim) |i, m,3 h, m1) 
dg, ys My Lf, azn) [is msg, ,)} da | 


Q 


dich : ; 4 
+ Ft 2 Glew diy dorms lel sms) F 


a, Ng 5 Ma 


(KA, m,; 7, m;|f(k, k) | h, m4; 7, m;) 


aay 9; ™4; t, mm; | f(k, k) l9; ™4; t, m)| 


abs 2 <g,m,|u,|h, m,)<h, ms] u,| 9, m3) 


MMs 


(<9 | ppl g)<h, mg; 9, m1 | f(K, k) |g, ms; h, m,) 


~ Chl pm | 0) <9, mgs, ma | F ) |, m3; a,-m,)) | f 
(4.6 a) 


(the last two f terms are not necessarily equal, since f need not be 
Hermitian). This will, in general, be a complex quantity, so that we 
would expect collisions to shift lines as well as to broaden them. It is 
quite easy to show that if E;<kT 


shift O K; 
breadth — (i) 
Hence in the microwave region, the shift will be a relatively small effect, 
and it is not surprising that it has not yet been detected. 

The equation (4.6) is believed to be more accurate than that used by 
Anderson. However, using Anderson’s equation, it is quite easy to 
‘calculate o,, from a known V(r), while using this equation, it is very 
difficult because of the difficulty of determining the f’s. It is therefore 
hoped to show, in a later paper, that Anderson’s equation is a first 
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approximation to our (4.6), and to estimate the error of this approximation 
(which will presumably be mainly due to neglect of inelastic collisions) 


T am much indebted to Professor M. H. L. Pryce for the advice and 


encouragement he has given me during the course of this work, and to the 
D.S.I.R. for a maintenance grant. 
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ABSTRACT 

The electrical resistivity of cold-rolled copper, low carbon steel, brass 
and sterling silver has been found to vary with direction of measurement. 
The magnitude of the anisotropy may be qualitatively accounted for by 
considering the effect of oriented dislocations. The temperature 
coefficients of resistivity of these rolled metals are isotropic to within 0-5 per 
cent. These results indicate that deviations from Matthiessen’s rule are 
to be expected. 


$1. INTRODUCTION. 


THE effect of cold-work on the electrical resistivity of cubic metals and 
alloys has been studied mainly after wire-drawing (Boas 1947). However, 
the possibility of development of anisotropy of electrical resistivity by 
* cold-work has generally been overlooked. Matsuda (1925), in a general 
investigation of the effect of cold-rolling on the physical properties of 
copper and aluminium and their alloys, found small but not systematic 
differences between the resistivities of longitudinal and transverse 
specimens. Becker and Born (1925) showed that anisotropy of electrical 
resistivity in cold-rolled tungsten could not be greater than 5 per cent. 

In the present paper experimental evidence is given for the existence of 
anisotropy in cold-rolled cubic metals, and an explanation of this effect is 
put forward. 


§2. EXPERIMENTAL. 


(i) Materials. 

The six cubic metals selected for study were aluminium of commercial 
purity, tough-pitch copper, fine silver, low carbon (0-07 per cent) steel, 
«+f brass, and sterling silver ; their analyses are given in Table I. 

Annealed strips of each metal, about 60 mm. wide and 2-5 mm. thick, 
were reduced 95 per cent in thickness (to 0-125 mm.) by cold-rolling. 
Specimens of about 60 x 2 mm. were cut from the strip at 0°, 45° and 90° to 
the rolling direction. Some of these specimens were then recrystallized 
by vacuum annealing. 


* Communicated by Dr. W. Boas. 
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(ii) Determination of Resistivity. 

A Kelvin Double Bridge was used for all measurements of electrical 
resistance and the conventional current and potential contacts were made 
using a system of four needle points on each specimen in turn. The 
resistance of the 20 mm. length between the potential needle contacts was 
measured while the specimens were immersed in a bath maintained at 
1:0--0:3° C. The usual precautions were taken to eliminate errors due to 
parasitic E.M.F.s. Some resistance measurements were also made at 0° C. 
in melting ice and at — 183° C. in boiling oxygen in order to. determine 
whether the temperature coefficients of resistivity were isotropic. 

For the calculation of resistivities the average cross-section of each 
specimen was found by determining its weight and length and by assuming 
a value for the density*. As the thickness of the strips from which the 
specimens were cut was not uniform across the width, the average cross- 
section found for specimens cut at 45° and 90° to the direction of rolling was 
not that of the central 20 mm. over which the resistance measurements were 
made. 

The difficulty in measuring the actual thickness of the central portion of 
specimens was avoided by remeasuring the resistance of the specimens of 
each orientation after they had received a suitable vacuum annealing 
treatment. The resistivity of annealed cubic metals is isotropic, so that 
differences in the results obtained from the annealed specimens of different 
orientations were attributed to the effect of non-uniformity of thickness. 
The value of resistivity found for specimens at 0° to the rolling direction 
was assumed to be free from this error. The test of anisotropy was whether 
the fractional change of resistivity on annealing varied with orientation. 
For this, a statistical treatment was developed. 


(iii) Statistical Treatment. 


In each of the three orientations, 10 specimens were taken and the 


resistivities 2,; (#=1, 2, 3; fea U2. cee is 10) were measured. After 
annealing the specimens were found to have resistivities %;; (=, Day; 
Ree lp2epp ict 10). Leta, q; a, (i=1, 2, 3) denote the mean resistivity for 


each set of ten specimens. We require to test whether the change in 
mean resistivity on annealing is dependent on orientation. 
Assuming that the variation within each set of 10 measurements is due 

solely to random fluctuations, we have, as a pooled estimate, S”, of the 
error variance, 

3 10 oe ae 

EPs (%;;—%i—)?+ 2 a (%,;—%;9)” 

Oe 4=1 j=1 t=1j=1 ' 
60—6 

a or te te SS 


ities i i istivi re: aluminium, 2°70 ; 
* The assumed densities in evaluation of resistivity were : @ am, 2°70 5 
copper, 8-94; silver, 10:60 ; 0-07 per cent carbon steel, 7:80; 64/36 brass, 8 5G; 


sterling silver, 10-31 gm./em.’. 
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The estimated error variance of the mean of 10 specimens is S*/10 and as 
the variance of the difference of two variates is the sum of the variances of 
both variates, the difference («;,—2;,) is seen to have an estimated error 
variance of 82/5. To compare the resistivity differences for any two 
orientations we therefore employ the standard error, \/2S. In the present 
case, referring to the t-table, we have ¢=2-01 for 54 degrees of freedom and 
probability P=0-05. Hence, if the resistivity differences at two orienta- 
tions differ by greater than 2-01 times the standard error, the effect of 
orientation is significant. 

It was found convenient to express the difference between means of 
cold-worked and annealed specimens of any one orientation as a fraction 
of the “annealed” resistivity in that direction and accordingly, the 
standard errors were calculated as fractions of the annealed resistivity. 


TABLE I. 
Analyses of Materials Used. 


Aluminium .... 0:5% Fe; 0:2% Si; traces of Mn, Mg, Cr and Cu. 
Copper 0-12% 0,(1:06% Cu,0) ; traces of Ag, Ni and Fe. 
Fine Silver Traces of Mg and Cd. 


Mild Steel 007% C ;0-21% Mn. 
Brass 63-5% Cu ; 36-59% Zn ; traces of Ag, Mn and Fe. 
Sterling Silver .. 92-8% Ag ; 7-2°% Cu; traces of Mg and Cd. 


§3. RESULTs. 


No changes in resistivity of the fine silver specimens were found after a 
‘“ high ’ temperature annealing treatment, and it was assumed that the 
cold-rolled strip had previously recrystallized at room temperature. No 
results for silver are therefore reported. 

The resistivities of the annealed metals used, determined at 1° C. on 
specimens at 0° to the rolling direction, are given in the second column of 
Table II. The fractional changes in resistivity, 4p/p, on annealing (and 
their standard errors) relative to the resistivities of the annealed specimens 
are also given in Table IT. for each orientation. The last column gives the 
anisotropy of resistivity of the 95 per cent cold-rolled metals used, 
calculated as 

1+(Apl/p)o _ (pplpa)o 

1+(Ap/p)oo — (pn/padao’ 
where pp, p, are the resistivities of deformed and annealed specimens, 
respectively. 

Table ITI. gives the values of the ratios of the resistances of specimens at 
—183° C. to those of the same specimens at 0°C. The ratios must be 
independent of orientation for the annealed specimens and the variation 
found gives an estimate of the experimental error. It will be seen that the 
ratios for the rolled materials are also independent of orientation within 
this error (~0-5 per cent). 


/ 
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§4. Discussion. 

The average values for the increase of resistivity accord well with pub- 
lished results, except in the case of copper. Table IT]. shows that for steel, 
brass and sterling silver, the transverse resistivity was significantly less 
than that in the direction of rolling. The reverse was found for copper. 


Taste II. 


Fractional Increase of Resistivity on Cold-Rolling (95 per cent Reduction 
in Thickness). 


Resistivity Ap/p x 10? (at 1° C.) 


Anisotropy 
of 
resistivity 


Orientation to rolling direction 


0° 45° 90° 


— - — - — 


1-58-+-0-41] 2-07-L0-41] 1-86-+-0-41 | 0-998+0-008 
4-03-+-0-46| 3-93+0-46] 5-04--0-47 | 0-990-+-0-009 
3°63-L0-40 | 2-32-40-40|] 2-65-+-0-40 | 1-010+-0-008 
25-13-+0-31 | 23-46+0-31 | 23-07-+-0-32 | 1-018-+-0-006 


20-19-+0-45 | 19-34-L0-44 | 17-88-L0-44 | 1-020+-0-008 


Taste III. 
Change of Resistance with Temperature. 


R_ isg/fo 


Condition Orientation to rolling 
direction 


0° 


Aluminium... . Annealed 
Rolled 
Annealed 
Rolled. 
Annealed 
Rolled 
Annealed. 
Rolled 

Sterling Silver.. | Annealed 
Rolled 


ay be attributed to the high oxide 
f the high value of the increase in 
lling direction were intermediate 


This anomalous behaviour of copper m 
content, which may also be the cause 0 
resistivity. The resistivities at 45° to the ro 
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within experimental error, between the longitudinal and transverse 
resistivities. No significant effect of orientation was found for aluminium, 
but it should be noted that the error was high, relative to the increase of 
resistivity, and an effect of the magnitude found for the other metals 
would have been undetected. 

The isotropy of the temperature coefficients of resistivity (Table ITI.) is, 
at first, surprising. However, Griineisen and Goens (1924) found that the 
temperature coefficients of zinc and cadmium single crystals were isotropic, 
despite the pronounced anisotropy of resistivity. 

The present results indicate that deviations from Matthiessen’s Rule 
can be expected for a worked cubic metal in which the resistivity is 
anisotropic. Rutter and Reekie (1950) have recently found deviations 
which they interpreted in terms of an increase, due to cold-work, in the 
thermal component of resistivity. It is suggested that a part, at least, of 
the deviations found by them was due to the dependence of resistivity on 
the direction of measurement. From considerations given below it will be 
seen that anisotropy of resistivity is likely to be more marked in drawn 
wires (as used by Rutter and Reekie) than in rolled metals. 

Three possible causes of the anisotropy of resistivity can be considered ; 
the effect of (i) a second phase, (ii) general lattice strains, and (iii) oriented 
dislocations. 

(i) The effect of a second phase. 

None of the metals used was single phase. However the anisotropy 
found appears to be independent of the estimated quantity of second 
phase (compare ¢.g. aluminium and steel). Again, the anisotropy 
developed in brass and sterling silver was similar, although in the case of 
brass, the second phase has higher hardness and lower resistivity than the 
matrix and in the case of sterling silver these properties are very similar for 
both phases. From this indirect evidence it is probable that the presence 
of second phases had little influence on the anisotropy of the metals used *. 

(ii) The effect of general lattice strains. 

Bridgman (1925) has studied the effect of tension, within the elastic 
range, on the transverse and longitudinal resistivities of annealed, flat 
specimens of several metals. The development of anisotropy was shown 
and found to vary approximately linearly with the applied load. 

This behaviour was presumably due to deformation of the cubic lattice 
under the externally applied stress. In cold-rolled strip the sum of the 
internal stresses must be zero, and any effect depending linearly on the 
stress is therefore also zero. This condition holds approximately in the 
strip with a high degree of preferred orientation of the crystals. However, 
if the distribution of crystal orientations is such that the sum of the 
strains is not zero, owing to the anisotropy of elastic constants, an effect 
on the electrical resistivity may be expected. 


* 'These considerations apply only when no cracks develop between the second 
phase and the matrix. It is tentatively suggested that in the present experiments 
the copper, which contained 1-1 per cent Cu,O0, had a higher transverse 


resistivity than longitudinal owing to the formation of longitudinal cracks. 
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A gross overestimate of the magnitude of this effect may be set by 
considering the application of a stress equal to the yield stress of the 
material used. Foriron, Bridgman gives the values of +1-42 x 10-8 for the 
longitudinal and -+0-5410-° for the transverse fractional changes of 
resistivity per kilogram persq.cm. If we assume a value of 8000 kg /em/” 
for the yield stress, we have 1-007 for the maximum value of the anisotropy. 
For copper, assuming a value of 4000 kg./em.? for the yield stress, and using 
Bridgman’s values of +1-7510-* and —2-410-8 per kg /em.? for the 
longitudinal and transverse coefficients respectively, we find the maximum 
anisotropy to be 1-015. 

These values are comparable with those found in the present experiments, 
but it should be emphasized that the effect of internal stresses will be very 
much less than that of a uniaxial stress equal to the yield value. The effect 
of long-range internal stresses can therefore be regarded as unimportant. 

(iii) The effect of oriented dislocations. 

Recent calculations of Mackenzie and Sondheimer (1950) on the increase 
of resistivity produced by cold work lead to a prediction of anisotropy. 
They consider the scattering of conduction electrons by dislocations and 
show that the increase in resistance parallel to a dislocation axis is zero 
while the increase in the slip direction exceeds that normal to the slip 
plane by a factor of about 2, depending on Poisson’s ratio. 

During the rolling of polycrystalline metals, the dislocations introduced 
will not be oriented at random as slip will occur preferentially along 
favourably oriented slip directions. When cold-rolling has produced a 
preferred orientation, a complete calculation of the anisotropy would 
require an estimate of the relative abundance of dislocations perpendicular 
to each slip direction in each of the slip planes in the ideal preferred orienta- 
tion produced by cold-rolling. Such a calculation is very sensitive to the 
distribution of dislocations among the various slip planes and of crystal 
orientations so that an estimate is not given here. However, it is 
clear that the present results for anisotropy of increase of resistivity 
(Ap/p)o/(Ap/p)9o (brass 1-195, silver 1-130, steel 1-370) can be explained 
by the above considerations if appropriate assumptions (suggested by the 
mode of deformation) are made. 

Recent developments of dislocation theory by Cottrell and Bilby (1949) 
have led Bhatia (1949) to predict anisotropy of resistance change during 
strain ageing in iron single crystals. Interstitial solute atoms are assumed 
to form “atmospheres ” around non-randomly oriented dislocations owing 
to the presence of local lattice strains. In such a case, the effect of dis- 
locations alone, as considered above, would be enhanced. It is interesting 
to note that in the present experiments the highest value of anisotropy of 
resistivity increase (1:370), was found for steel. This suggests that the 
diffusion of carbon atoms to dislocations, as considered by Bhatia, had 
occurred and accentuated the purely dislocation effects. | 

From the above considerations it is concluded that the major cause 
anisotropy of resistivity of cold-rolled cubic metals is the presence 0 


oriented dislocations, 
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ABSTRACT. 


Estimates for the energy and momentum lost by particles which 
_ produce stars when passing through the nuclei in a photographic emulsion 
support the conclusions reached in a previous paper that the majority 
of disintegrations are not produced by neutrons which are captured in 
a single nuclear encounter. The ratio of these two quantities, which 
have been obtained from the detailed statistical analysis of the various 
particles produced in nuclear disintegrations, suggests that nuclei are 
partially transparent to high energy neutrons which at 3600 m. have 
sufficient energy to produce a small number of stars. 


INTRODUCTION. 


In a previous paper (Harding 1949 a) the nature of the radiation 
responsible for the production of cosmic ray stars has been discussed 
by considering the momentum and energy transfer from the incident 
particle to the excited nucleus. An estimate for the transfer of momentum 
was obtained from the angular distribution of the residual nuclei, or recoil 
fragments, which are left when the original nuclei have cooled down by 
particle emission. However, this analysis was carried out for stars 
formed in C2 ‘“‘ Nuclear Research’ emulsion which is not capable of 
recording tracks left by protons of more than 50 MeV. It is therefore 
necessary to re-examine the former experimental data and to include 
further results obtained with electron sensitive plates in order to investigate 
the validity of the previous conclusion, viz., that the majority of nuclear 
disintegrations at 3600 m. are produced. by nucleons, but that these are 
not absorbed catastrophically in a single encounter. 


EXPERIMENTAL RESULTS. 


(a) The angular and range distributions of recoil fragments are given in 
the previous paper, as is the size distribution ofstarsexamined. In this paper 
the size of a star is denoted by the number of evaporation (black) tracks. 


* Communicated by Sir George Thomson,-F.R.S. 
+ Present address : Atomic Energy Research Establishment, Harwell, 
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(b) The variation of the «/p ratio with star size has been given by 
Perkins (1950) and Page (1950). The ratios assumed in this paper for 
stars in heavy elements (Ag, Br) are : 


Star-size 2, 3 and 4 5-14 15 
a/p ratio 0-4 0-5 0-45 


(c) The particles responsible for the “ grey ” tracks in these stars are 


mainly protons of energies 30-350 MeV. (Harding 1950). In plates’ 


exposed at high altitudes, Camerini et al. (1950) have shown that about 
5 per cent of such tracks are produced by mesons. In the following we 
assume that all ‘‘ grey’ tracks are produced by protons. (It would be 
easier to maintain our original conclusions were we to assume a certain 
percentage to be produced by mesons, as for a given momentum transfer 
the energy transfer would be greater than for protons.) 

From an analysis of these “ grey’ tracks (see Appendix) we obtain 
the following : 


Average energy (including binding energy) per particle=90 MeV. 
Average downward momentum per particle=110+20 MeV./c. 


(d) It is intended to apply momentum and energy consideration to 
the heavy nuclei, but in the previous paper a certain number of stars 
arising through the disintegration of the light elements (C, O, N) were 
also observed. In order to allow for such stars we recall an experiment 
(Harding 1949 b) designed to give information on these light stars, 
i.e. on stars formed in layers of gelatin. 

Having now observed 38 such stars, we estimate that 35-.6 per cent 
of a random sample of stars in ©.2 emulsion arise through the 
disintegration of the light elements. 

Assuming a cross-section for star production proportional to A232 
(which is approximately correct), we estimate that from these 38 stars 
we should expect a charge of 260e to be emitted. We have seen 

62 protons (deuterons or tritons): charge 62e, 
54 «-particles (some may be due to heavier fragments but 
identification uncertain) : charge >108e, 
3 particles identified as having Z>3: charge ~10e. 


Estimated charge due to tracks of length greater than 10 y, these tracks 
remaining in gelatin layers, 7. e. “lost”, charge ~10e. In the Appendix 
we show that in a random sample of small stars there are three-quarters 
‘“ grey’ proton per star. We interpret these as ** knock-on’ protons, and 
they should therefore be produced in light nuclei, from which they should 
stand more chance of escaping than from heavy nuclei. We are therefore 
justified in assuming that at least three-quarters ‘“ grey’ proton are 
produced per light star, 7. e. emitted charge >? . 38e—29e. 

Therefore total emitted charge is >219e. 

Therefore missing charge per star is <41le 38, 
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This must all go into tracks of <10,, and these must remain in the 
gelatin layer. It seems unreasonable, and is certainly not expected 
from evaporation theory, to suppose that we get one short proton per 
light star. Rather we asswme that one-third of these stars have a 
“recoil” of average charge 3e. As will become apparent later on, it 
will be easier to maintain the original conclusions, viz., stars are not 
produced by neutrons which are captured at the first nuclear collision, 
if we were to assume one-quarter or one-fifth of the light stars to have 
“ recoils ’’ of 4 or 5e. 

In the previous paper a random sample of 411 stars were examined, 
and hence 144-++-25 are estimated to arise from light elements. A further 
202 stars with six or more prongs, 7. e. stars arising from heavy elements, 
were also included in the previous statistics. 279 recoils were observed, 

-of which we assume (144-+25)/3=48+8 are fragments from light nuclei. 


FORMULATION OF THE SCHEME FOR STAR PRODUCTION. 

We consider a star-producing particle (“primary ”’ for shortness) of 
initial energy and momentum E,, P, respectively travelling at an angle y 
to the downward vertical. It emerges from the nucleus with energy and 
momentum E,, P, in the laboratory system, and is accompanied by a 

~*small number of ‘“ knock-on ”’ (grey) protons and neutrons which carry 
away energy and momentum E, and P,. In escaping from the nucleus, 
the knock-on particles will undergo further collisions which leave the 
nucleus, now of mass m, in an excited state and moving with a small 
velocity. Let V, be the downward component of this velocity. We can 
now apply the laws of conservation of energy and momentum : 


E,—E, =E,+ & (=H, say), 


net, 2 ae) (1) 
(P7—P.) n=—P, nm 4 (=P, say). | 


where & is the excitation energy of the nucleus—its kinetic energy is 
negligible—and 7 is unit vector in the downward direction. 

While moving, this excited nucleus will evaporate «-particles, protons 
and neutrons. We assume that this process is isotropic in the centre of 
mass system of the excited nucleus. In general this excited nucleus does 
not possess sufficient energy to break up completely, and a residual nucleus, 
or recoil fragment, will be left. If this possesses sufficient energy it will 
be recorded as a short track. In fact we define a recoil track as one of 
less than 10, and it is now necessary to estimate the distribution in 
mass of recoil fragments by considering the numbers of “ knock-on ” 
and evaporated particles. 

From the 613 stars analysed, 469 are estimated to arise through the 
disintegration of silver and bromine. We assume that the size distribution 
of the 144 light stars is similar to that of silver stars of two to five tracks, 
This is suggested by the size distribution of the gelatin stars. We have 
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assumed silver and bromine stars to have similar size distribution. (The 
final result is very insensitive to the exact assumptions made about these 
size distributions. ) 

Tt can then be shown that for the heavy stars 1:6 “ grey’ protons 
are knocked on per star. Since a heavy nucleus contains about 14 neutrons 
per proton we assume that 1} 1-6 neutrons are also directly knocked on 
by the primary particle.. Therefore 3} fast nucleons are knocked on per 
star. (From the recent Berkeley results it would seem that the n—n 
scattering cross-section is smaller than the »—p cross-section, in which 
case we might expect fewer knock-on neutrons. As will become apparent 
later on, such an assumption would make it easier for us to maintain 
our original conclusions.) 

Since we observe a f-ray in less than 5 per cent of the stars, we conclude 
that the recoils are generally stable fragments. As the average «/p ratio 
of ejected particles is ~0-5 while the original neutron/proton ratio in a 
heavy nucleus is 1-25, it follows that 1-5 neutrons must be evaporated 
for each proton. 

If we now assume that the ratio of silver to bromine stars depends on 
the ratio of their geometric cross-sections and their relative abundance in 
the emulsion, we can now estimate the distribution in mass of the heavy 
recoil fragments. Two peaks are observed in this distribution since we 
start with two groups of nuclei. From the point of view of further — 
calculations it will be sufficiently accurate to consider two “ average ” 
recoil groups, one of mass 93, charge 4le, and the other of mass 67, 
charge 30e. These average recoils occur with relative frequencies 216 
and 253 respectively. Range-velocity curves for such fragments were 
obtained by the method discussed in the previous paper. 


ESTIMATION OF THE ENERGY TRANSFER, EK. 

We have on average 34 nucleons knocked-on by the primary. Therefore 
from the experimental results our best estimate for E, is 90 x 34=315 MeV. 
It has here been necessary to assume that the angular and energy dis- 
tribution of particles producing “ grey ’ tracks is the same in light and 
heavy nuclei, for we use the distribution which was obtained for a random 
sample (~% light stars) to estimate E;, and later P; (see Appendix). 
For stars of medium size (~6 tracks) the distribution of energy of protons 
evaporated from an excited nucleus can be represented by 

P(E) dE= aa exp {(E—V)/T} dE, 
where V=5 MeV. and T=4 MeV. For neutrons we assume the same 
distribution, but with V=0 MeV. A similar distribution can be fitted to 
the «-particle spectrum with V—7 MeV. 

Allowing 4 MeV. as binding energy of an «-particle and 8 MeV. for 
protons and neutrons, we estimate that, for the heavy stars examined, 
the average energy per star taken by «-particles, neutrons and protons 
is 40, 100 and 90 MeV. respectively. Thus &=230 MeV. 

Therefore our best estimate for E is 315-+-230=—545 MeV. 


Origin of Cosmic Ray Stars 67 


ESTIMATION OF THE DownwaRp Momentum TrAnsFsr, P. 


For the three and a half knocked-on nucleons our best estimate of P, is 
3% x 110=385 MeV./c. 

We have now to estimate the downward velocity V, of the excited 
nucleus, and for this we will use the angular distribution of the recoil 
fragments as given in paper I. 

Assume that we can observe all recoils of range greater than 3. This 
assumption is justified below. Our previous results gave 77 recoils 
observed to travel upwards and 133 downwards. We now believe that 
48 of these could be produced from light elements, but unfortunately 
we have no information about their angular distribution. The greater 
the downward excess of recoils the greater the downward velocity of the 
excited nucleus. It is unreasonable to suppose that the light recoils 
are collimated upward. Therefore if we wish to determine the downward 
velocity of the heavy nuclei by subtracting from the observed angular 
distribution of recoils those due to light nuclei, we will obviously obtain 
the greatest downward excess of heavy recoils (and hence greatest 
velocity) by assuming that the light recoils are isotropic, i. e. if 


__ No. of heavy recoils travelling upwards 
_ No. of heavy recoils travelling downwards 


(17477) —(244/24) 
ee (A334 7133) (24-4724) 


R 


—0-49-+40:1. 


(This, then, is the most unfavourable case to take if we are trying to 
confirm our original conclusion.) 

In any given star the velocity of a recoil in the laboratory system will 
be the vector sum of the downward velocity component Vand a velocity Vg, 
where 


hig vassal yl bene. pee ha a f0 (2) 


mp, being the mass of the recoil fragment and P, being the resultant 
momentum of the evaporated particles in the centre of mass system of 
the excited nucleus. In practice there is a distribution of Py which can 
be calculated according to “ Random Walk” theory. Let My be the 
number of stars with N prongs. If .all particles carried off equal 
momentum P’ then the differential distribution p,(P_) dP, is given by 


Chandrasekhar (1943) : 
bs ( ae { peXp (—3PE/2NP) Pp 
pi(Py) al p= KPR{ 2 —PaNP ya VS dP, 


In practice there is also a distribution in the momentum carried off by 
the individual particles, 7. ¢. a distribution p,(P’) dP’ which is calculated 
from the energy spectra of the particles. Therefore we obtain 


aie) ges (— SPZIZNE*) Ms} po(P’) dP’. 


p(Py)=K'PS | = (@nNP?/3)9? 


FZ 
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From this distribution we can therefore obtain the velocity distribution 
p(V.) dV, of the two groups of recoils in the centre of mass system of the 
excited nucleus by using equation (2). 

If a recoil is observed to travel upwards with velocity V it must have 
been ejected upward in the centre of mass system with a velocity greater 
than V. Averaging over all recoils observed with velocity V in the top 


hemisphere we have that the average velocity (Vo)y,=V/fi(Vi/V 
Similarly, in the lower hemisphere, recoils with velocity V in the laboratory 
system must be ejected with velocities V, in the centre of mass system 


smaller than V. On the average. (Vo)gown=Vfo(V1/V). The functions 
f, and f, can be obtained in an analytical form and their numerical values 
obtained for any given value V,/V. It can now be shown that if V4, Vz 
are recoil velocities in the laboratory system corresponding to ranges 
of 3 and 10 p, then 
g Veh 
(1—V,/V)o(V2) dV, 


= FS Vahh 


No. of recoils up _R 


7Vele 


(1+ Vi/V2)p(V2) dV2 


No. of recoils down 


Vt 


We assume a value for V,, calculate the limits to the integrals and hence 
determine R by numerical integration. We can therefore plot a graph 
against V, and hence determine the value of V, for the experimental 
value of R=0-49. The value so obtained is 0-9 x 108 em./sec. 

Averaged over the silver and bromine we obtain mV,—270 MeV./c. 

Hence our best estimate for the downward momentum transfer P (for 
the assumption regarding the light recoils) is 270+385—655 MeV./c. 
(Were we to assume that the light recoils were collimated downwards 
to the same extent as the heavy recoils, this analysis would give 
V,=0-7 x 108 cm./sec., giving an estimate of 595 MeV./c. for P. Taking 
the other limiting hypothesis, in which all the light recoils are collimated 
downwards, our estimate for P would be 425 MeV./c. 


ESTIMATION OF RANGE DISTRIBUTION OF RECOIL FRAGMENTS. 

We have assumed that we can observe all recoil fragments of range 
greater than 3 p.. To test this we calculate the expected range distributions 
of recoil fragments and compare them with the observed distributions in 
the upper and lower hemisphere. ; 

It can be shown (see paper I.) that the differential velocity 
distributions N(V) dV for recoils in the laboratory system are given by 
(a) (Upper hemisphere) 

V+V; 


N(V)=— | AV) ay, 
2ViJyvvuve Va ar 
(b) (Lower hemisphere) 


Vi cvveevat (V7 
N(V)=— = a 2) aV,. 
1 Viens 2 
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These distributions have been worked out (taking V,—0-9 x 108 cm./sec.) 
for the two groups (A—93, Z=41, and A=67, Z=30) of fragments. Using 
theoretical range-velocity curves, calculated according to the method 


given in the previous paper, the range distributions were obtained and 
are shown in fig. 1. 


Discussion oF ERRORS IN P anp E. 


It is here convenient to discuss the errors arising from statistical 
consideration. We assume first that all energy measurements are correct, 
but because of the uncertainty of the angular distribution of the “grey ” 
tracks, the estimated momentum of 110 MeV/c is only correct to 
+20MeV./c (standard deviation—see Appendix). Due to the finite 
number of stars examined, there is an uncertainty in R (for any given 
assumption about light recoils). For R=0-49--0-1 we obtain limits for 
the downward momentum of the excited nucleus of 205 and 330 MeV/c. 
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Range in microns. 
Assuming these uncertainties in momentum can be represented by error 
distributions, it can be shown that the uncertainty in P can also be 
represented by an error distribution with a standard deviation of 75 MeV/c, 
i. e., P=655+75 MeV/c. 

Due to uncertainties in the grain-counting technique, there must be 
an uncertainty in our estimate of E of the order of 10 per cent. However, 
if we overestimate E, we overestimate P and vice versa, and thus the error 
in the ratio P/E is not so important. Therefore, if in the further analysis 
we assume the 10 per cent error in FE, the resultant error m the ratio 
will be slightly overestimated. Thus 
655475 
~ 645-55 


(In these system of units c=unity.) 


P/E = 1:20-0:18. 
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EFrrect of AN ANGULAR DISTRIBUTION OF THE PRIMARIES. 


So far we have estimated the downward momentum transfer. The 
true momentum transfer will be greater. We must now take an angular 
distribution of the primaries, P(y)dy, in order to estimate the true average 

‘magnitude of PISPA 

Assume P(y)=cos? y sin y, where y is the angle with the downward — 
vertical. This distribution has been found for the hard component of 
cosmic rays by counter technique, and seems to hold for the angular 
distribution of the ionizing primaries of stars formed in electron sensitive 
emulsions (Coates 1949, Page 1950). 

We will now investigate the hypothesis that the star primary is captured 
by the nucleus. In this case the resultant momentum transfer will be 
along the direction of the primary particle. Assuming that the energy 
distribution of the primaries is independent of their direction, it may be 
shown that for the above angular distribution, the downward component 
of the average momentum is three-quarters of the magnitude of the true 
momentum transfer. 

Hence, if the primaries are captured, the ratio of the momentum to 
energy transfer is 4/3 x (1-20-+0-18)=1-6-+_0-2,, and it is on the basis 
of this ratio that we test assumptions about the primary particles. 


DIscUSSION OF RESULTS. 


80 per cent of the stars at 3600m. are formed by non-ionizing 
radiation (Brown et al. 1949). We will therefore assume neutrons of 
average energy 545 MeV. produce one star. In this case the momentum 
transfer will be 1150 MeV./c, and thus P/E=2-1. The probability of 
experimental results being consistent with this value is 0-015, and we 
therefore see that, even with the most unfavourable assumptions about 
the angular distribution of recoil fragments from light elements, the 
hypothesis that neutrons are captured in a single encounter with a heavy 
nucleus is untenable. Taking P/E=1-6 for a neutron defines its energy 
as 1200 MeV. The results are therefore consistent with neutrons producing 
a small number, ~1200/545= 2, of stars. 

In 9 per cent of the stars we observe one fast renee (if a proton its 
energy is above 350 MeV.) to be emitted in a downward direction. It 
is reasonable to suppose that this is the primary particle escaping from 
the nucleus, possibly after change of charge. A further 5 per cent of 
stars show showers of fast particles, the majority of which are known to 
be mesons (Camerini et al. 1950). From the point of view of the above 
momentum and energy considerations the ejection of fast mesons 
could take place as an alternative to the re-emission of the primary 
particle. Thus in 14 per cent of stars fast particles are observed to leave, 
and hence it is necessary to assume that in about 50—(9+5)—36 per cent 
of stars fast neutrons capable of further star production must be emitted. 
Hence we should expect the absorption cross-section for the primaries 
to be less than the collision cross-section. Taking the latter quantity 
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Fig. 2. 


Number of tracks per 10° interval in 6—— 


0 20 3840 60 160 


80 100 120 
Down +-6—> Up 


140 180 


Distribution of the values of 0, zenith angle for “orey”’ tracks assumed. to be 
produced by protons of energy 

(2) 30— 50 MeV. 

(b) 50-100 MeV. 

(c) 100-200 MeV. 

(d) 200-350 MeV. 
The shaded portions represent tracks which are likely to have been produced 
by the star primary rather than by a particle ejected in the disintegration. 
The dotted lines represent corrections which have been, applied for loss of 


steeply dipping tracks. 
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as determined by the geometrical nuclear cross-section we do in fact 
obtain a smaller value for the absorption cross-section for lead, aluminium, 
carbon, air and water (see Bernadini, 1950). 
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APPENDIX. 


A. The Angular Distribution of the “‘ Grey” Tracks. 

A random sample of 295 stars formed in electron sensitive emulsion 
exposed at 3600 m. have been analysed with the co-operation of 
Messrs. P. E. Hodgson and T. T. Li. The plates were calibrated by grain- 
counting on p-meson tracks and then the energies of the particles, assumed 


Fig. 3. 


No. of “ grey ” tracks per star. 
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Star size (evaporation tracks). 


Average number of ** grey” tracks as a function of the number of evaporation 
tracks. For comparison, x indicates results of Brown et al. 1949. 


to be protons, which produce the “ grey ” tracks were obtained from grain- 
density measurements. In figs. 2 a, b, c, d we show the distribution in 
zenith angle of the tracks produced by protons 30-50, 50-100, 100-200, 
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200-350 MeV. The variation in the number of such “grey” tracks with 


star size is shown in fig. 3. Good agreement is obtained with the results 
of Brown e# al. (1949). 


The statistics for the angular distributions are somewhat poor, and. 
thus, if we take the energy measurements as correct, we estimate 


(a) Average energy (K.E.+binding energy) per ‘grey’ proton= 
90 MeV. 


(6) Average downward momentum oper “grey” proton= 
110+20 MeV./c. 


REFERENCES. 


BERNADINI, G., CoRTINI, C., and Manrrepint, A., Phys. Rev., 76, 1792. 

Brown, R. H., e¢ al., 19—, Phil. Mag., 40, 862. 

CaMERINI, U., e¢ al., 1950, Phil. Mag., 44, 413. 

CHANDRASEKHAR, S., 1943, Rev. Mod. Phys., 15, 1. 

Coatss, A. C., and Herz, R. H., 1949, Phil. Mag., 40, 1088. 

Harpine, J. B., 1949 a, Phil. Mag., 40, 530; 1949b, Nature, Lond., 163, 
440; 1950, Nature, Lond., 165, 273. 

Pag, N., 1950, Proc. Phys. Soc. A, 63,, 250. 

Perxins, D. H., 1950, Phil. Mag., 41, 138. 


VII. On the Velocity of Second Sound in Liquid Helium I 


By H. N. V. TEMPERLEY, 
King’s College Cambridge *. 


[Received October 25, 1950.] 


ABSTRACT. 


Some recent experiments by Atkins and Osborne (1950) in the Mond 
Laboratory give values for the velocity of second sound at very low 
temperatures in close agreement with Landau’s (1941) predictions. It is 
shown that such results cannot be obtained from any theory of Tisza— 
London type, unless it is postulated that other states besides the very 
lowest one contribute to p,, the density of superfluid. Landau’s (1941) 
method of calculating p, can be applied to any model and does imply such 
a postulate. 

It is concluded that the experimental results only throw light on the 
correct method of calculating p, and cannot be regarded as evidence in 
favour of Landau’s quantum hydrodynamics. On the other hand, the 
fact that liquid He® does not show superfluid properties can no longer be 
regarded as evidence against Landau’s theory, as the predictions of the 
Landau and Tisza—London types of theory are identical for both the perfect 
Bose-Einstein and the perfect Fermi—Dirac gases. 


§1. INTRODUCTION. 


RECENT measurements, Pellam and Scott (1950), have confirmed that the 
velocity of ‘‘ second sound ”’ passes through a minimum value at about, 
1° K. and begins to rise rapidly at lower temperatures. Recent work 
by Atkins and Osborne in the Mond Laboratory (1950) confirms the 
existence of this rapid rise, and shows further that, at very low temper- 
atures, the velocity approaches a value of the order of 150 metres/sec. 
This is in agreement with Landau’s (1941) prediction that the ratio of the 
velocities of propagation of first and second sound should at very low 
temperatures tend to the value ./3. In this paper, we shall show that 
such a result cannot be obtained from any theory which treats liquid 
helium IIT as a condensed Bose-Einstein assembly, the lowest state 
constituting the superfluid, but that this must not be taken to mean that 
the hypothesis of the essential connection between Bose—Einstein conden- 
sation and the properties of liquid helium IJ will have to be abandoned ; 
nor can it be regarded as evidence for the soundness or otherwise of 
Landau’s quantum hydrodynamics (the theoretical validity of which has 
been questioned by some workers). We shall also show that the difference 
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between the two types of theory is not nearly as great as is usually thought, 
because both theories give precisely the same result when applied to the 
perfect gas model, and the absence of superfluidity in pure liquid He? can 
be understood on the basis of either theory. 


§ 2. THe Two-FLuID CONCEPTION. 
The velocity of “second sound” can be shown; see, for example, 
Landau (1941) to be given by the expression 


c= ped” 
Pr(PS/oT) ’ 

where S is the entropy per unit mass. 

This result follows according to either theory and can be derived 
_ independently of any physical model, making the assumptions that 
helium IT consists of a solution or intimate mixture of two fluids, each 
having its own density and velocity field, the normal fluid with density 
p, resembling an ordinary liquid, the superfluid with density p, being 
unable to transport any entropy. This set of assumptions has been found 
to be capable of accounting for a large number of the experimental facts, 
and it is usually referred to collectively as the two-fluid theory. Now 
@S/éT can be found directly from specific heat measurements, and 8 can 
be obtained from them by numerical integration following a reasonable 
extrapolation to absolute zero, so that a measurement of c. may also be 
regarded as a measurement of p,/p,. A quite independent measurement 
of this ratio can be made, Andronikashvili (1948), by measuring the 
angular momentum of a pile of disks on a common axis, oscillating in 
liquid helium II, and identifying the relative density of the superfluid with 
that fraction of the liquid not “ carried round ”’ by the disks. Agreement 
between these two methods is satisfactory down to 1-3° K. 


§ 3. THe DIrFERENCE BETWEEN THE Two THEORIES. 


There seem to be two main differences between the two theories. The 
first difference, which is not a fundamental one, consists of the use by 
Landau (1941) of the Debye model, in which attention is fixed, not on the 
atoms themselves, but on the quantized sound waves that represent their 
relative motion. Tisza (1950), on the other hand, assumes that a satis- 
factory model of liquid helium can be obtained by modifying the energy 
spectrum appropriate to a gas of non-interacting helium atoms in a vessel, 
Evidently the true state of affairs in a liquid involves some pomptenue 
between these two models, one of which we may describe as * solid-like 
while the other is “gas-like”. They appear to have quite different 
statistical mechanical consequences, because the total number of material 
particles is fixed, while the total number of sound quanta may bah 
this difference is only apparent. Either type of model may beapenite 
by the density of energy-levels as a function of energy, or, which ae 
to the same thing, by the relationship between the energy and tota 
momentum of an elementary excitation. Given a particular model of one 
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kind, it seems to be almost always possible to devise a distribution of 
energy-levels so that a model of the other kind may give the same partition 
function, and the two types of model may, in a sense, be regarded as 
two different approximate descriptions of an assembly, rather than as 
embodying different physical assumptions. 

The really important difference arises from the methods of calculating 
p,/p. The most natural assumption to make in a theory of Tisza—London 
type is simply to put this ratio equal to the fraction of atoms occupying 
the very lowest state of the assembly, so that this method is obviously 
not directly, applicable to a model of Debye type. Landau’s (1941) 
method, which consists of imagining the assembly to be in statistical 
equilibrium with a body moving with a very small velocity and then 
calculating its mean momentum statistically, is applicable to a model of 
either type (irrespective of whether his quantum hydrodynamics is sound 
or not) so that a direct comparison of the two methods is possible. They 
give the same result for the perfect Bose—Hinstein gas, and for the perfect 
Fermi—Dirac gas both methods predict that p, should be zero, but, in 
general, the two methods give quite different results, and the experimental 
evidence thus indicates that Landau’s method of calculating p,/p is more 
correct than that suggested above for the Tisza—~London type of theory, 
which is tantamount to saying, that if liquid helium II can be represented 
as a condensed Bose-Einstein assembly other states of the material 
particles besides the very lowest one contribute to p,. Landau (1941) 
also considers, but does not actually use, an alternative method of calcu- 
lating p,/p by imagining the assembly set into rotation with a small angular 
velocity, calculating the mean angular momentum, and comparing it with 
the “rigid body”’ value. He assumes tacitly that the two methods give 
the same result. By considering the perfect gas model by both methods, 
we shall show that they are likely to give substantially the same results, 
though the “ rotation’ method is, logically, slightly preferable. There is 
a possibility that, at extremely low eee ND the two methods might 
not agree for all models. 


§4, THE CALCULATION OF ¢, ACCORDING TO THE LONDON-—TISZA 
TYPE oF THEORY. 


We shall consider only the behaviour of p,/p,, and c, at very low temper- 
atures (as it is already clear that there is no significant difference between 
the predictions of the two types of theory above 1° K.) for which purpose 
we need only specify the low-energy end of the spectrum relating energy 
and momentum of excitations. We assume for this region the following 
relation between the energy E and the length of the momentum vector p. 
K=Ap" where A and n are constants. Now an integration over phase 
space involves an integration over the three components of momentum, 
and the volume element in momentum space is p? dp sin 0 d0 dd. If we 
change the variable from p to E, we see that the model is equivalent to 


assuming a density of energy-levels 47A’E! where t= taal We make 
n 


Velocity of Second Sound in Liquid Heliwm IT We 


the calculation, first by calculating the occupation number of the lowest 
state, then by Landau’s method. For the assembly to be capable of a 
Bose-Einstein condensation, we must have t>+0, i.e. n<3. In practice, 
since the energy of a particle obeying Schrédinger’s equation will always 
contain a term p?/2m, we probably have n<2. 

For the total occupation of all the states but the lowest 


N—N,=47A’ { dK K4/[exp (E/kLT—1)] per unit volume. . (1) 
0 


For sufficiently low temperatures, we may neglect the difference between 
Pils and Pri P> so that 
v3 Am A’ 
ee rei elys, . . 
Ps 

by expanding the denominator of the integrand in powers of exp (— H/kT) 
and integrating term by term. ¢ is the Riemann zeta-function. The 
free energy per unit volume is given by 


Bo4rA'kT | dE Blog [l—exp (—E/kT)], . . . (3) 
0 
from which we obtain for the entropy the expression 

_—— SAAT RTI POEL G2). 5 8 2: (4) 


( 
Substituting for S and dS/OT in the expression for c, we get 
NET (t-+2)L(t+2) 


ee ae Dames ee ok Se ae YR Sala ee (by 
pep (FINE) 

It is interesting that this expression is independent of A’, and only varies 
quite slowly with torn. Inserting for N/p the value 1-6 x 1078 appropriate 
to liquid helium, we deduce that at very low temperatures 6% 45T? 
metres/sec. Thus, if the above method of calculating p,/p,, 18 correct, we 
should expect the behaviour of c, at very low temperatures to be practically 
independent of the distribution of energy-levels. Now it is known, 
Pellam and Scott (1950), that the curve for c, begins to turn upwards at a 
temperature just above 1° K. and they report a value of 33:9 m./sec. at 
some unknown temperature below this. This is not a high enough value 
to exclude the possibility that the law just mentioned is obeyed at very low 
temperatures but very recent work in the Mond Laboratory by Atkins ac 
Osborne (1950), in which velocities up to 150 m./sec. have been observed, 
seems to rule out this possibility quite definitely. 


§5. CALCULATIONS OF C, ACCORDING TO LANDAU. 


We now give the calculation for the same set of models according is 
Landau’s method (1941). “We shall apply it first of all to - ates 
assembly obeying Bose—Hinstein statistics, the Debye model an 3 t . oF 

gas then appearing as particular cases. (There is, hee * 2 <a 
- in applying it to assemblies obeying Fermi-Dirac statistics tor wie ; 
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complete energy-spectrum is known, knowledge of the low energy end of 
the spectrum being no longer sufficient.) The method is as follows :—- 
We suppose the whole assembly in statistical equilibrium with a body 
moving with a very small velocity, v, along the z-axis and then calculate 
the mean value of the momentum, the energy-levels of the assembly first 
being modified to the values that would be measured by an observer 
moving with the velocity v. Now there will be no change in the potential- 
energy if the whole assembly is referred to moving axes, while the change 
in kinetic energy associated with the rth state will be vp”, where p®” is the 
z component of momentum referred to fixed ones. 

The justification of this process is as follows: Schrédinger’s equation 
for non-interacting particles in a potential field V is 


Taking axes moving relative to the original ones with a velocity v in the 
z-direction, we have 
z=2'+00'," tt, 


and the equation becomes 
. : : 
Fig Ee vy, 
877m Qa Ot © 2 Oz 
which shows, from the definition of the momentum operator, that the 
energy is changed by vp, by the change of axes. 

A slight logical difficulty arises for the “ particles in a box’ model, 
because the Schrédinger equation cannot be solved in a compact form 
using the set of axes moving relative to the end wall of the box. It is 
necessary therefore, either to assume that the assembly is “ infinite ”’ in 
the z direction, or else that the z-ends of the “ box ” are joined up to make 
a torus, in which case the model becomes a limiting case of the “‘ rotation ”’ 
model considered below. Ina true ** box ’’ model, the introduction of the 
term vp, can have no effect on the distribution of momentum in the positive 
and negative z directions (in spite of the difference in energies) because the 
amplitudes of the forward and backward travelling waves are constrained 
by walls perpendicular to the z direction to be the same. With an 
“infinite ’’ or toroidal “‘ box ” this difficulty does not arise. Landau is 
therefore calculating the modifications in the occupation numbers of the 
wave-functions referred to the fixed axes that has to occur in order that 
the assembly may be in equilibrium relative to an observer moving with 
the new axes. The method can quite easily be generalized to a model 
requiring the many-body wave-equation. 


> 


Thus, for the assembly in equilibrium with a moving body, we must 


replace K,, the energy of the rth level, by E,—.vp”. The mean value of 
the momentum in the z-direction is then ; 


ede 
9 J oJ yt A. exp. (Heme) al ceal an wena ityc: allies 
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which; to the first order in v, may be written 


2% rm fo a 1 
ize ‘ / 0 ad ‘ see | ee ee eae . 
lk ft fi p COS & . pu cos 0 an (; ah EEN) p> dp sin 6 dé dd. 


(7) 


Since E= Ap” we have 0/0K = (1/nAp"-1)(@/dp) and we may now integrate 
by parts with respect to p, the integrated term vanishing at both limits 
(provided <5). The integration over the angular coordinates in momen- 
tum space leads to the factor 47/3 so that the total momentum of the 
assembly is then 


47 Bn pt” dp ae 

rs aban i, Nexp (E(p)/ET)—1 per unit volume, area (Sy 
and the effective mass, identified by Landau with the ‘‘ mass of the normal 
fluid ” is obtained simply by dividing this by v. For an assembly such as 
the Debye model, where the number of quanta present is not fixed, we 
have A=1, while for a degenerate Bose—Kinstein assembly d is very nearly 
unity and we have 


hacks _(e De ID de | 
Nat"), pep + 


the energy of the lowest level being taken as zero. For the perfect gas 
E=p?/2m so that n=2 and A=1/2m, in which case expression (8) is 
exactly equal to the integral in equation (9) multiplied by mv, which means 
that the atoms in the very lowest state do not participate in the motion, 
this state contributing to equation (9) but not contributing to expression 
(8) because the momentum of the very lowest state, referred to the fixed 
axes, is zero for a large assembly. Thus, for the perfect gas, both methods 
of calculation give the same result for p,/p and therefore for c, also 
whereas, if the energy is any other function of the momentum, as will occur 
if there are interactions between the atoms, the results of the two methods 
of calculation may be quite different. For example, if n equals unity, we 
conclude that c, tends to the constant value A/\/3 at low temperatures, 
just as was found by Landau for phonons for which A=, the velocity of 
ordinary sound, whereas equation (5), with t=}, gives the result to be 
expected using the London—Tisza type of approach, plus the assumption 
that the superfluid is identical with the lowest state of the assembly. 

_ Equations (6)-(9) can easily be modified to suit vie Fermi- Dirac 
assembly and, for the perfect gas, we conclude that the “ effective ~ mass 
of the assembly is precisely the same as its actual mass, so that nothing of 
the nature of a “superfluid component ” is to be expected according to 
Landau’s method of calculation. The non-appearance of superfluidity in 
liquid He? cannot therefore be regarded as evidence against Landau’s 

theory. 
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§6. Are Lanpau’s Two METHODS OF CALCULATING p/,p,, EQUIVALENT ? 


We return to the perfect Bose-Einstein gas, but we suppose it to be 
contained in a cylinder of radius a and length d. The possible wave- 
functions and energy-levels are 


laz 


; ar : 
bi, ny =Sin mF J, (2:7) exp (v6), 


E,,,—= [Pe ra ee 
ire ca > 


82r2m | d2 a 


where j,,, is the nth zero of J,(x) which is given by the asymptotic 
expression j,, ,~(n-+43v—})7 which is, in fact, fairly accurate even for 
small values of m andv. By forming the partition function, replacing the 
sums by integrals in the usual way, and changing the variables of 
integration from n and v to n+4v and n—}y, we can easily show that, 
as we should expect, the partition function is in the limit, the same for 
a cylindrical enclosure as for a rectangular one and that the assembly 
shows the phenomenon of “ condensation’ into the lowest level in 
exactly the same way. It is possible to deal with a toroidal box by 
introducing Y, as well as J,. However, if we imagine the whole assembly 
in equilibrium with a body rotating with a small angular velocity and 
make a calculation analogous to that of equations (6) to (9) we reach a 
similar conclusion, with the difference that all the states with v equal to 
zero are devoid of angular momentum, and can therefore contribute 
nothing to the mean angular momentum of the assembly. The “ number 
of superfluid atoms ” therefore now comprises not only the number in 
the very lowest state, but all those with v0. The number N, of atoms. 
in such states, apart from the very lowest, can be obtained by summing 
the ordinary expression for the occupation number over all possible 
values ot / and n, which gives, if the sums are replaced by integrals : 


ee [" fz dl dn 
OW J, Jy Nexp [OP/BmET) (PP pret 


This integral diverges if A is precisely equal to unity, but, since Ny=1/(A—1) 
expression (11) implies 
T\* F ones 
N,=N,+Ni=N | — -+- T O(N log No); . . (12) 
0. 0 

where T, is the transition temperature. Thus the difference between 
the two methods is actually negligible for the perfect gas model for a very 
large assembly unless the temperature is extremely low, but this 


calculation does at least show the possibility of the ‘‘ linear’? and 
“ rotational ” definitions of p,jp,, giving differing results for some models, 
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§7. CONCLUSIONS. 


(a) The assumption that Helium II is a condensed Bose—Kinstein 
assembly and that p,/p is proportional to the occupation of the lowest 
state is disproved by extending measurements of C. to lower temperatures. 

(6) This experimental decision does not necessarily mean that liquid 
helium does not show Bose-Einstein condensation or something similar. 
It is possible to retain this assumption, provided that one admits that 
other states besides the very lowest one contribute to p, in all cases 
except the perfect gas, which seems to be an “ accidental ” exception. 

(c) The fact that He? does not show superfluidity does not, by itself, 
mean that Landau’s ideas are wrong. As we have seen, Landau’s (1941) 
method predicts no superfluidity for the perfect Fermi gas. No 
reasonable analogue of the Debye model is known that will represent 
the excitations of a degenerate Fermi assembly, so a calculation using 
such a model cannot be made. 
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SUMMARY. 


Cosmic-ray stars in lead have been investigated by means of photo- 
graphic plates in close contact with lead foil. The energy distributions of 
protons and alpha-particles from lead stars have been determined and 
are compared with the prediction of the evaporation theory. The mean 


mate value of the rate of star production in lead deduced. 


ratio of doubly to singly charged particles is estimated, and an approxi- 


§1. INTRODUCTION. 

Tue photographic plate method has already been used (Harding, 
Lattimore and Perkins 1949) to investigate nuclear disintegrations in the 
emulsion caused by high energy cosmic-ray particles. Unfortunately, 
however, the disintegrations cannot be ascribed to any particular element, 
apart from specific nuclear reactions in the light elements of the emulsion, 
which can be identified in exceptional cases. Also, although stars from 
which more than eight charged particles are emitted must be due to the 
disintegrations of one of the heavy elements of the emulsion, there still 
remains the ambiguity between silver and bromine. 

In order to overcome this difficulty, Harding (1949) has used sandwich 
emulsions composed of alternate layers of emulsion and pure gelatine. 
Any star observed to start in the gelatine layer must be due to the 
disintegration of one of the constituents of gelatine, viz. carbon, oxygen 
or nitrogen. Attempts have also been made in this laboratory to study 
disintegrations of lead nuclei by incorporating layers of lead phosphate 
in the emulsion. However, the amount of lead which can be incorporated 
in this way without impairing the transparency of the emulsion is so small 
that to obtain a reasonable number of stars a prohibitively large plate 
area would have to be scanned. 

In the experiment described in this paper we have exposed plates with 
the emulsion in contact with lead foil, and analysed the stars originating 
in it. 

§ 2. EXPERIMENTAL PROCEDURE. 

Lead foil was rolled on to the emulsion of Ilford C2 100 micron nuclear 
research plates. It adhered well but, in order to ensure that there was 
no gap between the foil and the emulsion, the plates were tightly clamped 
together. 


* Communicated by Sir George Thomson, F.R.S, 
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After exposure to cosmic rays on the Jungfraujoch, the lead was peeled 
off the emulsion, which was then developed in the usual manner. The 
plates were scanned twice, at magnifications of 130 and 270, to reduce 
the possibility of missing stars originating in the lead. 

C2 emulsion was used instead of the more sensitive G5, firstly because 
the large number of tracks crossing the latter emulsion makes searching 
difficult and also increases the number of apparent two-track stars in the 
lead due to the chance concurrence of unrelated tracks ; secondly, C2 
emulsion gives good discrimination between proton and alpha-particle 
tracks, permitting the («/p) ratio to be found: thirdly, since only 
evaporated particles are recorded in C2 emulsion, it may be assumed 
that the particles observed were isotropically emitted from the lead nuclei, 
an assumption which greatly simplifies the calculation of their probabilities 
of observation. 

Each lead star was examined and the numbers of protons and alpha- 
particles from it which ended in or passed through the emulsion were 
recorded. The distance of the centre of the star from the emulsion surface 
was found by a simple geometrical method. The projected range of each 
track and the difference in depth between the two ends of its emulsion 
path was measured and hence, assuming the thickness of the emulsion 
to be exactly 100, the range in emulsion of each particle was found. 
A geometrical calculation gave the range of the particle in the lead before 
it reached the emulsion. Any deviation of the emulsion thickness from 
the assumed value of 100 produces a proportional error in all measure- 
ments perpendicular to the plane of the emulsion. Measurements on 
the undeveloped emulsion of similar plates led to the conclusion that such 
deviation can be of the order of 5 per cent. 

Grain counts were made on all tracks of protons which passed through 
the emulsion except those making an angle of more than 56° with the plane 
of the emulsion. The latter were rejected, as grain counts on such steep 
tracks are unreliable. The energy—grain-density curve for the emulsion 
was constructed in the usual manner, using the tracks of protons and 
mesons which ended in the emulsion. 


§ 3. SumMMARY OF RESULTS AND METHOD OF ANALYSIS. 

The range and grain density measurements were used to determine the 
energy distributions of the protons and alpha-particles in the way described 
in §4. The results are shown in figs. | and 2. The mean energy of 
emission of the protons was 14-2 MeV., and that of the alpha-particles 
27-7 MeV. 

The apparent proton size distributions of the stars observed in the 
first and second fifty micron layers of lead from the emulsion surface are 
shown in Table I. The term ‘proton size distribution ~ means the 
distribution of star sizes obtained by counting only the protons associated 
with the stars. The constants of the best exponential curves to fit 
these distributions were found by the method of least squares. The 


G2 
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justification for assuming that the size distribution is of exponential 
form is discussed in §7, Then, using the size distribution conversion 


Fig. 1. 
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Corrected alpha-particle energy distribution for stars in lead based on 28 tracks, 
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functions derived in § 8, the proton size distribution of the stars in lead 
was deduced. Hence, using the value of the («/p) ratio obtained in § 6, 
the total size distribution of the lead stars was found by an analogous 
use of the conversion functions. From the time of exposure and the 
area of plate scanned the rate of production of lead stars at the J ungfrau- 
joch was deduced. These calculations are summarized in Table IT. 

A lead star, from which four protons and two alpha-partarcles reach 
the emulsion, is shown in Pl, V. 


§ 4. DETERMINATION OF THE PROTON AND ALPHA-PARTICLE 
ENERGY DIsTRIBUTIONS. 


All the tracks due to particles from stars in the lead from which at 
least three particles reached the emulsion were examined. “Stars” 
having only two apparent tracks were excluded since many of these are 
due to the chance concurrence of unrelated single tracks (see § 10). 

The energies of emission of particles which ended in the emulsion were 
evaluated directly from their ranges in emulsion and lead using the 
respective range-energy curves (H. Bradner et al. 1950, Montgomery 
1949). 

‘The energies at the mid-points of their emulsion paths of those particles 
passing through the emulsion were evaluated from the grain densities 
of their tracks. Then, as before, the energies of emission of the particles 
were calculated. 

Kach track was weighted by the reciprocal of the probability that a 
particle of that energy, emitted from a nucleus at that distance from the 
emulsion surface, should enter the emulsion or pass through it, whichever 
is applicable. Allowance was made for the tracks rejected on account 
of their steepness. 

The proton and alpha-particle energy distribution histograms obtained 
are shown in figs. 1 and 2 together with evaporation curves calculated 
as follows : 

Assume that 

(a) The nucleus evaporates down the Heisenberg valley with 
ON/OZ=2°3. ; 7 

(b) The binding energies of the emitted particles are Q,=Q,=7'7 MeV., 
Q,.=2°6 MeV. 

(c) Only protons, alpha-particles and neutrons are emitted in the 
evaporation process. 

In addition we have, from the experimental results : 

(d) The energy distributions correspond to stars from anion on ne 
average about eight protons are emitted. This value was calculate 
from the true and observed proton size distributions and the distribution 
in depth of the stars in the lead. Abe Msc 

(e) From the energy distributions the mean energies of emission are 


E,,=14-2 MeV., E,=27-7 MeV. 
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(f) The («/p) ratio for lead stars is 0-34-+-0-07 (§ 6). Thus about three 
alpha particles are associated with a star from which eight protons 
are emitted. Also, to satisfy assumption (a), twenty-six neutrons must 
be emitted. 

Let the mean nuclear temperature be T. Then 

Energy removed by protons is (14-2--7-7) x 8 MeV. 

Energy removed by alpha particles is (27-7-+2-6) x 3 MeV. 

Energy removed by neutrons is (2T+7-7) x 26 MeV., 


7. €. mean excitation energy U=52T+466 MeV. 
It has been shown by Le Couteur (1950) that 


hs 
3 Riga lee ai /U. 
: ’ A ; ; 
Also, for silver plus bromine, A= war where A=atomic weight. 


—- 2 A2 
Ue a =37-5T?, whence T~4°3 MeV. 


From the evaporation formula (Harding, Lattimore and Perkins 1949) 


E—V E—V 
P(E) dE= ps °XP (- =) dk, 


where V is the height of the nuclear potential barrier, it can be shown 
that the mean energy of emission E=2T+V. 


V,»=E,—2T=14-2—8-6~6 MeV. 
and V,=E,—2T =27-7—8-6~19 MeV. 


where the suffixes p and « refer to the values of V, T and E appropriate 
to the proton and alpha-particle energy distributions respectively. 

The evaporation curve with T=—4, V=6 fits the experimental proton 
results quite well, which shows that the energy-temperature relationship 
in lead is consistent with that in silver plus bromine, within the accuracy 
of the experimental results. 

The value of the potential barrier height for protons is considerably 
less than the 9-6 MeV. obtained by interpolation from the results given by 
Bethe (1937), which, however, were determined from low energy scattering. 
This difference may be ascribed to (i) violent oscillations of the nucleus 
due to the high excitation causing local lowering of the potential barrier, 
(i) quantum mechanical penetration of the barrier, and (iii) broadening 
of the energy distribution by experimental error. For emulsion stars in 
silver plus bromine of a similar size to the lead stars examined, Harding, 
Lattimore and Perkins obtained a value V~4 MeV., which is less than 
the value for lead, as expected. The number of nucleonic collisions made 
by the primary in the target nucleus increases with A, but the resulting 
temperature for the same initial excitation decreases as A increases. 
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The temperature of 4 MeV. found for the lead disintegrations is less than 
the 53 MeV. found by Harding, Lattimore and Perkins for silver plus 
bromine, which indicates that the latter effect outweighs the former, 
although the experimental errors make comparison uncertain. 
Considering the alpha-particle energy distribution, we should have 


V,=2V,. ny a= Le, MoVe; ~. Ty=4(27-7—12)~8 MeV. 
The experimental results agree better with V,=12, 1l=—8 than with 
V,=19, T,=4. This indicates that the alpha-particles cannot all be 
considered as sharing in the normal evaporation, a conclusion also reached 
by Harding, Lattimore and Perkins from the examination of silver plus 
bromine stars. They suggested that the alpha-particles are evaporated 
off in an initial process of higher excitation and, therefore, higher nuclear 
temperature. Another explanation, which does not exclude this process, 
is that some of the alpha-particles are not evaporated at all, but come 
from the disintegration of highly excited nuclear splinters (Perkins 1950, 
Hodgson 1950). 

In these calculations no account has been taken of the deuterons and 
tritons which, since they form about 15 per cent of the particles from 
silver plus bromine stars, are probably emitted also from lead stars. 
There is little advantage in increasing the accuracy of theory and experi- 
ment until the proportion of such particles has been evaluated. This can 
only be done by a series of grain counts on long tracks parallel to the 
surface of the emulsion. From the very nature of this experiment such 
tracks are not available. 


§ 5. DETERMINATION oF P,, AND P,,. 

Let P,, and P, be respectively the probabilities that a proton or an 
alpha-particle from a star in the lead is observed in the emulsion. 
Both are functions of f, the distance from the star to the emulsion surface. 
Since C2 emulsion, in which only low energy particles are recorded, was 
used, it may be assumed that the particles were isotropically emitted from 
the lead nuclei. 

Particles whose directions of emission make a small angle with the 
~ emulsion surface have less chance of being observed than those making 
a large angle, both because of their limited range and the increasing likeli- 
hood of their being missed due to the limited field of view of the 
microscope objective. Therefore 

Bey mw sin 8 cf, ayD(f, 8) dd, 
Bo ed 
where 6 is the angle the track makes with the perpendicular to the 
emulsion surface, and R(f,4) is the probability that a track has range 


greater than f sec 9, 7. €. 


o EV E-V\ 
R(f, a)= ibe 2 exp & “r) dks, 
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where T and V are the constants of the energy distribution and K, is the 
energy of a particle of range f sec 0 in lead (Montgomery 1949). D(f, 0) is 
the probability that a track of inclination @ passing right through the 
emulsion from a star f microns in the lead is found. Since it depends on 
the distance from where the particle enters the emulsion to the foot of 
the perpendicular from the star to the emulsion, it may be written 


D(f, 0)~D(f tan 6). 
D(f tan @) is unity for values of f tan @ from 0 to about 100 microns, after 
which it drops quite steeply to zero. It depends, of course, on the care 


with which the plates are scanned, and on the field of view of the objective 
used. 


Fig. 3. 
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The distribution in @ of the proton tracks from stars from which at least 
three particles entered the emulsion was measured for several ranges of 
f and compared with the curve of } sin @ R(f, 6).. This curve represents 
the distribution expected if D(/,#)=1, that is, no tracks missed, and 
it was calculated from the appropriate energy distribution (see § 4). 
From the amounts by which the experimental results fell short of those 
expected, D(f tan @) was calculated for several values of ftan@. The 
mean value of D(f tan @) was used to find P,, and P, as functions of f by 
graphical evaluation of the above integral. The results are shown in 
fig. 3. 

These curves show that it is best to ignore all stars for which f> 100 ;1, 
as the probability of particles from them being detected is small. 


§ 6. THE «/p Ratio ror Luap Srars. 
The numbers of alpha-particles and protons from stars in various ranges 
of f were counted and divided by the corresponding P. This gave a 
mean value for the «/p ratio, averaged over all star sizes, of 0-34-L0-07. 
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aia as ae the x/p ratio means, more strictly, the ratio of the 
ce iui Sie charged to the number of singly charged particles 

- # shows the variation of mean «/p ratio with nuclear ch: aed 
on the results of several authors. . sce bana 


Fig. 4. 
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§ 7. Tue Size DistRIBUTION ASSUMPTION. 


Several authors (George and Jason (1949), Page (1950) and unpublished 
work at Imperial College) have shown that the size distribution of stars 
formed by cosmic rays in C2 emulsion at 3500 m. may be represented 


by the formula 
N,»=Np exp (—én), 
where N,, is the number of stars of n tracks, N, is a constant and 
ED'S for 3<n<8 
and ¢-~0-33 for n> 8. 


: 
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These figures for the differential size distribution have been calculated 
from those given by George and Jason for the integral size distribution. 

Emulsion stars can be formed either in the light elements carbon, 
oxygen and nitrogen, or in the heavy elements silver and bromine. The 
largest star that can be produced in the light elements is one of eight 
tracks in oxygen, but, since the stars in light elements have about the same 
numbers of alpha-particles as protons, very few light element stars have 
more than five tracks. 

These observations strongly suggest that the emulsion size distributin 
is composed of the superposition of two exponential curves, one of the 
form exp (—0-33n) for the heavy elements, and another, of greater slope, 
for the light. The slope of the light element curve may be calculated 
from the figures given above, and is found to be ~0-6. 

It therefore seems reasonable to assume that the size distribution for 
lead stars has exponential form also. 

The assumption is only approximate, since it gives a finite number of 
stars with n>Z, whereas the true number is zero. 


§ 8. Size DisTRIBUTION CONVERSION. 


Suppose the stars produced in the lead have a proton size distribution 
given by 


N?=N, exp (— €,”) (n>2), 


where n—=number of protons evaporated from the nucleus, 
N?=number of stars with n evaporated protons, 
N, and €, are constants. 

N, depends on the intensity of the star-producing radiation, the number 
of lead nuclei and the time of exposure, and &, on the material in which 
the stars are produced and_on the energy distribution of the star-producing 
radiation. 

Then an n-track star in a layer distant between f and f+df from the 
emulsion surface will have a probability of being observed as an m-track 
star in the emulsion given by 


n! 


(n—m)imi-? | ekg (l—P, rn, 


where P,, is a function of /. 
The total number of m-track stars from the layer df observed in the 
emulsion will therefore be 


oe n! ; f 
M,, Ole * <n (n—m)!m! Seg =P) eeNG exp (— E,) df. 

N is the number of tracks in the largest possible star in lead. Since the 
number of stars with n>N predicted on an exponential law is extremely 
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small for lead (less than 110-4 per cent), N may, without appreciable 
error, be replaced by infinity in the above sum. 


ee Te! ; 
M,.= era) ml P@(1—P,)"-™N, exp (—E,n). 


=N,P; exp (— pm =(L—— Ps) exp (ein 1) 


ee NS { P,, exp (—é,) ia 
ES ae) exp (=£,) L(l Pe) exp (5) 


=M,; exp (—&m), 


where 


N exp (¢,)—(1—P,) 
ee and exp oe Be 
 1=(1—F,) exp (—é,) P (es) P, 

This is the observed proton size distribution of the stars produced in the 
layer df. Thus, provided P,, is known, the actual proton size distribution - 
of the stars in lead at any particular depth f can be calculated from the 
observed proton size distribution. 


§ 9. THe Tota Size DistRIBUTION. 


The method of determining the proton size distribution has been 
outlined in § 8, and that for the («/p) ratioin § 6. _‘It has also been shown 
(unpublished work at Imperial College) that the distribution of alpha- 
particles and protons among the stars is nearly random in the sense that 
the emission of one particle gives no information on the probability of 
emission of another similar particle. Recent work (Perkins 1950, 
Hodgson 1950) has shown that close pairs or triplets of particles from the 
disintegration of unstable nuclear fragments are sometimes emitted from 
stars, but the frequency of such events is so small that their effect on the 
randomness of emission can be neglected. 

The total size distribution of the lead stars was deduced from the 
proton size distribution in exactly the same way as the latter was found 
from the observed proton size distribution. The formulae of §8 were 
used with é, replaced by &, & by &, No by N, and M, by No. P was 
replaced by the probability that a given track from a star is a proton,— 
that is, p/(«+p)=[(«/p)+1]*. 

The total size distribution of the stars in lead is then given by 


N,=Ny exp (— én), 
where n=number of charged particles in a star, 
N,,=number of stars with n charged particles, 


Ny, € constants as before. 
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§10. EXPERIMENTAL RESULTS FOR THE RATE OF STAR PRODUCTION. 


The proton size distributions of the stars observed in the lead are 
shown in Table I. The numbers of apparent two-track stars have been 


Taste I. 


Numbers of Lead Stars with Different Numbers of Protons, with Two-track 
Stars corrected for Chance Coincidences of Single Tracks. 


Number of protons 


Range of f 


0<f<50 wu 


50<f<100 p 


corrected for those due to the chance coincidence of unrelated single 
tracks. These are quite a serious source of error, being as numerous as 
the real two-track stars on the plates used, and hence the final number 
of two-track stars is subject to considerable uncertainty. The number 
of spurious three- or more-track stars is negligible. 

The best values of € and M, to fit these distributions, together with the 
proton and total size distributions calculated from the observed proton 
size distribution, are shown in Table II. The close agreement between the 


TABLE II. 


Results of Calculation of Total Size Distribution. To convert to stars per 
nucleus per day at the Jungfraujoch, multiply the numbers of stars 
obtained by substituting the above values in the exponential formula 
by 5x 10-*4, 


Apparent Real Total 
proton size proton size size 
Range of f | distribution distribution distribution 


fy | My Ep | No |[(x/p)+1]> 


0<f<50 p 0-45 | 22-8 0-21 | 12-1 0-75 0-17 9-5 


50<f<100 w | 0-71 | 52-5 0:24 | 21-4 0-75 0-18 | 16-8 


final values of € is probably fortuitous in view of the poor statistics. 
However, the slope of the lead size distribution is certainly less than that 
for silver plus bromine, which is expected from the considerations of § 7. 


Disintegration of Lead Nuclei by Cosmic Rays 93 


The variation of slope of size distribution with nuclear charge is shown 
in fig. 5 together with the results of other authors. Those of Barhaur and 
Greene were obtained with G5 plates, and so are rather lower than the 
other values of €, which refer to stars in C2 plates. The result of Hardin | 
obtained by the gelatine sandwich method, is based on a rather ean 
ae of stars, and therefore has less accuracy than the value anime 
in § 7. 

The final values of Ny, which should be equal, are only considered 
reliable to within about a factor of two. | 


Fig. 5. 
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§ 11. CONCLUSIONS. 


The method outlined here for investigating cosmic-ray stars is only 
suitable for materials of high nuclear charge. This is due to the relatively 
small probability (~0-3 on the average for stars with 0<f<100,) that 
a track from a star outside the emulsion will be observed, and to the 
variation of size distribution slope with nuclear charge. This results in 
stars in light elements appearing in most cases as two-track stars, if at all. 
It is very difficult to correct accurately for spurious two-track stars: 
indeed, it is probably best to neglect them altogether. Therefore only a 
very small fraction of the stars from a light element will be analysed, and 
thus a prohibitively large plate area would have to be scanned to find a 
sufficient number. 

For light elements the sandwich method of Harding is therefore superior. 
But for heavy elements, which cannot be incorporated in a sandwich in 
satisfactory quantities, the method described here is preferable, 


94 On the Disintegration of Lead Nuclei by Cosmic Rays 


Recently, Barbour and Greene (1950) have exposed metal foils with 
emulsion on either side. This doubles the value of P, and so the above 
objections do not apply to this method providing the foil is thin and is 
in good contact with both emulsions. These conditions are, however, 
very difficult to fulfil in practice. 

All methods will tend to give a low value for the relative cross-section, 
because stars outside the emulsion are very much more difficult to find 
than those in it. 
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Disintegration of a lead nucleus from which four protons and two alpha particles 
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IX. Landau Diamagnetism and Meissner Effect. 


By A. Papaperrou, 
Physics Department, University of Manchester * 


[Received October 26, 1950.] 


Summary. 
A magnetic field of the form 


Hy: 27H, cos (27%) ,  H,=H,=0, 


is considered, and its influence on free electrons is discussed to the first 
approximation. A field of the form 


H,=27H, . cos (2" ;) . COS (2- *) 


has been considered too. It is shown that, if in a superconductor the 
distribution of the electrons in momentum space differs a little from the 
spherical one, an internal field of the second type will be formed inside 
the superconductor ; the Meissner effect then follows for not too strong 
external fields. 


§1. InTRoDUCTION. 


LANDAU’S first investigation of the influence of a magnetic field on the 
orbital motion of free electrons has been followed by a number of similar 
calculations, for free electrons as' well as for real metal-electrons. All 
these calculations refer to a homogeneous magnetic field, this being 
justified by the very small magnetic susceptibility of the metals under 
ordinary conditions. 

An entirely different situation presents itself in the case of a super- 
conductor, in which an external magnetic field penetrates only to a very 
small depth (order of magnitude 10~° cm.). Instead of considering a 
homogeneous field now it is more reasonable to split the magnetic field 
into its Fourier components. In the simplest case of a superconductor 
in the form of an infinite plane slab (of constant thickness) we get the 
one-dimensional Fourier components 


H,=2rHy cos( 274), H,=Hj=0. “2... (I) 


(Oy-axis normal to the slab, Oz-axis in the direction of the field.) The most 
important components will have wavelengths of the order of magnitude 


* Communicated by the Author, 
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of the penetration depth, i.e. A10-5cem. For a superconductor of 
prismatic form the Fourier components will be two-dimensional. The 
most important components will be 


Bo=27H, cos (275) cos (2=), .| a 


again with A10~° cm. 

In the present paper we shall investigate the behaviour of the metal- 
electrons in a magnetic field firstly of the form (1), and then of the form (2), 
in both cases to the first approximation only. As we shall see later, the 
results of this investigation show us a possible explanation of the Meissner 
effect. 

Strictly speaking we ought to base our calculations on the properties 
of the real metal-electrons, 7. e. to consider not only the magnetic field (1) 
or (2), but also the periodic electrostatic field of the lattice. However, 
such a calculation would be extremely complicated, if possible at all. 
Hence, we shall restrict ourselves to the case of free electrons. The 
results will, in any case, have only qualitative significance. 

In the first section we shall consider the one-dimensional magnetic 
field (1); in the second section the two-dimensional field (2). In the 
third section we shall discuss the results of those calculations in connection 
with the Meissner effect. 


§2. THe ONE-DIMENSIONAL FIELD. 


We assume that the electron-gas occupies the whole space and we 
impose on the wave-functions the usual periodicity condition with a 
period-length a. The unperturbed wave-functions are then 


Hane =@re exp d= (Ee-+my tte), et eee 


€, , ¢ being integer numbers (positive or negative). The unperturbed 
energy is 


Koen (S sh 02)a cs aie aah oes 
with 
h2 
“1 Sma . . . . . . . ° e (4a) 


The vector-potential of the field (1) is 
A,=—H,A sin (20% ALA 220 
x 0 iN > os a . ‘ . . (5) 


In the general formula for the perturbing potential 
e2 


v= 
2mc 


zeh 
a (A grad)-+ A2 
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we shall omit the second term, since we restrict oursleves to the first 
approximation. We find then immediately 


ehHyo€é 
Vibe Rea 2imeny (be, H+, ce N—No; t)s . . . : (6) 
a aay 
phere N= > = positive integer. oe oe aie 8 
The perturbation method gives then for the perturbed wave-function 
we 1eH EA? iF i 
Pen, tHE, 0, or he Caa7 eS In—No 1 lite .) 5 Ae (8) 


‘This formula shows that we shall have single states only if we assume 
ie —oanimieper, ~~. . |. 2. . . (Ta) 


since in this case there will be no vanishing denominators in (8). The 
final results of the calculations are the same for n»=even or odd. We 
ean, therefore, introduce the assumption (7a)—which simplifies the 
summation over the occupied electron states—without imposing any 
unnecessary restriction. The magnetic energy will then be zero for all 
states, 

DUO eames eee wc tos) ees 2s (9) 


7. e. in this approximation the energy of the electron state reduces to (4). 
The state described by the wave function ¥ is associated with a current 
density s given by the general formula 


ee rgd Pade) 
~ 4. 4me st ae le, me 


OU A. 


Introducing in this formula the expression (8) we find 


A ra ie y 2@/ 1 I heé 
26 stan es SN Gh ee = ait | 
Sa MC?ANy sin ( 74) j a N (ae a)  inea > 
The last term, which does not contain Hy, can be disregarded because it 
gives a vanishing contribution when we sum over all occupied states. 


The components s,, 8, contain only the term corresponding to the last 
one in (10), and consequently we may put 


$,=8,=0. 


Taking into account (5) we may write (10) in the form 


eA, 22/1 1 
= — — — : he (10 
fe mera° [1+ 1 ea sea) oat 


Before going further it is worth seeing up to what values of Hy our 
formuls will be valid. The condition for the validity of the perturbation 
method is that the coefficients of #- ,,.,,,¢ in (8) be small compared with 1. 
For the large majority of occupied states we have |€| Emax and 
] 2n7% |~max='max: Therefore, the condition for the validity of 
our formule will be 


eH,A” 
he 
SER. 7, VOL. 42, NO. 324.—JAN. 1951 H 


eles 


100 A. Papapetrou on the 


Now we see that the quantity 1/n(w{(y)) changes slowly even in the 
neighbourhood of y=0, and consequently we can evaluate the sum 


appearing in (15) Py; integration : 
y + 00 dy \ 
Zi wiid=—4e | Peep ee) 


We are interested in the case 
éy> kT, 


in which this integral ean be evaluated by a well-known method*. The 


result is 
kT\?2 
12, 1/2 es 
2 a(n) ae E ae *) | 


Introducing this in (15) we get 


2a e'n2zA, Eo a (kT\2 ‘ 


The dependence of S,, on the temperature is not shown directly in (15 a) 
because ¢, itself is dependent on T. The dependence of «, on T can be 
determined by calculating the number N of electrons in terms of e) and T. 


"The result is 
N= or Es ‘|=2 ca ae 
3 éy ? 


if we denote by éo) the value of «, Be c= eae Eo In (15a) 
and taking into account (4 @) we finally find 


Qn eA, kT 
$,=- = wanen VOM) [1B ( aN Rema) 


£00. 
This result can be brought to a more familiar form. Let us calculate 
a distribution of magnetic moment m equivalent to the distribution of 
current S by means of the general formula 


S=curl m. 


If then we introduce the magnetic susceptibility x, 


m= yH, 
we find (provided that x is independent of x, y, z) 
S=y-.curkH,-). aj ae (17) 
Hence in our case ; 
0H 
x8. 
Considering (16) we find 
_ __ &ao/(2mMe qo) kT 
x= 6amczh [1- 5 (— ) |, thin? arena Ls) 


* See e.g. Fréhlich, Elektronentheorie der Metalle (Springer 1936), p. 357. 
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which is identical with the Landau susceptibility of free electrons in a 
homogeneous field. This is not surprising at all, since e.g. for 
A~1 cm. we should have demanded that our calculations lead to the 
Landau susceptibility (18); our conclusion is now that this is generally 
true, if only 1% <apax*. 

Coming back to equation (10 a) we stress an important feature of it. 
The current s, produced by the magnetic field (1) acting on any given 
electron state is not of the order of S,,|N, but very much larger. The smallness 
of S, (or of x) is the consequence of a very complete cancellation of much 
larger terms, this cancellation being due to the spherically symmetric 
distribution of the electrons in momentum spacet. If we disregard the 
states with ¢°<nj—which are relatively few because of the assumption 
% <Nmax—=€max—, the second term in the bracket of (10a) will be 
predominant. We see then that we shall have a strong paramagnetic t 
current. if || <m)/2, or a strong diamagnetic current if [7 | >n,/2. We 
shall make use of this fact in the last section of the present paper. 

Instead of the metal extending infinitely in all three dimensions we 
may equally well consider a slab of finite thickness in the Oy-direction, 
e.g. the one filling the space 0<y<a. The calculations are very similar 
to those given above. But there is one characteristic difference in so 
far that we cannot now start with the field H, given by (1). The reason 
is that, as can be seen immediately, the quantity Vy cannot in this case 
be developed in a series of the unperturbed wave functions . The 
difficulty disappears if we consider the field 


H,=27H, sin (20 i) ; 


v. €. a field vanishing on the surface of the metal. The perturbation method 
gives then for s, a formula very similar to (10a); for S, and y we get 
again equations (16) and (18). 


§3. THE Two-DIMENSIONAL FIELD. 
The calculations are similar to those for the one-dimensional field. 
Therefore, it will be sufficient if we just indicate the main points and give 


the final results. ) 
Since we wish to arrive again at a relation of the form soc A, we must 


impose on A the condition 
div A=0. 


*In the case %)<max OF %>7)max the evaluation of the sum in (11) would 
have lead to a final result different from (16) ; but these cases are of no physical 
importance, and, therefore, we do not discuss them here. 

+ The same is true also in the case of a weak homogeneous field ; compare 
Boe aes that the first term in the bracket of (10 a) by itself would 
lead to a diamagnetic current as strong as that which is necessary for the 


explanation of the Meissner effect. 
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We find, then, that the vector-potential of the field (2) is 


A,=—4H,A cos (2 5) sin (2- *) ‘ 


A y=tHoa sin (2- i) cos (2= ‘) 5 
ASU: 


We shall consider the case of an infinitely large piece of metal, in which 
the unperturbed wave functions will be given by (3). The perturbation 
method leads, then, to the following expression for the perturbed wave 
functions : 


Pe n He, n,& 


ieH,r* 
Sho {(é- n) (aa Ves ne, n+, an ean rae —— No Pe No, 1) — nat) 


1 
(6-1-7) (<a: ES wag i Z jmcceell :)} Sep 


The 
where again n,=a/A. This formula shows that we shall have degenerate 
states when €+7n=-+7. These degenerate states, for which a separate 
calculation would be necessary, cannot be avoided by an assumption like 
(7a). But again we could avoid them if we replaced the field (2) by a 
field of the form 


H,=27H, cos (276, ;) . COS (201 t) 


with e. g. Es—yy=1; if &, n> 1, this field is physically equivalent to the 
field (2). The degenerate states occurring in the case €)=ny=n) have 
only the effect of making the calculation of the total current S (by 
summation of s over the occupied states) essentially more complicated. 
But the final results are exactly the same in both cases, 7. e. with and 
without degenerate states. Hence, it will be sufficient if we consider 
here the non-degenerate states described by (20). 

The current s associated with the state described by (20) is found to be 


Re eh ted es nab 
s= waras| +31 eae Ewa t ena aah | dy ou Wha) 


The calculation of the total current S from s can be made in much the 
same way as in the one-dimensional case. The final result is 


4a = ePA cs 
el es a Fag V(2mE 0) [1-5 (=) |.  eeral ea) 


differing from the corresponding one-dimensional formula (16) by a 


factor 2. If we calculate the equivalent susceptibility from (17), we find 
again the Landau value (18). 


(19) 
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: As in the one-dimensional case, here also the small value of S or x 
is due to a cancellation of much larger contributions of the individual 
states, again because of the spherical distribution of the electrons in 
momentum space. If we disregard the (relatively few) states for which 
&, 1? Sz, we see from (21) that the states for which either (E€—n)? <n? or 
(€+n)*<n@ are strongly paramagnetic; those are the states which in 
momentum space occupy positions inside the two layers of thickness )/2 
along the planes €--y=—0. The states situated outside those two layers 
are strongly diamagnetic. 

Instead of the metal extending infinitely in all three dimensions we 
could consider a prismatic piece of the metal, e.g. one filling the space 
0<x, y<a. When we try to solve the perturbation problem, we see 
again, as in the one-dimensional case, that this cannot be done with the 
field (2) ; instead we must take the field 


H,=27H, sin (2 i) sin (20 ‘) ‘ 
A A 
which vanishes on the surface of the metal. The calculation leads, then, 
to a formula for s which is essentially the same as equation (21). We 
notice for the subsequent discussion that the paramagnetic states have 
momentum vectors whose projections on the €, y-plane are directed at 
45° to the normal on the surface of the metal. 


§4. Discussion. THE MEISSNER EFFECT. 


Equation (11) can be written in the form 


eA P 
S= — mea eff > . . . > . ° ° ° . . ° (23) 

: oe i 1 ) x 
i = SS pS | —— SS : ee ee CA 
is Steal ee No lane 2n+N co 


A similar relation holds in the two-dimensional case. It has been stressed . 
by London a long time ago that the Meissner effect demands for Nog a 
value of Nog N. 

- The results of the previous calculations suggest a very simple way for 
increasing N,g to a value of Noga N. This will be the case if the effect 
giving rise to superconductivity produces a certain deformation of the 
otherwise spherically symmetric distribution of the electrons in 
momentum space, without influencing appreciably the wave functions 
and the currents s of the individual states ; 7. e. if we assume that in the 
superconductor the distribution of the occupied states in momentum 
space is not exactly spherical, but that there are certain misuse es 
on top of the Fermi sphere. The alternative independent explanation 
would be by assuming that in the superconductor the distribution of the 
electrons remains essentially spherical, but the wave functions and the 
currents s of the electron states are changed as much as is needed in 
order to account for the necessary increase of N,~- Though both these 
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factors most probably contribute to the actual increase of Nog, we shall 
restrict ourselves to the discussion of the first one, for the simple reason 
that even a qualitative discussion of the second factor would be extremely 
difficult. After all, it is almost certain that the first will be present in any 
case ; and its effects will, most probably, be there even in the presence 
of the second factor. 

By closer examination one sees that it is not possible to explain the 
Meissner effect by simply assuming that the prominences will be oriented 
in such a way as to contain strongly diamagnetic states. The reason is 
obvious. There are other orientations in which the prominences would 
contain paramagnetic states. And since diamagnetism corresponds to 
an increase of energy, while paramagnetism leads to a decrease of energy, 
the prominences will actually take the orientations corresponding to the 
strongest possible paramagnetism. 

If the prominences are not too small, this tendency towards 
paramagnetism may lead to the formation of a self-consistent internal 
magnetic field. By this we mean a field H either of the form (1) or (2), 
which gives rise to a current § so large that, in its turn, this current 
produces just the field H. It is possible to determine the order of 
magnitude of the size of the prominences required for this purpose. 
From the general relation 


AA= —47S8 


we find eff — PE (Vas): 
On the other hand, it follows from (10 a) or (21), if we denote by N’ the 
number of electrons contained in the prominences*, 


eff ~~ 


A 2 nar? 
pe erecta Mca 
0 
Combining the last two relations, and taking in account that &,,,.2N, 
we find 
Res, 2/3 e 
Tv ~w tn Foe EE Se ST eee en) eS (25) 
with n=N/V, n’=N’/V. Introducing the values A~10-5 em. and 
n ~~ 10" cm~? we get from (25) n’ ~10-4n. 
It seems to-day certain that the effect leading to superconductivity 
is the interaction of the metal-electrons with the thermal vibrations of 
the lattice, as suggested recently by Fréhlich (1950). It is significant 
that in Fréhlich’s theory prominences of the size required above are 


* If there were one prominence only, we should have to consider the last 
term in the expression of s given by (10): a state with one prominence only 
is associated with a permanent current independent of the magnetic field 
We avoid such a state—which, according to the general theorem of Bloch, 
cannot be the state of lowest energy—by assuming that prominences occur 
in pairs directed in opposite directions. 


ee 
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energetically possible. We, therefore, introduce the following postulate : 
The superconducting state is associated with an internal magnetic field of 
the form (2)*. In view of the very peculiar magnetic properties of the 
superconducting state the assumption that this state has an intrinsic 
magnetic structure is quite plausible. The only general restriction 
which we might, at present, impose on this structure is that it must be 
of a microscopic nature, since in the absence of an external field the 
superconductor shows no macroscopic field at all. This condition is 
fulfilled by the field (2). 

The explanation of the Meissner effect is now straightforward. 
According. to the remark at the end of §3 the prominent directions 
will be at right angles to the Oz-axis and at 45° to the normal on the 
surface of the metal. As regards the external magnetic field, which 
declines exponentially inside the superconductor, it is to be expected that 
the prominences will act on it in the same way as on the one-dimensional 
field (1) (with the axis Oy along the normal to the surface). Since Ny <Nmax? 
the states contained in the prominences will have |n|>m,/2 and, therefore, 
we shall have diamagnetism for the external field. Of course, this will 
be true only for not too strong external fields H’: Because of the 
diamagnetism of the superconductor as regards external fields, the field 
H’ will increase the energies of the electron states inside the prominences. 
Hence, it will tend to decrease the size of the prominences ; and when H’ 
reaches a certain critical value, the prominences will disappear, i. e. the 
superconducting state will be destroyed. 

The present considerations leave the value H, of the internal field 
entirely undetermined ; we might only guess that this value should be 
of the same order of magnitude as the critical magnetic field of the 
superconductor. For a theoretical determination of H, it will be necessary 
to pursue the calculations to a higher approximation. 
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X. CORRESPONDENCE. 


Calculation of the Density of States Curve for the 3d Electrons 
in Nickel. 


By G. C. Frercuert and E. P. WoHLFARTH, 
Department of Mathematics, Imperial College, London f. 


[Received November 23, 1950.] 


Irv is the purpose of this letter to present some preliminary results of 
calculations of the density of states, N(E), for the 3d electrons in a face- 
centred cubic lattice, using the approximation of tight binding. The 
calculations. are applied in the main to nickel. The final results have 
an important bearing on the physical properties of transition metals, 
such as magnetic, electrical and thermal properties. Comparison may 
be made with the results of analogous calculations by Krutter (1935) and 
Slater (1936), who employed the cellular approximation. 

The details of the tight binding approximation are well known (see, 
for example, Mott and Jones 1936). The atomic wave functions which 
have been used, and which correspond to the fivefold degenerate 3d level 
of the isolated atom, are 


di=a,0yf(r)/1?, dPo=agyef(r)/r, d3=agzaf(r)/r?, 
Py=a(V—Y)f(r)/?,  bs=As(822—r?) f(r)/r?, (1) 
where f(7) is the normalized radial function of the isolated atom and the 
a, are normalizing constants for the angular functions. Using (1) Bloch 


type functions may be constructed and the energy calculated from the 
determinantal equation 


| i ene LO eB ste 


The elements H,,,, may be expressed in the general form 
3 
Hin= 2 Anntimn(k) ee es (8) 


Here the 7, are functions of the three Cartesian components of the 
reduced wave vector k, expressible in terms of circular functions of 
3ak,, 3ak,, sak,, where a is the lattice constant. The coefficients AQ, 
have the general form 


AQ = [6#@—R)AV}(e—R')dr, Meese he ee) 


+ Now at Department of Mathematics, University College, Exeter. 
{ Communicated by Professor H. Jones. 
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where R, R’ refer to the position of a central atom and one of its nearest 
neighbours, and where AV is the difference in the potential due to an 
atom when isolated and when in interaction with its nearest neighbours 
in the metal. Owing to the symmetry properties of the angular functions 
in (1) some of the coefficients vanish and the others are interdependent, 
such that only six independent coefficients had to be calculated. For 
this purpose the most suitable available values of f(r) in (1) and of V(r), the 
potential due to an isolated atom, were taken to be those of Hartree 
and Hartree (1936) for Cu+ with exchange. These could be approximated 
to closely by simple algebraical expressions, involving exponentials. 
The atomic polyhedron was taken to be spherical and the assumption 
made that 4V—0 inside the central polyhedron and that outside it 
AV=V(|r—R’|). The integrations in (4) were carried out using prolate 
spheroidal coordinates. The resulting algebraical expressions, for the 
coefficients were evaluated using the appropriate value of the interatomic 
distance. The values of the six independent coefficients range from 
0-033 eV. to 0-260 eV., and these values, as that for the total width of 
the 3d band, are of such an order of magnitude as to justify the use of the 
tight binding approximation in this case. 

With the H,,,, known at all points in the Brillouin zone the energy 
may be calculated by solving (2). The functions 7, are, however, 
quite complicated, and hence it was decided, in the jirst place, to neglect 
arbitrarily certain of the non-diagonal elements such that the quintic 
equation (2) breaks into a cubic (arising from the atomic functions 
¢;, ¢2 and ¢,) and a quadratic (arising from ¢, and ¢5). pe 
approximation also underlies the calculations of Jones and Mott (1937) 
on the body-centred cubic lattice and is implicit in much of the work of 
Brooks (1940) on ferromagnetic anisotropy. Calculations are now 
proceeding which aim at. solving the complete quintiv in the energy 
range of interest for nickel, near the top of the band. Certain results 
may, however, be obtained at once, and, in addition, such comparison 
as may be made gives very close agreement with experiment using the 
curve which has been calculated so far. It is believed that although 
inclusion of the remaining H,,,,, may distort the energy surfaces derived 
from the two independent equations, the effect on the final N(E) curve is 
less serious. The highest energies are found, in any case, to correspond 
to regions of the zone where the neglected E,,,, are small. Finally the 
energy values which have already been computed may, in many cases, 
be used as first approximations to the roots of the complete quintic. 

The two resulting determinantal equations were solved in the jae 
ky, ky, k,>0 at all points of a lattice formed by taking ee . i 
7/6a, Cardan’s method being used to solve some 500 cubis. oe yi 
cellular approximation it was possible to calculate the siete a ok 
very much fewer points in the zone). Constant energy pas ie 
constructed in the planes k,=v7/6a(v=0... 12) for each a chet bi 
solutions. Finally N(E) curves were calculated by graphical integration. 
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A check on this stage of the calculation is provided by the fact that the 
area under each of the five curves corresponds to one electronic state 
per atom. The resulting provisional N(E) curve, the sum of the five 
constituent curves, is shown in the figure. 

The following deductions of particular interest may be made : 

(1) The total width of the 3d band is 2-70 eV. This result is independent 
of the neglect of the non-diagonal elements mentioned. The value is 
about half that given by Slater (1936), with whose N(E) curve the present 
curve may be compared. 

(2) The value of the degeneracy temperature for nickel calculated for 
0-3 vacant states per atom, corresponding to the right-hand vertical line 
in the figure, is 2-04 x 103 deg. K., in very close agreement with the value 
deduced by Wohlfarth (1949) from the magnetic and thermal properties. . 


Ferromage. Ni 


0 ! E (electron volts) e 3 


Provisional density of states curve for nickel. 


EL, electronic energy in eV., measured from zero at the bottom of the band. 
N(E), density of states, such that the total area under the curve corresponds 
to five states per atom of either spin, Fermi limits for paramagnetic, 
ferromagnetic nickel as shown. 


(3) The linear low temperature electronic heat coefficient, y, for 
ferromagnetic nickel (left-hand vertical line) is deduced to be 
1-15 10~calmol.-1deg-2. The observed value is 1:74 10-8, but as 
this includes a contribution due to the conduction electrons, the significant 
value for comparison is somewhat lower, possibly 1-5 10-3, 

(4) A completely general result concerning the degeneracy character 
at the centre of the zone, k=0, states that here the five energy levels 
split. into a threefold and a twofold degenerate level, separated, in general, 


———EEOEOEOEEeEeEeEEeEeEeaEeEeEe 
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by a finite gap. This result is completely verified in the present case, 
the gap being found to be 0-88 eV., again independently of the neglect of 
the non-diagonal elements mentioned. In calculations using the cellular 
method, on the other hand, the degeneracy at k=0 is found to be almost 
invariably fivefold—a result of the approximations used in the method. 

It is hoped to publish fuller details when the effect of the approximation 
indicated above has been more completely investigated. The necessary 
calculations are, however, so troublesome that it was felt that a 
preliminary account of work already completed may be of interest. 

We are indebted to Professor H. Jones for his most helpful advice, 
and one of us (G.C.F.) to London University for a postgraduate 
studentship. 
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Stars produced in Nuclear Emulsions by 150 MeV. Neutrons. 


By E. W. Tirrerton *, 
Atomic Energy Research Establishment, Harwell +. 


[Received November 6, 1950.] 


GARDNER AND PETERSON (1949) and Gardner (1949) have used the 
Berkeley cyclotron to carry out an investigation of stars induced in 
nuclear emulsions under deuteron and «-particle bombardment. Several 
authors, notably Harding e¢ al. (1948) and Page (1950) have investigated 
cosmic ray induced stars and recerttly there have been several papers 
describing stars induced by 7~ meson capture both for cosmic-ray mesons 
(Menon et al. 1950) and artificially created mesons (Chester and Goldfarb, 
Adelman and Jones 1950). ; 

In the present experiments a fast neutron beam, obtained from the 
Harwell Cyclotron by bombarding a beryllium target with 170 MeV. 
protons, has been used to irradiate C,, G; and NT, plates and the neutron 
induced stars have been studied. The neutron spectrum (Cassels et al. 
1950) is indicated in fig. 1 and gives a maximum yield at 150 MeV. The 
neutron flux was determined by integration and includes neutrons from 
60 MeV. to the maximum of the spectrum. At the plate location the 


flux was 2 x 105 n/cm.?/sec. 


* Communicated by the Author. ; nating 
+ Now on the staff of the Australian National University, Canberra, 
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Fast neutron spectrum from Beryllium target bombarded by 170 MeV protons. 
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Percentage of stars having N prongs plotted against prong number. 
Full circles : Present neutron induced stars. 
Open circles : Deuteron stars (Berkeley). 
Triangles : «-particle stars (Berkeley). 
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Five thousand stars have been examined to date and preliminary results 
were reported at the Harwell Conference (Titterton 1950). The stars can 
be classified into three main groups :— 


1. Those having a single energetic proton in addition to other low 
energy prongs. The proton is apparently knocked-on by the 
incoming neutron and the remaining excited nucleus breaks up to 
give further fragments. The spectrum of the knocked-on protons. 
has been examined and extends to 55 MeV. ; 

2. Stars in which the prongs are collimated in the forward direction 
but contain low energy tracks. These stars are presumably caused. 
by disintegration of the light elements 0, C, N present in the emulsion. 


Fig. 3. 
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Angular distribution of star prongs measured in the plane of the emulsion. 
Full circles : Neutron induced stars. 
Open circles : Deuteron induced stars. 
Triangles : «-particle induced stars. 


3. Stars in which the prongs are distributed more or less isotropically 
about the incoming neutron direction. These probably result from 
the Ag and Br atoms in the emulsion and are true evaporation stars. 
Evaporation stars are well known in plates exposed to the cosmic 
radiation and have been treated theoretically by Le Couteis (1950). 

The prong spectrum of the stars is plotted in fig. 2 racks it os sei 

with spectra obtained in Berkeley for deuteron and oil dee = peas 
stars. The distributions are similar except that more two-prong stars are 
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found in the present experiments. This difference is probably due to the 
fact that the emulsions employed were more sensitive than those used 
in the Berkeley experiments where protons of energy greater than 17 MeV. 
were not recorded. The angular distribution of the star prongs about 
the beam direction in the focal plane of the objective is shown in fig. 3 
where the fraction of prongs falling in 60° sectors is given as a function 
of the angular position of the sector. Comparison with Berkeley results 
shows good agreement with the deuteron distribution which might be 
expected, since in this case, many of the stars are probably formed by 
the capture of a proton or neutron in a stripping process. The greater 
predominance of tracks in the forward direction in the «-particle 
bombardment is compatible with the greater momentum transfer in this 
case. 

As the relative contributions to the star population by the heavy and 
light constituents of the emulsion are unknown the cross-section for star 
production is averaged over all the elements in the emulsion, excluding 
hydrogen which cannot give rise to a multiparticle star. The value 
obtained is 54 mB which is considerably smaller than the 210 mB found 
for 170 MeV. «-particles by Gardner (1949): 

If it is assumed, as seems reasonable, that in this energy range the 
cross-sections are in the ratio of the geometric cross-sections then the 
cross-section for star production in the light elements (C, N and QO) is 
25 mB and in the heavy elements (Ag and Br) is 95 mB. Because of. 
uncertainties in the measurement of the neutron flux these values are 
not accurate to better than +30 per cent. 


The plates were searched by Miss E. M. Chellingworth, Miss M. Hart, 
Miss K. B. Treacy and Mrs. B. Naish. The help of Mr. M. M. Cheshire, 
of the Cyclotron Group, in operating the machine is gratefully acknowledged. 
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Hammer Tracks in Neutron and Proton Induced Stars. 


By E. W. Trrrerton f, 
Atomic Energy Research Establishment, Harwell f. 


[Received November 6, 1950.] 


In the course of the experiments described in the preceding letter, 21 of 
the 5000 neutron stars were observed to be characterized by “ hammer ” 
tracks. Of these stars 18 were found in C, and 3 inG, emulsions. Such 
events are common in emulsions exposed to the cosmic radiation 
(Franzinetti and Payne 1948) but were not reported by Gardner 
and Peterson (1949) or by Gardner (1949) in the work with deuterons 
and «-particles. The events are usually attributed to Li nuclei which, 
after coming to rest in the emulsion, undergo decay according to the 


scheme *Li > SBe*, ®Be* > 24He. An alternative possibility has been 
suggested by Alvarez (1950) following his discovery of the nucleus 8B 


which undergoes a similar decay scheme °B se ®Be*, SBe* > 24He. 

The hammer tracks observed in the present experiment, whether they 
are due to *Li or *B, almost certainly result from the break-up of light 
elements since the emission of such fragments from Ag or Br nuclei at 
these excitation energies would be an improbable process due to the 
impenetrability of the Coulomb barrier. The prong spectrum and 
angular distribution for the 21 stars are compared with the equivalent 
distributions obtained in the case of neutron induced stars without 
hammer tracks in fig. 1 (a) and fig. 1 (6). The comparison is compatible 
with the view that the events are due to the disintegration of the light 
elements in the emulsion. 

Accepting this view application of the principles of charge and, to a 
lesser degree, momentum conservation, then allows conclusions to be drawn 
about the reactions responsible for these stars. Thus of the three-particle 
stars the majority (7) appear to result from the reaction 


C124+—Li8+Het+H!—22:7MeV. . . . . .« (1) 
and three from the reaction 
016 +n=Li8+ Be?+H!+n—49 MeV. A oP ES: 


the proton being knocked-on with high velocity as indicated in the 
micrograph of fig. 2. 

Of the seven four-particle stars only one can be due to carbon as all 
others have, in addition to the hammer track, at least one other prong 
of charge two. Four appear to be described by the reaction 

Ni44=Li8+Het+2H!+n—35 MeV. . .. . (3) 


Communicated by the Author. inbd 
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Fig. 1(a). 
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No. of prongs. 
Per cent of stars with N prongs plotted agaist the prong number N. 
Histogram : Stars with hammer tracks. 
Circles : All other stars. 


Fig. 1(b) 
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Angle from beam direction. 
Per cent of prongs per 60° sector as a function of angle between, centre of 
sector and the beam direction. 
Open circles : Stars with hammer tracks. 
Closed circles : Stars without hammer tracks, 
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and two by the reaction 


0+ n—Li®+He!+Het+H!—30 MeV. ... . (4) 
The five-prong star can only be due to oxygen and may be described by 
O%+n—Li8+Het+3H!+2n—58MeV. . . . . (5) 


a the 21 ous only the two-particle stars could reasonably involve 
the nucleus B. One of these comprises the hammer track plus a fast 
proton and on charge conservation grounds is therefore 


C™+n'—B*+H1+4n1—54MeV. . . . . . (6) 


the other shows a gross momentum unbalance about the incoming neutron 
direction and has a second track of Zs2. It could be either Li’ or B8. 


Fig. 2. 
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Micrograph of star assigned to reaction (2). 


These results suggest that Li® is a much more likely residual nucleus 
than B8 and, accepting the above assignments, the cross-section for Li® 
production, averaged over C, N and O atoms in emulsion is found to be 
0-6-.0-2mB. The energy spectrum of the Li® nuclei is indicated in 
fig. 3 (a) and the energy release in the break-up of the resulting Be*® 
nuclei is given in fig. 3 (b). The latter histogram is in agreement, within 
the statistics, with that found for Be® resulting from the reaction Be®ypLi® 
(Titterton 1950). ; 

Experiments are now in progress to study stars produced by 170 MeV. 
protons and the preliminary results show a similar percentage of stars 
with hammer tracks. This agrees with Wright (1950) who finds that *Li 
is formed when a variety of elements are bombarded with 340 MeV. 
protons and also with 190 MeV. deuterons. Furthermore, in a study of 
3000 7- meson induced stars Adelman and Jones found 11 cases of hammer 
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track emission. This frequency compares well with that observed in 
present experiments with neutrons and protons as might be pa 
since the nuclear excitation in this case would be about the same value 


Fig. 3(a@) 
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Energy spectrum of ®Li nuclei emitted from the stars. 
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Energy of «-particle pair MeV. 
Energy release in the break-up of the 8Be nucleus forming the 
head of the hammer. 
Full lines, present 21 events. 
Dotted lines, events resulting from the reaction Be® yp Li8. 


(142 MeV.). Combination of these results suggests that 8Li should be 
found as a residual nucleus in the break-up of highly excited light nuclei 
and it should therefore have been found in the high energy «-particle 
bombardments. 
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The mechanism proposed by Wright (1950) to explain Li® production 
from light elements requires the incoming particle to be captured and the 
compound nucleus to boil off the necessary number of neutrons and 
protons leaving a residue consisting of ®Li. This is not in general 
agreement with the present observations. As shown above an excited 
nucleus is formed by inelastic scattering or by capture of the neutron 
and then breaks up into charged particies of which a large proportion 
are «-particles. In general very few neutrons are emitted. 


Thanks are due to Miss M. Hopkins for microscope work and Mr. T. A. 
Brinkley for the photomicrograph. 
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XI. Notices of New Books and Periodicals received. 


Progress in Biophysics. Edited by J. A. V. BuTLER and J. J. RANDALL, F.R.S. 
(Butterworths.) Price 50s. ; 

i i i to define, but this well- 
Tur boundaries of biophysics have never been easy 5 
produced volume helps both to indicate the fringe of recent ne, one 
consolidate the ground already won, As the editors point out, ie ae ry 
and physiology have their own domains and mere application of a p Se: (o) 
chemical technique to a biological problem does not necessarily constitute a 


passport to this mandated territory. Perhaps the main justification for a 


ubdivision of science such as biophysics is that it provides opportunities for 
choad fertiliea ttl which may Share offspring more bold and vigorous than 
nts. 
pra soa the volume, with separate chapters by a number of experts, 
there is noticeable a healthy tendency to stress what may be called the dynamic 
spects of the various subjects. s 
i There is a good author ides but the subject index is rather brief for a volume 
covering so wide a field. For most readers too, it would be more aie ee if 
the separate bibliographies were arranged in alphabetical order. It would also 
be helpful if in future editions all the contributors could be induced to provide 
summaries of their material and personal views. jen 
Dr. Butler and Professor Randall have performed a valuable service. in pro- 
ducing a work in which there is so much information, some wisdom and no little 


wit. W. G-W. 


Cosmical Electrodynamics. By H. Aurvin. (Oxford.) Price 25s. 

TuE subject of which this book treats is a new one. It deals with the electro- 
magnetic behaviour of an ionized compressible gas in hydrodynamic motion, 
and with the applications of this study to problems in astrophysics. 

At present we are very far from having a complete solution of the 
electrodynamic-hydrodynamic problem in its most general form. Unfortun- 
ately a proper investigation of the astrophysical effects requires nothing short 
of such a complete solution. So we are forced to makeshift methods. These 
have mostly followed one or other of two lines. There is the cautious school 
of thought that accepts only those electromagnetic effects that seem unavoidable. 
The second school of thought, on the other hand, accepts every electromagnetic 
effect whose existence cannot be disproved. For example, this school will 
regard an electric field as being present in a particular problem so long as its 
existence cannot be disproved, whereas the first school will not accept the 
existence of the electric field unless its absence can be disproved. Dr. Alfvén 
belongs to the second school of thought, and to the reader of his book this 
carries an important advantage. For Dr. Alfvén’s applications of electro- 
magnetism to astrophysics are so extensive that many of the ideas set forth 
in his book are reasonably certain to have permanent value. Which these are, 
however, the reader must discover by his own efforts. 

One of the most important current issues concerns the relation between the 
electric current density in an extremely diffuse ionized gas, the electromagnetic 
field, and the motion of the gas. This affects not only the structure of the 
solar atmosphere, but also whether or not there can be a galactic magnetic field 
of appreciable magnitude. Alfvén quotes a formula due to Cowling without 
apparently noticing that Cowling’s result depends on the assumption that the 


gas is not being appreciably accelerated. This places an important restriction — 


on the use of Cowling’s formula which may be stated briefly as follows :— 
To an observer moving with a particular element of gas the electric and 
magnetic fields at the position of the element are substantially parallel. 
Without this result it seems impossible to attempt, as Alfvén does, the 
investigation of very diffuse gases. 
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Another question that is now open to doubt is whether under normal 
conditions the general solar magnetic field can extend much above the 
photosphere. Certainly the observed variations of angular velocity from one 
part of the Sun’s atmosphere to another are inconsistent with the simple dipole 
field assumed by Alfvén. This has serious repercussions on the sections dealing 
with the structure of the chromosphere and corona. 

Tn reading this book one can hardly fail to notice the author’s tendency to 
claim a displacement of a problem as its solution. For instance, it would not 
be reckoned a complete answer to a question such as: ‘‘ Why did you arrive 
at work so early this morning ? ”’ to say : “‘ Because I left home an hour before 
my usual time.” Yet this is what Alfvén does on numerous occasions ; for 
example, in the problem of the reversal of the magnetic polarities of sunspots, 
in that of the darkening of sunspots, and in postulating the electric field 
necessary to accelerate the auroral particles. 

But in spite of all criticism, and the author can hardly expect some of his 
ideas to escape criticism, this book cannot fail to be most useful to all who are 
interested in the problems of solar physics and of cosmic rays. It should be 
on the book-shelves of every astrophysicist. F. Hoye. 


A Hundred Years of Physics. By Witu1am Wison. [Pp. 519.] (Duckworth 
and Co., London, 1950.) Price 21s. 

THoszE who hope to find in this volume a history of the revolutionary changes in 
physical ideas which have taken place during the period covered (1850-1950) 
may be somewhat disappointed. The book is rather a kind of supplementary 
text, suitable for revision purposes, in which rather more than the usual 
emphasis is laid on the chronology of the various topics considered. An honours 
graduate in Physics—actual or potential—could read it with pleasure and profit ; 
but the author’s claim in his preface that the reader needs very little previous 
knowledge of physics and very slight mathematical equipment must be taken 
cum grano salis. For the rest, the book is well written, well produced, with 
many references to the original publications and a satisfactory index. N.T. 


Contributions to Mathematical Statistics. By R. A. Fisozr. Wiley (Chapman 
and Hall.) Price 60s. net. 

Tuis volume is a collection of reprints of 43 papers by an eminent statistician 
selected by him as the more important of his published papers. They cover a 
wide variety of subjects and are accompanied in many cases by author’s notes 
added for the purpose of the volume. Though photographically reproduced 
from the original journals, the volume is attractively prepared on good paper. 
Many will find in it a useful work of reference and study. 


_ Methods of Mathematical Physics. By H. and B.S. Jerrreys. Second Edition. 

[Pp. 708.] (Cambridge University Press.) Price £4 4s. Od. net. 
THE value of this important work is shown by the appearance of a second 
edition only four years after the original. The first edition was reviewed in 
this journal (Vol. 38, p. 226, 1947). It has now been considerably revised, 
in some cases by the addition of new material or the expansion of old, as in 
the account of matrices and the treatment of relaxation methods. In other 
cases the subject-matter has been mainly re-written, as in the chapter on 
multiple integrals. Where possible the proofs have been either replaced by 
shorter ones or generalized. , 

An extraordinarily wide range of topics is covered, as indeed must be the case 
in achieving the authors’ aim to give an account of those parts of pure 
mathematics that are most frequently needed in physics. 


[The Editors do not hold themselves responsible for the wews 
expressed by their correspondents. | 


NOTICE TO CONTRIBUTORS 


UNIT OF HEAT 


The Council of the Royal Society, having considered a 
memorandum entitled “The Unit of Heat’’, recommends that 
quantities of heat and all other dependent concepts such as 
specific heat, latent heat, entropy etc., should in future be 
expressed both in joules and in calories. 

The Editors express the hope that, in submitting papers for 
publication in the Philosophical Magazine, contributors will also 
adopt this precedure. 


i: ERRATA, 


An Interferometric Study of some Optical Properties of Evaporated Silver 
Films, by R. C. Faust, Phil. Mag., 1950, 41, 1238. 


P.1241, line 5, read ‘‘ glass substrate ”’ for “ glass substitute”. 
P.1243, line 30, read “rate of evaporation”? for “rate of absorption ”. 


P.1254, Insert the reference 
BrossgEL, J., 1947, Proc. Phys. Soc. Lond., 59, 224. 


